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FLOWFIELD MEASUREMENTS IN A SEPARATED AND REATTACHED 
FLAT PLATE TURBULENT BOUNDARY LAYER 


William P. Patrick 


SUMMARY 


An experimental investigation of turbulent boundary layer separation with 
reattachment on a flat plate test surface was conducted in the United 
Technologies Research Center Boundary Layer Wind Tunnel. The objective of the 
study was to provide a comprehensive set of data for use in assessing the 
accuracy and assisting the further development of computational procedures for 
predicting boundary layer separation phenomena. The approach was to conduct a 
large scale, two-dimensional flow simulation of an airfoil turbulent separa- 
tion bubble at low subsonic Mach number by imposing a strongly adverse and 
subsequent favorable pressure gradient. A combination of laser Doppler veloc- 
imetry, hot-wire anemoraetry, pneumatic probing techniques, and flow visualiza- 
tion were used as diagnostics. 

The boundary layer separated and reattached over an axial distance of 
approximately 55 cm. The maximum bubble thickness, measured to the height of 
the mean dividing streamline, was 17 cm, which was twice the thickness of the 
inlet boundary layer. A strong backflow region in which the flow moved 
upstream 100 percent of the time was measured near the test' surface over the 
central 35 percent of the axial extent of the bubble. A laminar backflow 
boundary layer having pseudo-turbulent characteristics including a log-linear 
velocity profile was generated under the highly turbulent backflow. Velocity 
profile shapes in the reversed flow region only matched Simpson's universal 
backflow profile at the upstream edge of the separation region where the back- 
flow was intermittent. A smoke flow visualization movie and hot-film measure- 
ments revealed low frequency nonperiodic flapping at reattachment. However, 
forward flow fraction data at reattachment and mean velocity profiles in the 
redeveloping boundary layer downstream of reattachment correlated with 
backward-facing step data when the axial dimension was scaled by the distance 
from the maximum bubble thickness to reattachment. 

Other principal observations concerning the separation bubble flowfield 
were that an outer inviscid rotational flow can be defined which essentially 
convects over the blockage associated with the inner, viscously dominated 
bubble recirculation region. In the outer flow region, velocity profiles were 
similar over the axial extent of the bubble while significant profile altera- 
tion occurred within the inner flow region, including the development of the 
reversed flow boundary layer at the base of the recirculation zone. The 
identification of this three-tiered structure is believed to have important 
implications relative to the numerical modeling of the flowfield. 
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CHAPTER 1 


INTRODUCTION 


Development of reliable numerical procedures for predicting compressor 
and turbine airfoil boundary layer separation is important to the gas turbine 
design process since the potential for occurrence of separation is high, when 
it occurs its effects are usually large, and present empirically based predic- 
tion methods have significant shortcomings. 

The adverse pressure gradient responsible for boundary layer detachment 
can arise from local surface curvature effects or be imposed on the boundary 
layer by the outer flow (flow deceleration or incident shock waves). Depend- 
ing on the strength of the adverse pressure gradient and the boundary layer 
characteristics (whether laminar or turbulent and the boundary layer shape 
factor) separation can occur on either the airfoil pressure or suction surface 
and at various axial positions. The boundary layer may remain deattached or 
reattach after a short or long distance (as a percent of chord) and it may 
undergo transition to turbulence if initially laminar. 

A schematic of a turbomachinery blade showing several of the possible 
separations in terms of location, source and downstream development is given 
in figure 1-1. A closed leading edge bubble (termed a transitional bubble) is 
shown on the suction surface where typically deceleration of a laminar 
boundary layer causes separation, and transition to turbulence in the bubble 
shear layer fosters reattachment. Under a sufficiently adverse pressure gra- 
dient imposed by the outer flow, reattachment may not occur and the 
phenomenon of massive stall will result. A second form of suction surface 
separation encountered in the outboard regions of supersonic fan blades is 
shown further aft. This shock-induced phenomenon usually involves separation 
of a turbulent boundary layer and arises from interaction of the boundary 
layer with a bow shock from a neighboring blade. Depending primarily on shock 
strength and the degree of downstream diffusion, reattachment may or may not 
occur. Another region of surface separation is the airfoil trailing edge 
where the thick, blunt geometries required for structural and durability 
reasons invariably result in a surface curvature-induced separation. A fourth 
separation phenomenon is shown on the pressure side where surface curvature 
can induce a turbulent boundary layer ("cove") separation which subsequently 
reattaches under the influence of a favorable pressure gradient. Complicating 
these isolated phenomena is the possibility of interaction where, for example, 
a boundary layer weakened by negotiating the pressure rise of a shock may 
separate upstream of the trailing edge under the influence of diffusion and 
hence alter the trailing edge separation. 
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It is therefore clear that the potential for occurrence of separation is 
high and that once separated, the subsequent downstream development can have a 
significant impact on the aerodynamic performance (and for the turbine appli- 
cation, also the surface heat transfer). While the prediction of the incep- 
tion of separation is important, and in some cases reasonably well treated by 
existing methods, the ability to predict the subsequent downstream flow devel- 
opment is equally important. Existing airfoil design procedures typically 
rely upon an inviscid transonic cascade calculation to predict the overall 
flowfield and a coupled boundary layer calculation to account for near surface 
viscous effects. In the absence of separation the procedures are generally 
satisfactory for design purposes. When separation occurs, however, lack of a 
reliable calculation method requires the use of approximate methods derived 
from test data and past experience. This empirical approach, while successful 
in many applications, not only limits the designer to interpolation within the 
experience base but also presents problems in the application of the available 
codes. For example, in treating either a closed leading edge or shock-induced 
separation bubble, it may be possible to use correlations to define the 
approximate extent of the separation region. Uncertainty exists, however, in 
boundary layer characteristics aft of the bubble. Since these are required to 
initiate an attached boundary layer calculation over the remainder of the 
airfoil surface, the accuracy of the overall airfoil calculation procedure is 
impacted. 

The focus of the present study was the pressure side "cove" separation 
depicted in figure 1-1. The approach was to conduct a large-scale, two- 
dimensional flow, low Mach number simulation on a flat plate test surface. 

The flat plate geometry was selected since it provided the simplest test case 
for code assessment purposes. There have been no previous investigations of 
this turbulent separation bubble geometry using non-intrusive laser 
velocimetry although Simpson (refs. 1-4) has performed extensive research 
relative to flat plate turbulent boundary layer separation without downstream 
reattachment. Simpson* s work has provided significant information regarding 
both the mean and turbulent velocity fields within the separated flow region 
and some of his results apply directly to the leading edge region of the 
bubble investigated in this study. Although a number of other investigators 
have conducted closed separation bubble studies which have indirect relevance 
to the problem considered here, test geometries invariably involved rapid 
surface curvature and salient edges which produced fixed separation points and 
also affected the development of the reversed flow region within the bubble. 
Only Perry and Fairlie (ref. 5) and Cutler and Johnston (ref. 6) have reported 
investigations of flowfields comparable to the current study. Their studies, 
described in Chapter 2, used pitot and hot-wire traverses to obtain data 
within the separation bubble and, in addition. Cutler and Johnston used 
thermal tufts at the wall. 
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For the near future, computer codes providing numerical solutions of the 
Reynolds-averaged equations will continue to be the most useful engineering 
approach for the solution of flowfields having turbulent separation bubbles 
(ref. 7). Such codes use numerical correlations and zonal modeling to divide 
the flowfield into several regions, each dominated by a particular type of 
flow, and to analyze each region by the computationally optimum numerical 
technique for that region (ref. 8). Current codes use flow models which have 
been derived from open separation bubble data for the separation region and, 
primarily, from backward-f acing step data for the flow at and downstream from 
reattachment. Within the separation bubble Gosman (ref. 9) noted that insuf- 
ficient reliable data exists to define the structure of the backflow region. 
He suggested that future experiments should be on large-scale models using 
laser aneraoraetry or other quantitative directionally sensitive techniques. 

The present study was undertaken to provide such detailed data for the entire 
flowfield of a separated flow configuration relevant to gas turbine airfoil 
design. The overall objective was to provide a comprehensive data base from 
which the accuracy of existing models could be assessed and improved models 
could be developed. 


The author wishes to acknowledge helpful suggestions and comments from 
M. J. Werle (UTRC), R. W. Paterson (UTRC), A. J. Strazisar (NASA), J. E. 
Carter (UTRC) , M. F. Blair (UTRC), D. E. Edwards (UTRC), E. M. Greitzer (MIT) 
and the late W. D. NcNally (NASA). In addition, the author notes that the 
laser velocimeter error analysis in Appendix C was derived from a comprehen- 
sive study of LV system errors conducted as part of UTRC's independent 
research activities in gas dynamics. 
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Figure 1-1 Airfoil Surface Separation 
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CHAPTER 2 


REVIEW OF PREVIOUS INVESTIGATIONS 

Detailed measurements of turbulent boundary layer separation and 
reattachment on a rotating compressor or turbine blade are unavailable because 
of the obvious difficulties of obtaining high resolution measurements in an 
operating engine environment. As noted in the previous chapter, only two 
experimental studies of the more tractible problem of turbulent boundary layer 
separation and reattachment on a stationary flat plate at low subsonic speed 
have been reported in the literature (refs. 5 and 6). Each of these studies, 
which yielded partial data bases, will be described below. To date, flow 
modelers have had to rely on the extensive body of literature which exists for 
both reattaching and separating turbulent flows. 

Subsonic turbulent flow reattachment has been reviewed by Mueller (ref. 
10) , Bradshaw and Wong (ref. 11), and more recently by Eaton and Johnston 
(ref. 12). Eaton and Johnston noted that at the time of Bradshaw and Wong*s 
review in 1972 very few reliable data for reattaching flows were available. 
Although details of reattachment length and mean flow profile development 
downstream of reattachment were known with good accuracy, turbulent stress 
data required to develop turbulence models were sparse. Since that time the 
development of laser velocimetry and pulsed-wire anemometry has enabled 
several investigators (refs. 13-17) to obtain reliable, quantitative 
reattachraent data. Details of these investigations will not be reviewed in 
this section although pertinent studies such as those by Chandrsuda and 
Bradshaw (ref. 18), Driver and Seegmiller (ref. 14), and Westphal, et al . 

(ref. 17) will be discussed in Chapter 5, "Analysis and Discussion". 

Ota and his colleagues (refs. 19-21) have studied separation bubbles 
having turbulent reattachment which form at the leading edge of blunt flat 
plates and circular cylinders aligned with the flow. These flows have a 
negligible boundary layer thickness at separation and therefore are unaffected 
by upstream boundary layer history as can be the case for b ackward-f acing step 
flows. Ota, et al. studied the redeveloping boundary layer downstream of 
reattachment and found strong similarities between their data and the backward 
facing step data of Bradshaw and Wong (ref. 11), Mueller, et al. (ref. 22), 
and others. 

Between the pioneering efforts of Sandborn and Kline (ref. 23) in defin- 
ing turbulent separation criteria and the recent definitive experimental 
studies of Simpson and his colleagues (refs. 2-4), numerous studies of turbu- 
lent boundary layer separation have been reported in the literature. These 
studies have been reviewed by Simpson (refs. 8 and 24) and except for the 
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studies of Simpson and his colleagues (refs. 1-4) and a later study by Smith, 
et al. (ref. 25) will not be discussed in this section. However, as noted 
above for the turbulent reattachraent data, pertinent studies will be included 
in the data comparisons presented in Chapter 5, "'Analysis and Discussion". 

Little reliable quantitative data was obtained within turbulent separated 
boundary layers prior to 1974 due to a lack of proper instrumentation. 

Simpson, et als. (ref. 1) pioneering effort in applying a one-component 
directionally sensitive laser velocimeter to the study of separating turbulent 
airfoil-type boundary layers revealed several new features of these flows 
including the following: 

(1) Lawof-the-wall velocity profile scaling is valid up to intermittent 
s eparat ion. 

(2) Intermittent separation occurs when the wall static pressure grad- 
ient begins to diminish. 

(3) Normal stress terms in the momentum and the kinetic energy equations 
are important near separation. 

(4) The separated flowfield shows some similarity of the strearawise mean 
velocity, the streamwise turbulence component, and the forward flow 
fraction. 

Later studies by Simpson and his colleagues using a two-component LV 
system culminated in a series of three papers (refs. 2-4) which presented 
fundamental mean flow and turbulence structure data from which improved models 
of the backflow region could be developed. Several significant conclusions 
resulted from these studies including the following. 

(1) The backflow velocity profiles are not consistent with law-of-the- 
wall type scaling. Simpson (ref. 26) postulated a universal mean 
backflow velocity profile scaled on the maximum negative mean veloc- 
ity 1% and the distance, N, of the maximum negative mean velocity 
from the wall. 

(2) Within the backflow the turbulent velocities are comparable to the 
mean flow velocities although the turbulent shearing stresses are 
low. 

(3) Mixing length and eddy viscosity models are physically meaningless 
in the backflow and have reduced values in the outer region of the 
separated flow. Reynolds shearing stress must be modeled by 
relating it to the turbulence structure and not the mean flow 
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gradients because the mean backfow velocities are obtained by time- 
averaging large turbulent fluctuations and are not related to the 
source of the turbulence. 

(4) A three-tiered structure was apparent in the backflow. It consisted 
of a viscous layer closest to the wall, a wider flat intermediate 
layer, and an outer backflow region which was part of the large- 
scaled outer region flow. 

Later work by Simpson has been directed toward obtaining an understanding of 
unsteady separating turbulent boundary layers (ref. 27). 

Smith, Hastings, and Williams (ref. 25) studied the structure of an 
equilibrium separated boundary layer on a flat plate with a single component 
LV system. The boundary layer was first accelerated and then decelerated such 
that the boundary layer edge velocity varied with distance according to a 
negative power law (i.e., U = Cx* 111 ). Mean axial velocity and turbulence 
intensity profiles were measured at several axial stations and one pair of 
profiles of the mean transverse velocity and turbulence intensity were 
measured in the separated region. Static pressures were measured on the test 
surface centerline. Forward flow fraction and shear stress profiles were not 
measured. 

Equilibrium separated flow was maintained for 1.1 m before the boundary 
layer on the test surface and tunnel ceiling merged. Some flow asymmetry was 
noted and measurements indicated that the flow was not in strict two- 
dimensional equilibrium. However, boundary layer displacement and momentum 
thickness grew linearly as expected for equilibrium flow. 

Limited measurements of backflow velocity profiles were presented. Data 
scatter appeared to be considerable near the test surface although no estimate 
of measurement error was presented. Smith, et al. corroborated Simpson's 
finding that the backflow velocity profile could not be modeled using the 
universal 1 aw^o f-t he-wall. They also concluded that normal Reynolds stresses 
were considerably more important in the momentum integral equation than would 
be expected on the basis of equilibrium attached flow. 

In addition to the separated turbulent boundary layer studies of Simpson 
and his colleagues and Smith, et al., only two experimental studies of turbu- 
lent boundary layer separation with reattachment on a flat plate have been 
reported in the open literature. These studies were conducted as Ph.D. 
dissertations by B. D. Fairlie under A. E. Perry at the University of 
Melbourne in 1973 (ref. 5) and by A. D. Cutler under J. P. Johnston at 
Stanford University in 1984 (ref. 6). 
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Perry and Fairlie (ref. 5) separated an equilibrium turbulent boundary 
layer to generate a large separation bubble (41 cm long x 8 cm thick at the 
maximum height of the dividing streamline) on the floor of a low speed open- 
return wind tunnel. The test section had a rectangular cross section measur- 
ing 61 cm long x 11 cm high at the inlet. It was 2. 5 m long with a specially 
constructed flexible roof, allowing any desired streamwise pressure gradient 
to be obtained. The variable tunnel outlet height was set at approximately 60 
cm for the two test cases reported in ref. (5). A false roof mounted below 
the flexible roof was used to provide a starting point for a new roof boundary 
layer thereby inducing the flow to preferentially separate from the flat test 
surface. Smoke, tufts and surface flow visualization indicated moderate 
three-dimensionality in the flowfield. 

Mean velocity, flow angle, and surface static pressures were measured 
throughout the flowfield. Neither forward flow fraction, turbulence, nor 
shear stress data were reported. The surface static pressures obtained on 
the floor centerline showed a pressure plateau over a 37 cm long region under 
the bubble with pressure maxima at the separation and reattachment points. 
Velocity traverses were made normal to the plate on the duct centerline 
throughout the flowfield. Data was obtained with flattened pitot probes 
upstream of separation. Near separation and within the separation bubble a 
twin hot-wire probe was used to obtain mean-velocity and flow angle measure- 
ments. Flow angles were obtained by traversing a temperature sensitive wire 
through the wake of the velocity sensitive wire. As reported by Perry and 
Fairlie, this intrusive hot-wire probing technique is inaccurate in highly 
unsteady flows or regions having intermittent flow reversals due to signal 
rectification by the velocity sensitive wire. Reported velocity profiles 
showed discontinuities at the edge of the backflow region. The velocity 
profiles indicated strong backflow within the separation bubble. However, the 
strength of the backflow could not be quantified without forward flow fraction 
data . 


Cutler and Johnston (ref. 6) studied the effect of the state of equi- 
librium or disequilibrium of the inlet turbulent boundary layer on separation 
and reattachment on a flat plate. Tests were conducted in a low speed closed- 
return wind tunnel having a 30.48 cm wide rectangular tests section which 
diverged over aim length from an inlet height of 10.16 cm to an outlet 
height of approximately 40 cm. Boundary layer bleed slots on the upper 
diverging wall caused the flow to preferentially separate from the lower wall 
flat plate test surface. Measurements were made of wall static pressure 
distributions, forward flow fraction at the wall using a thermal tuft, total 
pressure profiles with pitot probes, and mean velocity and turbulence param- 
eters using normal and crossed hot-wire probes. 
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Complete detachment of the boundary layer was not achieved in the Cutler 
and Johnston test cases. The test surface boundary layer passed through 
incipient detachment, closely approached detachment, and then fully 
reattached. The minimum value of the forward flow fraction measured at the 
wall was 0.59. Careful conditioning of the flow ensured excellent flow two- 
dimensionality in the inlet flow. Nevertheless nonnegligible three- 
dimensional effects were evident near detachment and reattachment when the 
terms in the two-dimensional momentum integral equation were evaluated using 
exerimental data. 

Cutler and Johnston drew several significant conclusions from their study 
including the following: 

(1) A tripped boundary layer which had not recovered to an equilibriun 
profile separated further downstream than a boundary layer having an 
equilibrium profile at the same operating conditions. 

(2) Since mean streamline angles to the wall are not small near the 
separation region even though the wall is flat it becomes important 
to model flow quantities, particularly turbulent shear stress, in 
coordinates locally orthogonal to the streamline. 

(3) Mean streamline convex curvature becomes significant between the 
detachment and reattachment locations, even though the wall is flat. 
The effect of the streamline curvature is to reduce the shear stress 
and the entrainment of the freestream fluid into the boundary layer 
and to induce nonnegligible static pressure gradients normal to the 
flat wall. 

(4) To calculate flows downstream of a reattachment, the development of 
mean-velocity-profile quantities and of turbulence quantities must 
be modeled separately because there is no unique relationship 
between these quantities. 
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CHAPTER 3 


DESCRIPTION OF THE EXPERIMENT 

3.1 Experimental Arrangement 
3.1.1 Wind Tunnel Arrangement 

The experiment was conducted in the UTRC Boundary Layer Wind Tunnel which 
is shown schematically in figure 3-1. It is a low speed wind tunnel of recir- 
culating design and consists of a blower, settling chamber/ plenum, contraction 
nozzle, boundary layer bleed scoop, boundary layer test surface, downstream 
diffuser, coarse air filter, temperature control heat exchanger, and return 
duct. This tunnel has been designed to produce large scale two-dimensional 
boundary layers with Reynolds numbers, freestream turbulence levels and pres- 
sure coefficient distributions typical of turbomachinery blades. A thorough 
calibration of the flow characteristics of the tunnel (ref. 28) has shown the 
flow to be highly two-dimensional and steady. The freestream turbulence level 
can be varied between 0.23 percent and 10 percent with the addition of turbu- 
lence generating grids. 

For this test program several modifications were made to the tunnel as 
shown in figure 3-2. The major modification was the installation of a new 
test section, described below, which was designed to produce a large two- 
dimensional turbulent separation bubble. Another major modification was the 
addition of a subraicron filtration system downstream of the diffuser section 
to permit LV seeding in the test section without contamination of the fine 
mesh screening in the settling chamber. 

The 1.83 m x 1.83 m cross-section submicron filtration system consisted 
of nine filtration modules arranged in parallel in a 3 x 3 matrix. Each 
0.61 m x 0.61 m cross-section filtration module consisted of a prefilter, a 
bag filter, and an absolute HEPA filter. Tests of a single filter module 
indicated that the filtration system was extremely efficient. Under condi- 
tions of heavy LV seeding the air in the return duct was found to be two 
orders of magnitude cleaner than the ambient unseeded air circulating in the 
unmodified tunnel. The filtration system had a nominal 3 cm 1^0 pressure drop 
at the 283 SCMM (10,000 SCFM) tunnel operating condition used in the test 
program. 

Several other modifications were made downstream of the contraction noz- 
zle. The leading edge bleed scoop was covered by a contoured wall section 
(see inset in fig. 3-2) and the original boundary layer test surface was used 
as a boundary layer development section. The downstream horizontal variable 
angle diffuser was then replaced by the new diffusing test section. An 
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0.64 cm square bar was installed downstream of the contraction nozzle to pro- 
vide a two-dimensional boundary layer transition process. The axial location 
of the trip was 245 cm upstream of the test section inlet (i.e., x = 0). 

3.1.2 Test Section Arrangement 

The test section consisted of a constant area approach section which 
functioned as a continuation of the boundary layer development section, a 
diverging section to produce the large adverse pressure gradient required 
to separate the turbulent boundary layer from the test plate, a converging 
section to force its reattachment and a constant area exit duct. The fixed 
dimensions of the test section were: length, 152 cm (60 in.); width, 86.4 cm 

(34 in.); inlet height, 21.3 cm (8.38 in.); exit height, 15.2 cm (6 in). 

The upper wall test surface was fabricated from 1.27 cm thick aluminum 

tooling stock. The plate was machine ground and hand buffed with steel wool 

to ensure surface smoothness. The lower wall was fabricated from 1.27 cm 
thick aluminum plates joined by flexures. The joints between the flexures and 
the aluminum plates were fabricated to eliminate steps in the surface contour. 
The lower wall was moveable by means of jack screws. During shakedown testing 

the lower wall geometry was adjusted over a wide range until a large separa- 

tion bubble having a minimum amount of three-dimensionality and unsteadiness 
was achieved on the test surface. Photographs of the final test section con- 
figuration selected for the experiment are shown in figures 3-3 and 3-4. 

3 . 1 . 2 . 1 J*o\indar^_L .aye r__B le j^d_Jl£S tem 

Boundary layer bleed on the test section diverging wall was used to pre- 
vent the preferential separation of the turbulent boundary on this surface. 
Test surface corner bleed was also used to improve two-dimensionality of the 
test section flowfield. The bleed system was driven by natural suction caused 
by a static pressure differential between the test section and the inlet to 
the submicron filtration system. A 22.9 cm 1^0 pressure differential was 
created by installing two 40 percent open perforated plates in series just 
upstream of the turning vanes ahead of the vertical diffuser section to raise 
the static pressure in the test section. 

The diverging lower wall was equipped with two bleed scoops which spanned 
the width of the test section as shown in the photograph in figure 3-3. The 
scoops were designed to force the lower wall boundary layer to remain attached 
throughout the strong adverse pressure gradient field required to produce 
two-dimensional separation on the test surface. The first scoop, which had a 
fixed height of 2.87 cm, could be moved axially, parallel to the lower wall, 
over a 22.9 cm. range. The second scoop which could be moved axially and 
vertically as well as rotated was found to be unnecessary and remained closed 
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during the experimental program. Approximately 30 percent of the flow 
entering the test section was removed by the operational lower wall bleed 
scoop. 

The upper wall bleed system consisted of two .32 cm wide, 152 cm long, 
suction slots on either side of the test surface at the corner between the 
test plate and side wall. The purpose of these bleed slots was to improve the 
two-dimensionality of the flow near the test surface by removing low momentum 
fluid caused by the sidewall boundary layers. Each bleed slot was connected 
to a plenum located above the test section (fig. 3-5) by fifteen 2.54 cm dia. 
holes drilled through the test plate support structure. Partitions in the 
bleed slots and the plenums divided each suction slot into three axial seg- 
ments each having independent control of suction strength by means of valving. 
Bleed flow was ducted from the plenums to the filter inlet by a 10 cm dia PVC 
pipe. Less than one percent of the flow entering the test section was removed 
by the upper wall bleed system. Prior to installing the corner bleed slots, 
sidewall fences were used in an attempt to improve flow two-dimensionality in 
the test section. The fences, made of 20 cm x 70 cm float glass, were mounted 
approximately 3 cm from the sidewalls on standoffs which permitted adjustment 
during rig operation. The fences were found to be only moderately successful 
in reducing flow three-dimensionality. 

3 . 1 . 2 . 2 ^e^t_Se^c^i^n_Ija s^rjmien_ta_t ijon 

The test section instrumentation consisted of 140 wall static pressure 
taps and 75 0.63 cm dia. probing ports. Of these 120 of the wall static taps 
and 55 of the probing ports were located on the test surface. As shown in 
figure 3-5, the static taps were arranged in three rows of forty located on 
the test section centerline and 15.2 cm to either side of the centerline. In 
addition 35 of the probing ports were located on the centerline and 10 on 
either side of the centerline (figs. 3-6 and 3-7). The remaining 20 wall 
static taps and 20 probing ports were located along the centerline on the 
lower wall of the test section. The probing ports on the test surface were 
used for traversing total pressure probes and hot-film probes normal to the 
test surface to obtain flowfield data. Specific locations of the ports, which 
were machined into the test surface after the position of the separation 
bubble was established during shakedown testing, are described in detail 
below. Probing ports on the lower wall were used for traversing total 
pressure probes normal to the surface. These data permit calculation of the 
displacement thickness development and hence define the lower bounding stream- 
line for CFD codes. Additionally, one sidewall of the test section was fabri- 
cated from float glass to facilitate LV measurements. 
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3.2 Test Program Definition 


3.2.1 Selection of Measured Variables 


Directly measured and derived variables were selected based upon code 
assessment requirements. A viscous-invisc id calculat ional procedure for pre- 
dicting two-dimensional boundary layer separation and reattachment on a flat 
plate requires input data to permit calculation of both the potential flow- 
field within the test section and the viscous flow along the test plate. The 
freestream velocity at the test section inlet and the location of the bounding 
streamline along the lower wall opposite the test surface were required as 
boundary conditions for use in the potential flow solution. The integral 
boundary layer parameters (displacement thickness and momentum thickness) and 
the Reynolds number based on momentum thickness at the test section inlet were 
required to define the initial station for the viscous calculation. Depending 
on the complexity of the turbulence model used, the Reynolds stress distribu- 
tion at the test section inlet may also be needed to initialize the model. 

The desired output of the calculat ional procedure is a complete specification 
of the mean velocity field and static pressure distribution throughout the 
test section, and depending on the turbulence model used, some turbulence 
field properties. 


Based on these and experimental considerations, axial and vertical mean 
velocity components, total pressure, wall static pressures, and inlet total 
temperature were selected as the mean variables to be measured directly. 
Turbulence quantities measured directly were: the axial component of turbu- 


lence intensity •V?: the vertical component of turbulence intensity and 

the Reynolds stress, uv. Turbulence spectra were taken to document the 
unsteady nature of the separation bubble. In addition, the forward flow frac- 
tion, Yp^> the fraction of the time the flow moves downstream, and the trans- 


verse flow fraction, , 


the fraction of the time the flow moves away from 


the test surface, were selected for measurement to define the zones of separa- 
tion and reattachment and to delineate the regions where hot-film and total 
pressure measurements would be reliable. 


3.2.2 Tunnel Operating Conditions 


The test conditions for the experiment were selected based upon two cri- 
teria: the turbulent separation bubble should be large enough to provide high 

resolution of bubble details and the Reynolds number of the turbulent boundary 
layer on the test surface at the test section inlet should be in the range 
applicable to gas turbine airfoils. Establishment of a large separation 
bubble was achieved by providing a diverging-converging test section of 150 cm 
axial extent and a suitably scaled test section height. Reynolds number 
simulation was established by producing a turbulent boundary layer having a 
Reynolds number based on momentum thickness of 15,000 at the test section 
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inlet at the inlet velocity of 27 m/sec. This value is typical for aircraft 
turbine blades at the onset of cove separation. 

Separation bubble reproducibility ( s t at ionarity) throughout the duration 
of the test program was ensured by adjusting the operating conditions to 
achieve the same Reynolds number for each test run. This was achieved by 
adjusting the tunnel dynamic head, Q, to compensate for changes in barometric 
pressure and tunnel temperature. The dynamic head was determined from the 
difference between the total pressure measured with a kielhead probe located 
in the core flow of the boundary layer development section and a wall static 
pressure measured at the same axial location. Tunnel air temperature was 
measured with a Cr-Al sheathed thermocouple inserted through an insulated plug 
in the tunnel sidewall. The tunnel dynamic head was adjusted mainly for 
changes in barometric pressure and to a lesser extent for tunnel temperatures. 
The tunnel temperature varied from 22 to 27 deg C between winter and summer 
testing but on any given day the tunnel temperature could be maintained con- 
stant to within ± 0.5 deg C. 

3.2.3 Selection of Measurement Locations 


To meet the code assessment objectives, a "benchmark" study of two- 
dimensional turbulent boundary layer separation and reattachment on a flat 
plate must, as a minimum, provide the input data required to perform the 
computation and sufficient flowfield data at relevant locations to determine 
prediction accuracy. Supplementary data should also be obtained to assess 
experimental accuracy and to document flow two-dimensionality. 

To perform the calculation, data must be obtained at an initial measure- 
ment plane to start the computation and at a sufficient number of stations to 
define a bounding streamline. A plane normal to the upper wall, located in 
the boundary layer development section 24 cm upstream of the test section 
inlet, was selected as the initial measurement plane. The wall static pres- 
sure on the test surface and a 70-point total pressure traverse normal to the 
test surface were used to define the one-dimensional mean velocity profile at 
this location. A 34-point hot-film traverse was also used to determine the 
axial turbulence intensity profile. To determine the coordinates of a bound- 
ing streamline, wall static pressures and total pressure profiles normal to 
the wall were measured at 20 stations along the lower wall on the tunnel cen- 
terline (fig. 3-6). The measurement increment normal to the wall was chosen 
to yield a minimum of 20 points in each total pressure traverse. 

To assess the accuracy of the code predictions and to acquire a detailed 
understanding of the flowfield upstream, within, and downstream of the separa- 
tion bubble, a dense array of static pressure taps and traverse locations was 
provided on the test plate as shown in figure 3-7. Wall static pressure taps 
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extended the length of the test plate from the upstream zero-pressure gradient 
region to the corresponding region downstream of separation. Thirty-six 
traverse stations were located on the tunnel centerline (Table 1). Twelve 
stations were located upstream of separation, five within the intermittent 
separation zone, eight within the reversed flow region, five within the 
intermittent reattachment zone, and six downstream of separation. Ten 
additional traverse stations were located on a line 15.2 cm to either side of 
the tunnel centerline to document flowfield two-dimensionality. Two of the 
off-centerline traverse stations were located in each of the five distinct 
flow regions on the test plate. 

On the tunnel centerline, total pressure traverses were made at each 
of the 36 measurement stations, single-element hot-film probe traverses were 
made at 30 stations, and LV traverses to measure axial and vertical velocity 
components were made at 16 stations, as shown in figure 3-7. Along the line 
located 15.2 cm to either side of the tunnel centerline, total pressure 
traverses were made at each of the 10 measurement stations, single-element 
hot-film probe traverses were made at 8 stations, and LV traverses were made 
at four stations. Except for off-centerline LV data taken at the x = 26.7 cm 
(x = 10.5 in) station, each LV traverse station coincided with a hot-film 
traverse station, and each hot-film traverse station coincided with a total 
pressure traverse station. An additional 20 LV traverses were made at 14 
centerline stations and six off-centerline stations, as shown in figure 3-7, 
to measure Reynolds stress. 

The measurement increment normal to the test surface was varied based on 
resolution requirements. Fine resolution was necessary for defining the 
velocity field in the near-wall region since this is critical for code assess- 
ment. Additional detail was considered desirable to demarcate the edge of the 
reversed flow region. Since more than 300 traverses were planned in this 
experiment, regular, fine, and extra fine increment schedules were defined to 
limit the number of data points required to produce the desired resolution. 

The regular increment schedule was sufficiently dense to provide a minimum of 
10 LV and 20 hot-film or pitot measurements per traverse. It was used for all 
hot-film measurements and LV Reynolds stress measurements and for pitot and LV 
velocity measurements at off-centerline locations to document flow two- 
dimensionality. The fine increment schedule, which provided sufficient detail 
in the near-wall region to produce a 30- to 40-point traverse, was used for 
all LDV velocity measurements on the tunnel centerline and most off-centerline 
pitot measurements. The finest increment schedule, which produced a 50- to 
80-point traverse, was used for pitot traverses to define the boundary layer 
throughout the test section along the test plate centerline and at off- 
centerline locations near the test section inlet. 
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In selecting the measurement matrix described above, it was recognized 
that pitot and hot-film measurements could be inaccurate in regions of high 
turbulence intensity and would not be meaningful in regions having intermit- 
tent flow reversal more than a few percent of the time. The LV measurement 
locations were selected to define such regions and therefore permit the pitot 
and hot-film measurements to be interpreted correctly. In regions where the 
pitot measurements were valid, finer resolution of the near-wall velocity 
field was obtained. In conjunction with the LV mean velocity measurements, 
pitot data provided a means to estimate the static pressure field within the 
flow. In regions where the hot-film measurements were valid, turbulence 
spectra were acquired to contribute to an understanding of the unsteady nature 
of the turbulent separation bubble. In addition, the measured components of 
mean velocity and rms turbulence intensity provided a redundant set of data to 
check the LV measurements. 

3.2.4 Flow Two-Dimensionality and Steadiness Criteria 

Objectives were identified at the outset of the experiment regarding 
spanwise uniformity of the flow at the initial measurement plane and at three 
measurement planes in the test section: at the test section inlet, just 

upstream of separation, and within the reversed flow region. The spanwise 
uniformity criterion at the initial measurement plane was defined as deviation 
of less than one percent in the measured pitot pressure at tunnel midheight 
during a traverse across the tunnel (excluding a 7.6 cm region at each side- 
wall). Within the test section the criterion for spanwise uniformity of the 
flow at each measurement plane was agreement to within five percent of the 
boundary layer integral properties measured on the test plate centerline and 
at stations 15.2 cm to either side of the centerline. 

Criteria were also defined relative to the two-dimensionality of the 
separation line and for steadiness of the flow in the test section. A satis- 
factory two-dimensional separation was defined as a separation line deviating 
from a straight line by less than 5 percent of the test section span excluding 
a 7.6 cm region at each sidewall. The criterion for flow steadiness in the 
test section was defined as the absence of intermittent flow reversal on the 
contoured wall opposite the test plate and on the test plate outside the 
separated flow region as indicated by surface mounted tufts. 


3.3 Instrumentation 


3.3.1 Laser Velocimeter System 

The LV system (fig. 3-8) is described in detail in Appendix A. The sys- 
tem consisted of a 2W argon-ion laser, backscatter optical system, counter- 
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type signal processor, computer for on-line data reduction, and a hard disk 
storage device for subsequent off-line data reduction. All system components 
are commercially available. 

The LDV was operated in a dual beam or "fringe" mode in which light from 
the intersection of two incident beams having a wavelength of 0.5145 ym is 
heterodyned to detect the Doppler shift from an injected seed particle, at the 
local, instantaneous fluid velocity. In this mode, the LV measures the velo- 
city component in the plane of the incident beams that is perpendicular to the 
bisector of the beams. The effective shape of the resultant measurement 
volume is an ellipsoid with major axis~in the direction of the bisector of the 
beams. Sketch A below shows these features and the theoretical measurement 
volume dimensions assuming that the ellipsoidal surface is defined by the 
locus of points where Doppler signal amplitude is 1/e 2 of its centerline 
value. No direct measurements of the effective measurement volume dimensions 
were made; they can be affected by signal amplitude and the signal processor 
threshold level setting. To obtain maximum resolution, the major axis was 
aligned in the spanwise direction. Resolution of the LV system was high when 
compared to the relevant dimensions of the experiment. The minor axis of the 
measuring ellipsoid, d m , was 1/850 of the inlet boundary layer thickness and 
the major axis, l m , was 1/560 of the test section width. 

The optical system used a 3.7 5X beam expander, a relatively wide angle 
lens ( k = 4.9 deg), and 152 mm dia. optics to produce a sufficiently high sig- 
nal/noise ratio (SNR) to permit measurements to be made in the backscatter 
mode. A Bragg cell was used to eliminate directional ambiguity and to provide 
a 360 deg acceptance angle. 

The counter-type signal processor featured user selectable measurement 
mode, fringe count, validation accuracy, input gain, threshold level setting, 
and large amplitude signal rejection level. The processor was operated in the 
continuous (CONT) mode in which a measurement is made each time the minimum 
number of fringes is counted. The CONT mode was chosen to counter the effects 
of individual realization biasing (see discussion in Appendix C). Optimal 
histograms and fast data rates were achieved with the fringe count selector 
set at 32 and the accuracy limit for data validation set a 2%. The signal 
gain and large amplitude signal rejection level were adjusted to reject 10 to 
25% of the Doppler bursts to ensure that particle lag biasing (which could be 
caused by large seed particles in the tail of the seed particle size distribu- 
tion) was eliminated. Data rates of validated signals ranged from 20/ sec to 
more than 2000/sec depending on the region of the flow being measured. Based 
upon the uncertainty analysis in Appendix C and the large turbulence levels 
generated by the separation bubble, 2000 samples were acquired per data point 
to generate velocity histograms for mean velocity and turbulence intensity 
measurements. For Reynolds stress measurements the previously measured local 
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turbulence intensity was used to determine whether 2000, 4000, or 8000 samples 
per data point were required to obtain the desired accuracy. Additional 
details concerning LV data processing are given in Appendix A. 

Titanium dioxide powder having a nominal particle size of 1.0 jjm dia was 
dispersed into the tunnel airstream to provide highly reflective light scat- 
tering centers for LV measurements using the seeding system described in 
Appendix A. A mixture of titanium dioxide powder and Cab-0-Sil, a silica 
deaglomerating agent was fluidized in a low pressure seeder and injected into 
the tunnel just downstream of the aftermost turbulence suppression screen 
using a seeder probe designed to minimize the disturbance to the tunnel 
airstream. Particle lag calculations and LV velocity measurements across a 
normal shock using similar seeding material (see Appendix C) indicated the 1.0 
ym dia titanium seed particles would follow the mean flowfield variations with 
a lag of less than 0.3 percent throughout the test section. Less than 1 
percent error in turbulence intensity was anticipated for turbulence 
fluctuations below 2.9 KHz. 

3.3.2 Total Pressure Instrumentation 


The total pressure measurement system, shown in figure 3-9, consisted of 
total pressure probes described below, a probe traverse system, pressure 
transducers, signal conditioners and amplifiers, averaging voltmeters, a data 
logger, an on-line computer, and a micromanometer. The probe traverse system, 
described in reference 28, consisted of an ball/screw traverse drive suspended 
on a linear ball bearing track above the test section. Three 0-0.5 psid pres- 
sure transducers were used to measure total pressure, test section wall static 
pressure, and the tunnel dynamic pressure. Transducers were calibrated each 
day with the micromanometer. The measured pressures were averaged for 100 sec 
using averaging voltmeters and were recorded together with the tunnel tempera- 
ture on a data logger. 

Two sets of specially designed pitot probes, shown in figure 3-10, were 
used for total pressure traverses. The two sets, which differed only in the 
offset distance (1.27 cm or 1.91 cm) from the sensing face to the probe shaft 
centerline, each consisted of a boundary layer probe for near wall measurement 
and a kielhead probe for midstream measurements. The boundary layer probe was 
constructed from .081 cm dia stainless steel tubing flattened to an oval 
cross-section at the probe face having major and minor axes of 0.122 cm and 
0.043 cm respectively. A 0.64 cm dia. brass sleeve plugged the probing port 
during near wall traversing and was pushed out by and locked into the probe 
shaft for the remainder of the traverse. The kielhead probe was fabricated 
from a 0.159 cm dia. aspirated kielhead tube soldered into a 0.064 cm dia. 
probe shaft. The aspirated kielhead tube contained four 0.038 cm dia. holes 
equally spaced around its periphery at a distance of 0.191 cm downstream from 
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the probe face and a concentric inner sensing tube (0.038 cm ID). The four 
aspiration holes functioned in the same manner as the open back end of a 
conventional kielhead probe by permitting a through flow of air past the 
sensing tube. The kielhead tube was bent into alignment with the midrange of 
the mean flow directions anticipated during a traverse to ensure that the flow 
direction relative to the probe face was within the ± 30 deg acceptance angle 
determined during probe calibration. 

3.3.3 Hot-Film Instrumentation 


The hot-film measurement system consisted of a single element hot-film 
probe, a commercially available constant temperature anemometer system 
(including a fourth-order polynomial linearizer and a signal conditioner con- 
taining high- and low-pass filters), two rms voltmeters to time average the 
mean and fluctuating output signals, and a spectrum analyzer. The single- 
element hot-film probes (fig. 3-11) each used an 0.051 cm dia. cylindrical 
sensor to measure the streamwise mean velocity and turbulence components. 
Frequency response of the anemometer system was measured to be 65 KHz in a 
27.5 m/sec airstream at the 1.5 overheat ratio used in the test program. The 
unique design of the probe was necessitated by the desire to insert the probe 
through the same wall on which boundary layer velocity profiles were to be 
measured. A pin through the probe stem set the minimum stand-off distance at 
0.127 cm In the unbraced configuration supplied by the vendor severe 
‘’ringing 11 caused by vibration of support needles was noted during probe 
calibration. Nonmetallic braces were used to stiffen the support to eliminate 
ringing. 

3.3.4 Flow Visualization 


Flow visualization was accomplished during shakedown testing with sur- 
face-mounted tufts, ”tuft trees”, surface flow visualization, and smoke. 

Tufts were mounted on the test surface in a dense 10 x 14 array to monitor the 
size and two-dimensionality of the separation bubble, on the diverging wall 
opposite the test plate to check that the flow remained attached, and on the 
tunnel sidewalls. A 3 x 3 array of "tuft trees” was mounted vertically 
between the lower tunnel wall and static pressure taps in the test plate. The 
”tuft trees” were constructed from stiff 0.089 cm OD stainless steel tubing 
over which swivel-mounted tufts were slipped and secured in place by short 
lengths of shrink-fit tubing. Approximately 10 tufts were positioned on each 
’’tree” to visualize the flow in the vicinity of the separation bubble and to 
estimate the bubble thickness. Surface flow visualization was achieved by 
injecting red low-viscosity fluid through static pressure taps in the test 
plate. The smoke flow visualization system, described in reference 29, used a 
traversible probe upstream of the test section inlet to produce ammonium 
sulfite smoke by impinging gaseous jets of anhydrous ammonia and sulfur 
dioxide within the tunnel airstream. 
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3.4 Measurement Approach 


3.4.1 Total Pressure and Hot-Film 


Total pressure and hot-film traverses were aligned normal to the test 
plate. The total pressure traverses were conducted using boundary layer 
probes over the range 0.025 cm to 5.08 cm from the test plate and kielhead 
probes were used from 2.54 cm from the plate to the vicinity of the opposite 
wall. During the total pressure traversing phase, which was conducted before 
the flowfield had been defined by LV measurements, the flow direction near the 
reversed flow region was not known a priori at each point in space. There- 
fore, two sets of traverses were conducted to obtain a complete set of data at 
each measurement plane. One set of traverses was made with probes facing 
upstream and the second set made from the adjacent upstream traversing port 
with the probes facing downstream. Using probes with either 1.27 cm or 1.91 
cm offsets to match the 2.54 cm or 3.81 cm probe port separation distance 
ensured probe face coincidence in the measurement plane. Hot-film traversing 
was conducted in a similar manner in the vicinity of the separation bubble. 

After the hot-film anemometers were checked at the start of each day an 
additional linearizer adjustment was made to offset the effects of hot-film 
aging. The hot-film sensor was positioned at the x = 7.6 cm, y = 10.2 cm 
location in the tunnel where the streamwise velocity was known to equal 0.9925 
of the freestream velocity calculated from rig operating conditions. After 
the linearizer output gain was adjusted to yield the desired velocity, the 
probe was moved to the next scheduled traverse location. 

3.4.2 Laser Velocimeter 


The four LV beam orientations, designated 0, 90, and ± 45 deg, used to 

measure U, V, and uv are shown in figure 3-12. The direction of 

motion of the Bragg-shifted fringes is also indicated in the figure. In the 0 
deg orientation the vertical fringes measure only the horizontal, U, velocity 
component. In the 90 deg orientation the horizontal fringes measure only the 
vertical, V, component. In the ± 45 deg orientations, / 2 (U ± V), is measured 
from which uv can be obtained using the equations derived in Appendix B. 


Because of the 9.8 deg included angle between the laser beams, the 
optical axis was inclined 5 deg relative to the horizontal test plate surface 

for measurements of V, and uv near the test surface (y 2.54 cm). The 

inclination of the fringes relative to horizontal introduced a negligibly 
small uncertainty in the acquired velocity measurement due to the component of 
the transverse velocity. The 90 deg orientation measured V cos 5° + W sin 5°. 
The ± 45 deg orientations measured / 2 (U ± V cos 5° ± W sin 5°). The data 

reduction equations required to extract V, and uv from these measured 

quantities are derived in Appendix B. 
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Figure 3-1 United Technologies Research Center Boundary Layer Wind Tunnel 


CONTRACTION SECTION 



Figure 3-2 Modified Boundary Layer Wind Tunnel 
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Figure 3-3 Turbulent Boundary Layer Separation Test Section 
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Figure 3-4 Turbulent Separation Bubble Test Section Installed in UTRC Boundary 
Layer Wind Tunnel 
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Figure 3-5 Top of Test Section 


82 - 10 - 80-3 


29 






STATIC PRESSURE (120 TAPS) 



I 

I 


I 

I 


£ 

o 

CD 

LO 



• O < x I 




84 - 8 - 1-2 



i — to 


31 


Figure 3-7 Measurement Locations on Test Plate 
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Figure 3-8 Laser Velocimeter Data Acquisition System 
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CHAPTER 4 


RESULTS 


4,1 Data Presentation 

Velocity and pressure data obtained in this study are presented in both 
tabular and graphical format. Velocity data have been normalized by the free- 
stream velocity, U re f > at the tunnel reference location (x = -156 cm) in the 
constant velocity approach section upstream of the test section. u re f equaled 
27.05 m/sec (88.75 ft/sec) at nominal tunnel operating conditions (T re f = 
22.2°C, 72°F and P re f = 780 mm Hg, 15.06 psia where P re f was the sum of the 

tunnel static pressure relative to ambient (18 mm Hg, 10 in 1^0) and the 

barometric pressure (762 mm Hg )} . U re f was adjusted to maintain constant unit 
Reynolds number at the reference location when tunnel conditions deviated from 
nominal. Static pressures are referenced to the tunnel static pressure at the 
reference location. Total pressures were normalized by the freestream dynamic 
pressure, Q re f } at the reference location. The coordinates used for presenta- 
tion of results are x, distance downstream from the start of the test section 

(see figure 3-6), y, transverse distance from the test plate, and z, spanwise 
position relative to the wind tunnel centerline. 

Presentation of the results in the following sections starts with a 
description of the flow conditions at the initial measurement plane. Then, 
the results of smoke, tuft, and surface flow visualization are presented to 
provide an overall understanding of the test section flowfield. Detailed 
presentations are then made of the test surface static pressure distribution, 
the mean velocity field, turbulence velocity field, Reynolds stresses, and 
total pressures throughout the test section. 


4.2 Initial Plane Boundary Conditions 

To provide data required to start computational codes, an initial 
measurement plane was chosen well upstream of the region of influence of the 
adverse pressure gradient and the resulting separated boundary layer. 

Although such codes should be capable of proceeding downstream from an 
arbitrary initial profile, the approach geometry for the present study 
produced a boundary layer having characteristics generally similar to an 
equilibrium turbulent boundary layer, as discussed below. The initial 
measurement plane was located in the boundary layer development section 24. 1 
cm (9.5 in) upstream of the test section inlet. This location, which was 2.5 
duct heights upstream of the start of the diverging section, was chosen 
because test surface static pressure measurements, described subsequently, 
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showed that the adverse pressure gradient extended less than 1.5 duct heights 
upstream of the divergence. Pitot probe traverses, both transverse and span- 
wise, were taken to define the flow uniformity in the initial measurement 
plane. Hot-film anemometry was used to define profiles of the mean velocity 
and the axial turbulence component in the test surface boundary layer. 

The measured total pressure distribution is shown in figure 4-1. The 
spanwise pitot traverse taken at the tunnel mid-heighth showed the total 
pressure variation to be less than ± 0.5 percent of the dynamic head over the 
central 90 percent of the span which was outside the sidewall boundary layers. 
A transverse (vertical) pitot pressure traverse at mid-span was used to define 
the boundary layers on the test (upper) surface and the lower surface. The 
boundary layer edge was defined as the location at which the measured total 
pressure equaled 99 percent of the freestream total pressure (equivalent to 
6gg 5 determined from velocity measurements). The boundary layer on the test 
surface, which was tripped at the start of the boundary layer development 
section (fig. 3-2), measured 7.67 cm (3.02 in) whereas the untripped lower 
wall boundary layer had a thickness of 3.56 cm (1.4 in), leaving a 9.9 cm 
(3.9 in) thick inviscid core. Integral properties calculated for the boundary 
layer on the test surface at the initial measurement plane are tabulated 
below: 

Displacement thickness, 6 = .799 cm 

Momentum thickness, 0 = .631 cm 

Shape factor, H = 1.27 

Momentum thickness Reynolds no. = 11,100 

The measured shape factor of 1.27 was somewhat less than the value of 
1.32 quoted by Coles (ref. 30) for a constant pressure turbulent boundary 
layer at the measured momentum thickness Reynolds number of 11,100. The 
velocity profile of the test surface boundary layer was found to be in good 
agreement with the " law-o f-t he-wall." In figure 4-2a the profile is plotted 
in n law-of-the-wal 1 M (U + vs. y + ) coordinates. To obtain figure 4-2a the fric- 
tion velocity, U , was chosen to minimize the least squares fit to the 
equation 


U U * CL41 ln y + 5,0 


(4-1) 


+ yu T 

where y = — - — as recommended by Coles (ref. 30) over the log-linear region of 
the profile (50 < y+ < 500). Then the skin friction coefficient could be 


2 2*2 

calculated from the equation Cf = 2U^/U^. e £ since U T = x, 


Data in the outer 


region of the boundary layer deviated from the logarithmic portion as 
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expected. The maximum deviation, AU + , referred to as the "strength of the 
wake component" was approximately 1.4. This was less than the value of 2.7 
quoted by Coles for an equilibrium turbulent boundary layer having the same 
momentum thickness Reynolds number. 

The measured streamwise turbulence intensity (normalized by boundary 
layer edge velocity and plotted versus y/6g g ) is shown in figure 4-2b to be in 
good agreement with the turbulence profile presented by Klebanoff (ref. 31) 
for a flat plate equilibrium turbulent boundary layer. In addition to overall 
rms measurements, spectra (0-20 KHz) were taken routinely at each hot-film 
measurement point to ensure that the hot-film signal was free of discrete 
frequency spikes associated with sensing element vibration. 


4.3 Flow Visualization 

To optimize the two-dimensionality of the separation bubble, during 
initial flow visualization testing, surface-mounted tufts and tuft trees 
(installed in the test section and described previously in Section 3.2) were 
monitored while adjustments were made to the divergence angle of the lower 
wall, the boundary layer bleed scoops, and the axial distribution of suction 
in the test surface sidewall bleed slots. It became apparent that two- 
dimensionality could only be achieved using a large divergence angle and 
suction in both the lower wall bleed scoop and the corner bleed slots. This 
configuration resulted in a long, thick separation bubble having minimal 
three-dimensionality. 

Since the technique of achieving the bubble configurations employed in 
this experiment may be of interest, some initial exploratory studies will be 
described. By varying lower wall divergence angle and suction arrangements 
for the two lower wall slots, a range of bubble configurations could be 
produced. While the overall bubble length was effectively set by the axial 
separation of the lower wall divergence and convergence points, bubble thick- 
ness could be reduced by decreasing the divergence angle. This tended, 
however, to produce a more unsteady flow and greater spanwise variation in the 
separation line. This does not imply that thinner bubbles are inherently more 
unsteady, rather that a smaller adverse gradient will increase the axial 
length over which the random process of intermittent separation can occur. 

The configuration producing the highest degree of bubble two-dimensionality in 
the present experiment was found to be one in which the divergence angle was 
increased to the maximum value (a = 28 deg) for which slot suction could hold 
the boundary layer attached to the lower surface. Although the resulting 
thick bubble caused significant streamline curvature in the x-y plane with 
resultant transverse static pressure gradients, the strong-adverse-pressure- 
gradient, strong-suction configuration was considered superior for providing 
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benchmark quality data for code verification and development compared to a 
configuration producing a shorter, thinner bubble having unknown three- 
dimensional contamination. 

Tuft flow visualization of the test surface is shown in figure 4-3. The 
five zones of the separation and reattachment process (attached flow, inter- 
mittent separation, reversed flow, the reattachment zone, and reattached flow) 
are clearly evident to the eye but are obscured somewhat in the photography. 
The inset, taken from a movie made during flow visualization testing, shows a 
closeup view of the tufts within the reversed flow region to illustrate the 
strong reversed flow across the entire tunnel span within the separation 
bubble. The maximum thickness of the reversed flow region within the separa- 
tion bubble as indicated by the direction of the tufts on the "tuft trees" 
which extended from the upper to the lower surface was estimated to be 8 cm 
(3.1 in). Also indicated by the "tuft trees" is the curvature of the free- 
stream flow convecting over the separation bubble. Surface mounted tufts on 
the lower wall confirmed that the boundary layer along the lower wall remained 
attached throughout the test section. 

In figure 4-4 the five zones of separation are depicted more clearly. 

In the attached flow region all the tufts were steady and pointed downstream. 
Within the intermittent separation zone, tufts fluctuated in direction between 
upstream and downstream. In the reversed flow region, tufts in the central 60 
percent of the span were steady and pointed upstream. In the reattachment 
zone, the tufts fluctuated in direction between upstream and downstream as the 
reattachment point shifted in a random, intermittent manner. In the 
reattached flow the tufts again pointed downstream but jittered due to high 
turbulence levels in the wake of the separation bubble. 

Photographs of the smoke flow visualization obtained while traversing the 
single-headed smoke probe vertically at inlet station 0 (x = -24.1 cm, -9.5 
in) are shown in figure 4-5. For the first photograph (fig. 4-5a) the smoke 
probe was positioned 10 cm below the test surface (y = 10 cm) just outside of 
the 7.67 cm thick boundary layer at the initial measurement plane = The smoke 
filament flowed smoothly through the test section but was displaced downward 
by the separation bubble on the test surface. Note that the axial extent of 
the separation bubble as indicated by tufts mounted on the test surface is 
marked on the photographs. At y Q = 7.6 cm (a location corresponding to the 
edge of the boundary layer at the initial measurement plane where the smoke 
was injected) the filament continued to flow smoothly around the bubble (fig. 
4-5b) , indicating that it was in the low turbulence, inviscid freestream. 
However, at y Q = 6.4 cm (fig. 4-5c) the upper edge of the smoke filament 
became entrained in the turbulent shear layer surrounding the separation 
bubble. Near reattachment the smoke filament appears diffuse due to unsteady 
flapping of the filament lines in that region. When the smoke filament 
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impinged upon the test surface it was observed that the filament would 
alternately be ingested into the recirculation region at the aft end of the 
bubble or be convected downstream. With the smoke probe positioned near the 
upper wall (y Q = 5.1 cm, fig. 4-5d) the smoke filament became more entrained 
in the separation bubble. Much of the entrainment occurred in the free shear 
layer surrounding the bubble as shown in the photograph where a large eddy is 
being entrained into the bubble in addition to the flow being ingested into 
the bubble near reattachment. Once within the bubble, the smoke dispersed 
rapidly causing the smoke density within the bubble to be too low to be seen 
in the photograph although the smoke was visible as a diffuse haze to the eye. 

Flow directions determined from tuft trees, surface tufts, and smoke 
flow visualization were used to map the streamline pattern throughout the 
test section (fig. 4-6). The reversed flow region was quite large, being 
approximately 50 cm long with a maximum thickness of 8 cm. Streamlines 
upstream of the separation zone were steady and parallel. Near reattachment 
the streamline pattern for the mean flowfield could only be approximated 
because of unsteady flapping of the flow-field. 

To determine the degree to which a two-dimensional separation had been 
achieved, tufts were removed from the test section and a low-viscosity red dye 
was injected through static pressure taps on the test surface (fig. 3-7 shows 
tap locations). The measurements confirmed the straightness of the separation 
line indicated by the surface mounted tufts earlier. However, the determina- 
tion of a quantitative measure of the straightness of the separation line was 
difficult because the low wall shear in the vicinity of separation barely 
moved the injected dye. After testing, a sheet of exposed blueprint paper was 
rolled over the test surface to lift off the dye streaks for permanent storage 
and subsequent analysis. The dye patterns were interpreted later with the aid 
of LV measurements of the forward flow fraction (see Section 4.6.1 "Forward 
Flow Fraction Measurements") at centerline and off-centerline locations near 
the test surface. Good correlation between a forward flow fraction of 0.5 and 
lack of motion of the injected dye was obtained. The separation line (i.e. 
location of transitory detachment) was determined from the dye patterns and 
interpolation of forward flow fraction data to be straight. The angle of 
inclination of the straight separation line was 7 degrees with respect to the 
perpendicular to the tunnel centerline with separation occurring further 
upstream on the west side of the span (z > 0) than the east side (z < 0). 

This 7 degree inclination angle corresponds to a deviation from the line 
perpendicular to centerline of ±6 percent of the span which is reasonably 
close to the ±5 percent objective. 
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4.4 Test Surface Static Pressure Distribution 


The static pressure distribution on the test surface measured along 
the tunnel centerline is shown in figure 4-7. The results are presented in 
the form of a non-dimensional pressure coefficient defined as 


c p = 


Pw (x) - P Wl 


^ref 


(4-2) 


where p w (x) is the wall static pressure at any x location, p w is the wall 
static pressure measured by the first static pressure tap in the test section 
(x = 6.35 cm) and Q re f is the dynamic pressure. Also shown in the figure are 
streamlines determined from flow visualization testing and curves estimating 
the Cp distribution at separation and reattachment based on an analysis by 
Perry and Fairlie (ref. 5). 

The static pressure on the test surface began to rise downstream of the 
test section inlet (x > 10 cm) due to the upstream influence of the 28 deg 
lower wall divergence which began at x = 30 cm and created a strong adverse 
pressure gradient. At x = 40 cm C p reached a value of 0.3 6 and the boundary 
layer separated. The wall static pressure remained constant under the separa- 
tion bubble except near reattachment where the impingement and partial stagna- 
tion of the reattaching flow caused C p to increase to 0.50 on the test 
surface. Downstream of reattachment C p decreased rapidly as the flow 
accelerated in the converging test section leading to the outlet duct. Near 
the test section exit C p approached a constant value of 0.22 due to the net 
diffusion which occurred in the flow between the inlet and outlet of the test 
section. As shown in the streamline patterns drawn in figure 4-7, when the 
flow which was removed from the test section via the lower wall bleed scoop is 
excluded from consideration, the effective inlet height, h£ n , of the test 
section was less than the outlet height, h out , resulting in a net flow 
diffusion within the test section. 

Perry and Fairlie noted that the C p distribution near separation and 
reattachment could be approximated by parabolae which were defined by the 
separation (or reattachment) angle, <f> , and the nondimens ional vorticity, u> 
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1/2 


t an <f> = ± 2 



(4-3) 


The parabolae shown in figure 4-7 were estimated using measured values of the 
Cp distribution and separation and reattachment angles estimated from the 
streamline patterns calculated below from LV data (fig. 5-3). The vorticity 
was assumed to be constant along the mean dividing streamline and was 
estimated to be the value which resulted in the best parabolic fit to the 
measured Cp distribution at reattachment. 

Test surface static pressure distributions were also measured at off- 
centerline locations as one of the diagnostics for establishing the degree of 
flow two-dimensionality. Results of the off-centerline are tabulated in Table 
IV beginning on page 192. 


4.5 Mean Velocity Field 

The velocity field within the test section was measured with a laser 
velocimeter (LV) and a single-component hot-film probe. The LV measured axial 
and transverse velocity components separately whereas the hot-film probe 
measured the magnitude of the velocity vector in the x-y plane. Hot-film 
measurements were made to provide a degree of redundancy in the measurement of 
the velocity field and to provide information about the unsteady nature of the 
flowfield. Additional hot-f ilm measurements were made at stations upstream 
and downstream of the separation bubble where optical inaccessibility pre- 
vented LV measurements. All hot-film data, including data obtained in regions 
of intermittent flow reversal where the measurements are known to be inaccur- 
ate, are presented in Table V. The hot-film velocity data is presented 
nondimens ionally as u e ff/ u re f an< * u eff^ U ref w ^ ere U eff anc * u eff are t * ie mean 
and rms values of the effective cooling velocity measured by the hot-film 
probe which was aligned such that the sensor axis was parallel to the z-axis. 

Comparison of these data to the more accurate LV data could prove useful in 
evaluating separation bubble data obtained by other investigators using the 
hot-wire/f ilm technique. A comparison between hot-film and LV data obtained 
at selected stations during this test program is presented in Appendix D. 

Error bounds for the LV measurements determined using the error analysis 
developed in Appendix C are presented in Table C-5. The complete set of LV 
data including mean velocity components, turbulence components, forward and 
transverse flow fractions, and Reynolds stress measurements are presented in 
Tables II and III. Error bounds for pitot and hot-film measurements are given 
in Appendix D. The velocity data presented below will consist almost entirely 
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of LV data except for hot-film data obtained at stations where LV measurements 
could not be made due to optical inaccessibility. 

Velocity vectors determined from LV measurements taken on the wind tunnel 
centerline have been plotted in figure 4-8 to show the mean flowfield in the 
test section. Several features of the flowfield are notable. First, the size 
and location of the separation bubble are in agreement with the flow visuali- 
zation results. Separation (i.e., U < 0 near the wall) occurred at approxi- 
mately x “ 45 cm with reattachment (i.e. U > 0 at y > 0) at x = 100 cm 
resulting in a 55 cm long separation bubble. The strong recirculation within 
the separation bubble is readily apparent. The maximum thickness of the 
reversed flow region is almost 9 cm (i.e., slightly larger than the 8 cm 
thickness indicated from flow visualization) causing a significant 
displacement and curvature of the freestream flow. Second, near the leading 
edge of the bubble at 45 cm < x < 60 cm, a strong free shear layer is evident 
between the reversed flow near the surface and the displaced boundary layer. 
Finally, at 90 cm < x < 100 cm, the flow vectors impinge nearly vertically 
onto the test surface. This impingement caused the peak in the surface static 
pressure at reattachment which has been discussed in the previous section. 

4.5.1 Axial Velocity Measurements 

Axial velocity profiles nondimens ionalized by ^ re f are presented in 
figure 4-9. All data were obtained from LV measurements except the data at 
Station 1 (x = 2.54 cm) which are streamwise velocity components measured with 
a hot-film probe. 

The general mean flow characteristics of the separation-reattachment 
process are evident in this figure. Retardation of flow near the wall is 
evident downstream of Station 6 with detachment occurring at Station 13. 

Within the reversed flow region, shown bounded by a dashed line, peak 
velocities are observed near the wall although retardation of the reversed 
flow at the surface is evident in the forward half of the region. Outside of 
this reversed flow region significant (apparent) profile adjustments are 
observed over the axial extent of the reversed flow zone. Downstream of 
Station 26 reattachment occurs with the development of a new boundary layer 
downstream suggested by the near-wall data points. 

More detailed axial velocity profiles are presented in figures 4-10 ard . 
In these figures the axial velocity profiles are plotted in the non- 
dimensional form U/U e vs. y/6 where U e is the local boundary layer edge 
velocity at y - 6 determined from laser velocimeter measurements and 6 is the 
local boundary layer thickness determined from total pressure traverses as 
described below in Section 4.9. The data are plotted using semilogar ithmi c 
coordinates near the wall (y/6 < 0.1) and linear coordinates outside the near- 
wall region to permit the closely spaced data points in the near-wall region 
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to be presented clearly. The profiles are also displaced for clarity. 

Symbols adjacent to left and right borders show individual profile ordinates 
of zero and 1.0, respectively. By comparing the profiles in part (a) to the 
zero values the retardation of the near wall layer between Stations 0 and 11 
is evident. Part (b) of the figure shows that at Station 13 a maximum 
reversed flow velocity of 0.10 U e (0.075 U r(jf ) has been measured. The maximum 
reversed flow velocity at each station increases monotonically with distance 
downstream until Station 22 (located approximately 2/3 of the bubble length 
aft of the detachment location) where a value of 0.397 U e (0.202 U ref ) is 
reached. As shown in part (d) of the figure, reattachment and the resumption 
of downstream flow occurs between Stations 26 and 28. At Stations 15 through 
24 values of U/U e greater than unity were measured within the boundary layer 
(y/6 < 1). These results, caused by normal static pressure gradients produced 
by flow curvature in the test section, are analyzed in Sec. 5.3.2. 

4.5.2 Transverse Velocity Measurements 

No nd imensional transverse velocity profiles, y/6 vs V/U ref , obtained 
from LV measurements are presented in figures 4-11 a-c. Positive transverse 
velocities correspond to flow away from the test surface. Figures 4-1 la and 
4-1 lb show that at Stations 6 through 16 in the upstream half of the test 
section, the transverse velocity, V, increases monotonically with y/6 as the 
flow moves away from the test surface. Positive transverse velocities within 
the reversed flow region are consistent with the expected counter clockwise 
bubble recirculation pattern. Downstream of Station 16 the transverse 
velocity component decreases to a value near zero at Station 20. Downstream 
of Station 20 negative transverse velocity components are measured indicating 
flow toward the test surface. As shown in figure 4-llb and 4-llc the compo- 
nents exceed 40 percent of U re £ in the reattachment zone located in the vicin- 
ity of Stations 28 and 30. At Station 34 in the exit duct the transverse 
component magnitude has decreased to the 0-5 percent range. 

In summary, the measurements show that significant transverse velocity 
component magnitudes occur within the boundary layer regions defined by total 
pressure measurements. As in the case of the axial velocity measurements, 
profile alterations with increasing axial distance are observed to be signifi- 
cant. 


4.6 Intermittency Measurements 

The directionally sensitive laser velocimeter was used to measure the 
fraction of the time the flow was moving in the downstream direction (U > 0) 
and the fraction of the time the flow was moving away from the test surface 
(V > 0). Simpson et al (ref. 1) have referred to the former measurement as 
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the intermittency and/or the forward flow fraction. In this report the latter 
terminology will be adopted and "transverse flow fraction" will be used to 
denote the fraction of the time the flow moved away from the test surface. In 
practice these quantities were determined from the ratio of the number of 
valid samples having positive velocities to the total number of valid samples 

rnmnri cino oarVi T \J mpacnrPTTiPnf 
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4.6.1 Forward Flow Fraction Measurements 


Figure 4-12 shows the distribution of the forward flow fraction through- 
out the test section. Corresponding tabular information is provided in Table 
2. Upsteam of separation (Stations 5 and 8) the flow fraction was unity 
across the entire profile indicating that the flow was moving in the down- 
stream direction 100 percent of the time. The first indication of intermit- 
tent reversed flow occurred near the test surface at Station 11 where y p 

F U 


first deviated from 1.0. Further downstream the backflow region intensified 
and spread outwards from the test surface. At Stations 18 through 22, a 
strong steady reversed flow (y p = 0) extended several centimeters from the 


test surface. The region of the flow in which reversed flow existed for some 

fraction of time (y p < 1) extended more than 17 cm from the test surface. 

F U 


Notice also that the locus of points at which y p 


U 


= 0.50 corresponded closely 


to the boundary of the reversed flow region. The correspondence would be 
exact if the LV probability distributions where symmetric with respect to 
U = 0. However, near the edge of the reversed flow region, the velocity 
histograms were skewed because of the nature of the flow. Since the maximum 
backflow velocity was much less than the freestream velocity, the probability 
of a turbulent fluctuation having a large positive velocity was greater than 
the probability of it having a large negative velocity causing a skewed 
h istogran. 


More details of the development of the separated flow can be seen in 
figure 4-13 in which the forward flow fraction is plotted versus y/6 in log- 
linear coordinates for stations within and adjacent to the separation bubble. 
As noted above, occurrences of reversed flow appeared initially at Station 11 
near the test surface. Downstream of Station 11 these occurrences became more 
frequent (y p decreased) and penetrated further into the viscous layer as 


defined by the local boundary layer thickness, 6. At Stations 15 and 16 the 
forward flow fraction distribution formed a trough near the wall with the more 
intense regions of reversed flow located slightly away from the test surface 
(0.01 < y/6 < 0.05). 


Downstream of Station 16 the trough-shape disappeared as strong, steady 
reversed flow was measured adjacent to the test surface at Stations 18 through 
22 over almost 10 percent of the boundary layer thickness. These data are 
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consistent with the velocity profile data shown in figure 4-10 which showed 
the most intense reversed flow which exceeded 30 percent of the local boundary 
layer edge velocity occurred at Stations 18, 20 and 22. 

Near reattachment, trough-shaped patterns reappeared as the reverse flow 
weakened. The first reoccurrence of forward flow near the test surface at the 
reference near-wall measurement location (i.e. y - 0.10 cm) appeared at 
Station 24 where y p ^ = 0.012 at y/6 = 0.0034. The near-wall forward flow 

fraction increased to 0.144 at y/6 = 0.0042 and 0.714 at y/6 = 0.0050 at 
Stations 26 and 28, respectively. At Station 30 the reattachment process was 
completed and a steady flow profile was reestablished. The axial extent of 
the reattachment zone from the occurrence of forward flow (vt> > 0) to Yt> = 

y . u 

0.50 near the wall (i.e. = 0) has been determined by data interpolation to 

be approximately 12 cm. The complete reattachment zone extended from Station 
24 to Station 30, a distance of 22.8 cm. For comparison, the boundary layer 
thickness at Station 8 just upstream of separation was 11 cm. A final point 
which should be noted in figure 4-13 is the absence of occurrences of reversed 
flow in the outer 40 percent of the boundary layer at any measurement 
station. 

The trough-shaped pattern noted near separation was similar to the flow 
fraction distribution measured by Simpson et al (ref. 2) at the downstream 
stations in their turbulent separation experiment. They did not achieve steady 

reversed flow (y D = 0) anywhere in their test section. The minimum forward 
u 

flow fraction which they measured, y D = 0.05, occurred at y/6 = 0.023 at their 

F U 

final measurement station. Thus, the open separation bubble data of Simpson 
et al . would not be expected to correlate well with closed separation bubble 
data obtained downstream of the leading edge of the bubble. 

4.6.2 Transverse Flow Fraction Measurements 

The distribution of the transverse flow fraction throughout the test 
section is shown in figure 4-14. The transverse flow fraction data are 
consistent with the mean LV flow directions presented above (fig. 4-8). Near 
the test surface y p - 0.50 except at Stations 28 through 31 where the mean 

flow impinged strongly onto the test surface. The transverse flow fraction 
value of 0.5 occurs when the mean flow vector is approximately parallel to the 
test surface since transverse turbulence fluctuations cause the flow to have 
an equal probability of moving toward or away from the test surface. 

The direction of the freestream flow as indicated by the transverse 
flow fraction is also consistent with the LV velocity data. In the upstream 
half of the test section (Stations 6-18) the freestream flow moved away from 
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the test surface 100 percent of the time (v^ = 1), 

F V 


In the downstream 

the test section (Stations 22-31) the freestream flow moved toward the 
surface 100 percent of the time (y p ^ = 0). 


half of 
test 


4.7 Unsteady Flowfield 

Flow visualization had shown the flowfield in and around the 
separation bubble to be unsteady. Much of the unsteadiness, such as the 
streamline flapping at reattachment and the oscillation of the separation 
line, seemed to be of large scale and occurred at low frequency. It was not 
clear from flow visualization whether such large scale fluctuations were 
random (i.e. low frequency turbulence) or ordered (i.e. periodic unsteadi- 
ness). In addition, to random versus periodic considerations, the question of 
whether fluctuations should be considered "turbulence* 1 or "flow unsteadiness** 
is important to developers of computational procedures. This arises because 
measured unsteadiness levels can be interpreted in terms of turbulence level 
(and hence effect turbulence modeling) or interpreted as an unsteady flow with 
lower levels of imbedded turbulence. 

To obtain some information regarding the nature of the unsteadiness in 
the flowfield, the frequency content of the flow was measured with a hot-film 
probe during traverses at three critical regions of the flowfield. Frequency 
data were obtained at Station 13 in the intermittent separation zone, at 
Station 18 within the reversed flow region, and at Station 28 in the reattach- 
ment zone. Typical spectra of the hot-film signal taken over the 0-500 Hz 
range are shown in figure 4-15. Except for a single discrete frequency spike 
in the spectrum obtained in the freestream (y = 27 cm) at Station 18, all 
spectra were broadband indicating random unsteadiness. The single low ampli- 
tude spike at 460 Hz was caused by blower noise generated at the blade passing 
frequency. Additional 0-20 KHz spectra taken during each hot-film measurement 
to check for probe vibrations revealed no flow periodicity at high frequen- 
c ies . 

4.7.1 Streamwise Turbulence Measurements 


The broadband spectra showed considerable strength in the low frequency 
range (f < 50 Hz). To quantify this low frequency contribution to the turbu- 
lence signal, two rms measurements of the linearized hot-film signal were 
made. For the first rms average, which was performed for each hot-film 
measurement made in the test program, the prefilters on the rms voltmeter were 
set to average the signal over the wideband range of 5 Hz-50 KHz. In 
addition, the voltages obtained at Stations 13, 18, and 28 were averaged on a 
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second voltmeter to permit determination of the rms contribution in the low 
frequency 5-50 Hz range. 


Turbulence intensity profiles (y vs. V? ff/U re f) calculated from the rms 
voltage averaged over the 5 Hz-50 KHz bandwidth are shown in figure 4-16 as 
solid lines. The mean square turbulence level u^~ff can be considered the 
turbulence level in the x-y plane associated with the local mean flow direc- 
tion. Note that the turbulence intensities have been normalized bv U c 

. J ref 

rather than U e to provide absolute rather than relative turbulence intensi- 
ties. Upstream of separation the flow had a standard Klebanof f-like turbu- 
lence profile as shown previously. High turbulence levels exceeding 20 
percent were measured in the free shear layer surrounding the separation 
bubble. These peak turbulence levels decayed somewhat downstream of reattach- 
ment but remained above 10 percent at the test section outlet. 

At Stations 13, 18, and 28 the portion of the overall streamwise 
turbulence level attributed to the low frequency (5 Hz-50 Hz) band has been 
shaded. This low frequency contribution to the overall turbulence level was 
substantial across the entire profile. In the free shear layer the low 
frequency (5 Hz-50 Hz) unsteadiness accounts for approximately 40-45 percent 
of the overall turbulence level (5 Hz-50 KHz). Uncertainty in these measure- 
ments exists within regions of flow reversal. 


Turbulence data at Stations 13, 18, and 28 are listed in Table V. 
Measurement points where flow reversal occurs some fraction of the time (0 
< < 1) are indicated in the table. At those points the rms average of the 


hot-film voltage does not yield an accurate measurement of the turbulence 
level (see Appendix D). Thus, care should be taken when using the data in 
Table V. 


With regard to the issue raised above regarding interpretation of 
measured unsteadiness as turbulence or a basically unsteady flow with imbedded 
turbulence, the above measurements were not definitive. Some form of condi- 
tional sampling may be required to address this issue more directly. That 
large scale turbulence with low frequency content existed within and in the 
vicinity of the recirculation zone is clear based on the smoke flow visualiza- 
tion movie taken as part of the study. Overall flowfield unsteadiness was 
also apparent. This was particularly noticeable in the aft end of the bubble 
where a previously described non-periodic flapping of the reattaching flow 
occurred. This flapping phenomenon has been reported previously by Eaton and 
Johnston (ref. 32), for a reattaching backward-f ac ing step flow, Cherry et al. 
(ref. 33), for a separating and reattaching flow at the leading edge of a two- 
dimensional rectangular blunt plate, and Mabey (ref. 34) among others. Mabey 
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showed that reduced shedding frequencies from various experiments correlated 
well when normalized by the reattachment length. He suggested that a feedback 
mechanism, whereby the shedding of vorticity from reattachment interacted with 
the formation of disturbances near separation, controlled the separation 
process. In the present experiment a characteristic disturbance length can be 
estimated from the flapping frequency, f < 50 Hz, and the average freestream 
velocity in the vicinity of the bubble, U e = 15 m/sec, to be greater than 
0.3 m (i.e. half the length of the separation bubble). 

An additional degree of unsteadiness would be expected in the current 
experiment because of the intermittent nature of turbulent boundary layer 
separation on a flat plate. The random upstream and downstream wandering of 
the separation line due to the lack of a salient edge to fix the location of 
separation presumably added to the downstream unsteadiness. Thus, there is 
sufficient evidence to conclude that the unsteadiness observed was not an 
artifact of the wind tunnel test configuration employed here. 

4.7.2 LV Turbulence Measurements 


Accurate measurements of the turbulence levels throughout the test 
section, including regions of flow reversal, were made with the laser veloci- 
raeter. An overview of the axial and transverse turbulence intensity profiles 
determined from LV data is shown in figure 4-17. Qualitatively, the axial 
turbulence intensity profiles are similar to the streamwise turbulence inten- 
sity profiles from hot-film measurements presented in figure 4-16. The 
transverse turbulence intensity profiles have the same general characteristics 
as the axial turbulence profiles although lower in level. Peak transverse 
turbulence intensities were approximately 12 percent while the corresponding 
peak axial turbulence intensities exceeded 20 percent. 

4 . 7 . 2 . 1 A.x_i-al_ C_omp_on^en_t ^ujrbjal^n^e__Mje a_su_r emejit_s 

More detailed axial turbulence intensity profiles are shown in figures 
4-18 a-d where the data are plotted in nondimens ional coordinates, y/6 vs. 

Vu^*/U,-~.f: . At Station 0, located 24 cm upstream of the test section, the 

L ^ i. 

boundary layer had a K1 ebanof f-1 ike turbulence distribution with maximum 
turbulence intensity of approximately 8 percent as noted previously. As the 
flow moved downstream and approached separation, turbulence levels near the 
test surface increased rapidly. At Station 11, the start of intermittent 
separation, turbulence reached a level of 17 percent at y/6 - 0.10. Just 
downstream of the mean separation location (Station 13) peaks in the turbu- 
lence intensity profiles grew higher and wider, and moved away from the test 
surface due to turbulence being generated in the free shear layer and 
diffusing throughout the boundary layer. The peak turbulence intensity of 21 
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percent was measured at Station 15 at y/6 = 0.34, Downstream of Station 15 
the peaks in the profiles decayed, continued to widen, and shifted outwards to 
larger values of y/5. Over the upstream half of the separation bubble the 
peak turbulence intensity occurred approximately 0.1 6 farther from the test 
surface than the mean dividing streamline (see Table VI, Streamline No. 1 for 
coordinates of the mean dividing streamline determined from LV measurements). 
Over the downstream half of the separation bubble the peak turbulence inten- 
sity shifted further outside of the mean dividing streamline. 

At all stations the axial turbulence intensity in the freestream 
approached 2 percent or less except at Station 34 at the test section outlet 
where upper and lower wall boundary layers began to merge causing the free- 
stream turbulence level to rise to 4 percent. The freestream turbulence level 
of 2 percent is a quantitative measure of the overall unsteadiness of the 
freestream flow in the test section. 

4 . 7 . 2 . 2 ^r^n^v^r^e_C^m^oiie-nt__Tiirb^ul_eiice_ Mea_su r emen ts 

Detailed nond imensional transverse turbulence profiles are shown in 
figure 4-19 a-d. At initial measurement Station 6 the turbulence distribution 
for this component again exhibited a standard K1 ebanof f-1 ike profile with a 
maximum intensity of approximately 4 percent. Approaching separation, trans- 
verse turbulence levels increased near the test surface. At the leading edge 
of the separation bubble the transverse turbulence level peaked at 8.2 percent 
at y/6 = .17. Just downstream of separation at Stations 14 through 16 the 
peak transverse turbulence level remained constant at approximately 8 percent 
but widened and diffused outward from the test surface. 

The behavior of the transverse turbulence intensity at the next several 
measurement stations provides insight into the nature of the transverse 
velocity fluctuations. At Station 18 the transverse turbulence intensity 
peaked at 12.6 percent. The peak declined to 9.6 percent at Station 20 before 
increasing to 11.8, 12.6, and 12.7 percent at Stations 22, 24, and 26, respec- 
tively. The mean transverse velocity profiles plotted in figure 4-llb 
indicated positive transverse velocities at Station 18, negligible transverse 
velocities at Station 20, and negative transverse velocities at Stations 22, 
24, and 26. Thus it seems possible that a major contribution to transverse 
turbulence intensities originate as streamwise velocity fluctuations when the 
mean streamline had a nonaxial orientation. 

The high transverse turbulence levels at Stations 26, 28 and 30 near 
reattachment also resulted from unsteady streamline flapping. Downstream of 
reattachment transverse turbulence intensity decreased due to turbulence decay 
and a more axial orientation of the streamwise flow vector. Transverse 
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turbulence levels in the freestream approached 2 percent or less at all 
stations, 

4.7.3 Turbulence Kinetic Energy 

The turbulence kinetic energy per unit mass, k, is defined as 

k = (u 2 + v 2 + w 2 )/2 (4-4) 


Although the spanwise turbulence component, w, was not measured in this 
experimental program, Simpson (ref. 4) has shown that in separated flows 
v — w . Thus, k can be approximated as 

k = (u 2 + 2 v 2 )/2 (4-5) 


Turbulence kinetic energy profiles calculated using eq. (4-5) are presented in 
figure 4-20. The profile shapes are somewhat similar to the axial turbulence 
profiles discussed above. The profiles had sharp and growing peaks near 
separation (Stations 11-13), remained at nearly constant peak amplitude but 
broadened near the maximum bubble thickness (Station 18-22), and decayed 
downstream of reattachment. The profile peaks occurred just outside the 
location of the mean dividing streamline over the upstream half of the bubble, 
whereas downstream of the maximum bubble thickness, they moved outward from 
the mean dividing streamline into the bubble wake. 


4.8 Reynolds Stress Measurements 

Profiles of the Reynolds stress, uv, were measured with the LV system 
at 20 stations. Fourteen of the stations were located on the wind tunnel 
centerline and six were located at ± 15.2 cm on either side of the centerline 
at three axial locations (upstream of the separation bubble, within the 
bubble, and downstream of reattachment). Reynolds stress data from thirteen 
centerline locations is listed in Table III. Profiles obtained at twelve 
stations along the wind tunnel centerline are plotted in figure 4-21. Five 
features of the Reynolds stress profiles should be noted. First, at Station 6 
the Reynolds stress has a K1 ebanof f-1 ike profile (fig. 18.5 in ref. 35). 

Second, in the freestream flow around the bubble the Reynolds stress is 
negligible as expected. Third, the Reynolds stress peaks near reattachment 
and remains high downstream of reattachment. Fourth, the Reynolds stress is 
small in the reversed flow region. And fifth, just upstream of separation and 
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in the shear layer surrounding the bubble no meaningful trend is apparent in 
the Reynolds stress data. 

The almost random nature of the Reynolds stress profiles in the separated 
flow region is not particularly elucidative. However, the Reynolds stress 
data in this coordinate system can be used to calculate shear stress profiles 
in a coordinate system aligned with the local mean flow direction. This 
streamwise shear stress is more meaningful for developing flow structure 
models, since the assumptions used in the boundary layer equations which are 
used to calculate separating flows are more accurate if applied in streamwise 
orthogonal coordinates (ref. 6). It will be shown below in Section 5.7, 
"Streamwise Turbulent Shear Stress Measurements" that the large values of 
shear stress expected in the shear layer are masked by the effects of flow 
curvature and they do not become fully apparent until the streamwise shear 
stress is calculated. 


4.9 Total Pressure Measurements 

Total pressure traverses were made normal to the test surface at 55 
stations. Thirty-five stations were located on the wind tunnel centerline and 
twenty traverses were made at off-centerline locations as shown in figure 3-7. 
The total pressure measurements obtained at each centerline traverse location 
and at fourteen off-centerline locations are listed in Table IV. 

The total pressure profiles obtained on the wind tunnel centerline are 
shown in figure 4-22. For clarity, only 18 of the 34 profiles have been 
plotted. The measured total pressure has been normalized by Q ref , the free- 
stream dynamic pressure at the inlet station (x = -9.5 in, -24.1 cm). Several 
features of the total pressure profiles are apparent. 

The inviscid core, defined as the region over which 0.99 < P/Q re f < 1.00, 
which is approximately 10 cm thick upstream of the location of intermittent 
separation (x < 40 cm), remained intact as it was convected around the bubble. 
In the outlet channel (x > 120 cm) the thickness of the inviscid core 
diminished but did not disappear until the last measurement station (x = 148 
cm) where the upper and lower wall boundary layers merged. 

The deceleration of the test surface boundary layer in the region 
0 < x < 40 cm is also apparent. The shape of the boundary layer profile 
changed as the flow encountered the adverse pressure gradient until an 
inflection point was measured in the profile at the wall at x = 40.6 cm just 
upstream of separation. 
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Within the reversed flow region 3p/3y - 0 except very near the test 
surface where a reversed flow boundary layer developed under the separation 
bubble. The presence of the reversed flow boundary layer is not apparent in 
figure 4-22 but it can be discerned in the tabulated data in Table IV and it 
is discussed in detail below in Section 5.5.1, "Description of the Backflow 
Boundary Layer M . 

Downstream of reattachment a new boundary layer began to grow along the 
test surface beneath the remnants of the upstream boundary layer which had 
been convected around the separation bubble. Details of the boundary layer 
downstream of reattachment are discussed in Section 5.6, "Description of the 
Reattachment Process". 

4.9.1 Determination of Bounding Streamline 


Total pressure traverses were made normal to the lower wall of the test 
section at 20 locations on the wind tunnel centerline to determine the loca- 
tion of a free streamline required as a boundary condition for numerical 
computation codes. Using wall static pressures measured at each of the 20 
traverse locations and assuming that no normal static pressure gradients 
existed between the lower wall and the location of the bounding streamline, 
mass flow profiles were calculated from the measured total pressure profiles 
at each of the 20 axial stations. After accounting for the mass flow out of 
the lower wall boundary layer scoop, points of equal mass flux were connected 
to define the bounding streamline shown in figure 4-23. 
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Figure 4*1 Spanwise and Transverse Total Pressure Distributions in Inlet Plane 
(Station 0, x= -24.1 cm, -9.5 in.) 
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Figure 4-2 Boundary Layer Characteristics at initial Measurement Plane 
(Station 0, x= -24.1 cm, -9.5 in.) 
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Figure 4-4 Surface Flow Directions on Test Plate as Indicated by Tufts 
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Figure 4-5 Smoke Flow Visualization of Separation Bubble 
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Figure 4-5 Smoke Flow Visualization of Separation Bubble (Concluded) 
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Figure 4-6 Streamlines Determined from Smoke Flow Visualization and Tuft Trees 



AXIAL LOCATION, 



Figure 4-7 Test Surface Static Pressure Distribution on Tunnel Centerline 
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Figure 4-10a Non-Dimensional Axial Velocity Profiles. Note Log-Linear Abscissa and 
Displaced Ordinates 
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Figure 4-10b Non-Dimensional Axial Velocity Profiles. Note Log-Linear Abscissa and 
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Figure 4-1 Oc Non-Dimensional Axial Velocity Profiles. Note Log-Linear Abscissa and 
Displaced Ordinates (Continued) 
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Figure 4-10d Non-Dimensional Axial Velocity Profiles. Note Log-Linear Abscissa and 
Displaced Ordinates (Concluded) 
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Figure 4-1 1b LV Measurements of Transverse Velocity Profiles (Continued) 
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Figure 4-1 3a Forward Flow Fraction Profiles • Note Log-Linear Abscissa and 
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Figure 4*15 Spectra of Streamwise Turbulence Component Measured with Hot*Film 
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Figure 4-21 Nondimensional Reynolds Stress Profiles Determined with Laser Velocimeter (z-0) 
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CHAPTER 5 


ANALYSIS AND DISCUSSION 


Chapter 4 presented the results from direct measurements along with their 
discussion. In this chapter, additional quantities derived from the measure- 
ments will be analyzed and discussed. These quantities include flowfield 
streamline coordinates, static pressure, skin friction, and streamwise shear 
stresses. Boundary layer integral properties and law-of-the-wall plots are 
also presented. The measured and derived data are then analyzed and compared 
to the results and correlations developed in previous investigations. In 
doing this the flowfield has been partitioned into eight regions or zones 
because the flat plate separated and reattached boundary layer has not been 
studied previously in its entirety with LV. Data in each flow zone will be 
compared to available data obtained in other separation and reattachment 
studies . 


5.1 Flow Zones 

The flow zones are defined in figure 5-1. Zone I is the attached 
boundary layer upstream of separation. Boundary layer profiles in Zone I can 
be compared to studies concerning the development of turbulent boundary layers 
in adverse pressure gradients. The attached boundary layer flows against the 
adverse pressure gradient until the kinetic energy in the boundary layer at 
the wall is converted into potential energy at which point the boundary layer 
separates from the test surface. Since the boundary layer kinetic energy at 
any point in a turbulent boundary layer is proportional to the square of the 
sum of the steady velocity field and the unsteady turbulence velocity field, 
separation occurs over a zone rather than at a point as occurs in a laminar 
boundary layer. The intermittent separation zone. Zone II, has been studied 
extensively for open separation bubbles by Sandborn, Kline, Simpson, and many 
others. Hence, comparisons can also be made for this region of the flowfield. 
Zone III, the turbulent separation bubble, is defined as the region between 
the mean dividing streamline and the test surface. Strong backflows exceeding 
30 percent of the local freestream velocity occur within this region. The 
open separated flow data of Simpson et al. can be used as a basis for com- 
parison in the upstream half of the bubble and the downstream half of the 
bubble can be compared to separation data obtained downstream of a backward 
facing step. The strong recirculating flow within the separation bubble 
causes a backflow boundary layer, Zone IIIA, to originate near the 
reattachment point and to develop along the test surface. No basis for 
comparison exists in this zone because details of the backflow boundary layer 
in a closed separation bubble have not been previously reported. 
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The reattachment of the separated shear layer is an inherently unsteady 
process which occurs within Zone IV, the reattachment zone, on the test 
surface. In this zone the reattaching flow impinges steeply onto the test 
surface causing a static pressure peak which drives the strong backflow within 
the separation bubble. The flow is characterized by high turbulence levels 
and frequent flow reversals. Boundary layer integral properties, skin fric- 
tion, and forward flow fraction data in the reattachment zone can be compared 
to the data of Driver and Seegmiller (ref. 14), Westphal et al . (ref. 17), and 
others for reattaching flows behind a backward facing step. 

Downstream of reattachment the flow recovers the structure of a flat- 
plate turbulent boundary layer. This process occurs in two stages. In Zone 
V, a new boundary layer begins to grow along the test surface and quickly 
develops a log-linear velocity profile near the wall. The relaxation of the 
outer free shear layer in Zone VI occurs much more slowly, however. Through a 
complicated process of turbulence energy exchange between Zones V and VI the 
turbulence structure downstream of reattachment gradually obtains the outer 
wake-like structure of a turbulent boundary layer. Data in Zones V and VI 
have been correlated with the recovering boundary layer profiles of Chandrsuda 
and Bradshaw (ref. 18), 

When properly normalized the axial velocity profiles in Zone VII, the 
outer flow region, are found to be invariant throughout the test section. 
Results in this zone can be compared to separated flow studies of Perry and 
Fairlie (ref. 5) and Paterson and Weingold (ref. 36). 


5.2 Static Pressure Field 

The static pressure field within the test section was calculated at 36 
stations. Twenty-six of the stations were located on the test plate center- 
line and five pairs of stations were located at ± 15.2 cm off the centerline. 
The measured quantities used to determine the static pressure, p, were the 
total pressure, P, and the mean axial and transverse velocities measured with 
the LV system 



where Q re f(P) is the dynamic pressure at the reference location during the 
measurement of P and U re £(U) and U re f(V) are the freestream velocities at the 
reference location during the measurement of U and V, respectively. As noted 
in Section 4.1, U re f was adjusted during testing to maintain constant unit 
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Reynolds number at the reference location. The results are presented in terms 
of the difference between the local static pressure and the static pressure at 
the reference location in the boundary layer development section, p , 

• ^ r e f 

normalized by the reference dynamic pressure, For these calculations 

the density, p, which varied less than 0.3 percent throughout the test 

section, was assumed to be constant. Calculated values of the static pressure 

are listed in Table IV. Static pressure profiles at 15 stations on the test 

plate centerline are plotted in figure 5-2. 

The transverse static pressure gradient, 3p/8y, was essentially zero at 
the inlet and outlet stations. It was negative between Stations 3 and 15 due 
to streamline curvature. Within the separation bubble (Zone III) the static 
pressure was nearly constant except at the wall near reattachment (as shown 
previously in fig. 4-7) where the static pressure increased due to flow 
impingement onto the test surface. Outside the separation bubble between 
Stations 18 and 24, 3p/8y > 0 due to the concave curvature of the streamlines 
in the freestream flow associated with flow around the bubble and the converg- 
ing lower surface. Between Stations 30 and 33 the convex curvature of the 
lower surface at Station 31 caused 3p/3y < 0 as the flow accelerated into the 


The normal static pressure gradients within the test section were signi- 
ficant. At some stations the static presure difference between the upper test 
surface and the lower wall equaled 0,5 Q re f* The static pressure gradients 
were caused by the streamline curvature in the inviscid flow due to the 
combined effects of the tunnel geometry and the aerodynamic blockage induced 
by the large separation bubble. The effect of the transverse static pressure 
gradients on the calculation of the boundary layer integral properties is 
assessed below in Section 5.3.3. 

Spanwise variations in the static pressure were minimal at most stations. 
The maximum variation in static pressure between the centerline and the off- 
centerline traverse stations, which occurred near the boundary layer scoop and 
at Station 34 downstream of reattachment , was less than .01 Q re f/ cm - 


5.3 Mean Flowfield Characteristics 

5,3.1 Mean Flowfield Streamlines 
— — T 

Axial velocity profiles measured with the laser velocimeter were inte- 
grated to determine the value of the streamfunction ^ £ at any point y£ 
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The locus of points of equal values of determined by interpolation, were 

used to define the mean flovfield streamline pattern shown in figure 5-3. 
Coordinates for the mean dividing streamline and eight additional streamlines 
located 1.27 cm apart at initial LV measurement Station 6 are listed in Table 
VI. The bounding streamline (dashed line) determined from lower wall total 
pressure traverses (see Section 4.9.1) has also been plotted in figure 5-3. 
Coordinates of the bounding streamline are given in figure 4-23 for use in 
computing the flowfield. The displacement of the mean flowfield by the 
separation bubble is quite evident. The maximum thickness of the separation 
bubble, h, measured from the surface to the ip = 0 mean dividing streamline was 
17 cm (6.7 in.) as noted in Table VI. Separation bubble length was 64 cm 
(25 in.). Significant streamline curvature is evident in the vicinity of the 
separation bubble. 

At most points within the test section the mean flowfield streamlines are 
tangent to the measured mean velocity vectors as expected in a steady flow. 
Nontangency between the streamlines and the mean flow vectors near the leading 
edge of the separation bubble is unexplained but may arise from individual 
realization bias errors (ref. 37) in the transverse velocity component 
measurements in regions of the flowfield having spatially non-uniform seed 
particle distributions. Such non-uniformities are noticeable near the separa- 
tion bubble where the seed density in the shear layer around the bubble was 
higher than the density within the bubble. Thus, positive, downward velocity 
fluctuations would provide relatively fewer valid counts than negative, upward 
velocity fluctuations because of the lower seed density near the test surface. 
This individual realization bias would tend to reduce the magnitude of the 
measured transverse velocity component. 

5.3.2 Inviscid Outer Flow 


Analysis of mean velocity data acquired in the present experiments has 
shown a feature of the velocity profile development which appears to have 
important implications for the numerical modeling of the flow. When static 
pressure variations associated with streamline curvature are taken into 
account, congruence of velocity profiles outside of the recirculating region 
is found over the axial extent of the separation bubble region. This feature 
will be discussed in this section. 


Axial velocity profiles measured on the test section with the laser 
velociraeter have been shown in figure 4-9. When the profiles from the various 
axial measurement stations are overlayed, station- to-stat ion similarity is not 


evident. Because of the large displacement of the flow streamlines caused by 
the separation bubble, the axial velocity profiles have been distorted by 
streamline curvature in addition to viscous effects. 

It is possible, however, to separate viscous from streamline curvature 
effects by normalizing the measured axial velocity component, U, by the 
inviscid axial velocity component, , at each point on the profile. The 
inviscid axial velocity component was calculated from the freestream total 
pressure, P re f, the measured transverse velocity component, V(y) , and the 
local static pressure, p(y), discussed in the previous section. 

u i<y) = ^ ^ P ref _ P ( y ) " \ P (5-3 ) 

This procedure for normalizing the axial velocity profiles by the local 
inviscid velocity has been used previously by others to calculate boundary 
layer integral properties in flows over curved surfaces and in shock boundary 
layer induced separations where normal pressure gradients are also significant 
(see Section 5.3.3. 1 "Definition of Integral Properties"). 

The normalized axial velocity profiles are shown in figure 5-4. These 
profiles have been overlayed in figure 5-5 using y-6 as the ordinate where y 
is the vertical distance from the flat plate test surface and 6 is the the 
local boundary layer thickness. For clarity, profiles at stations 18, 20, 22 
and 26 have been omitted from figure 5-5 although their profiles lie within 
the envelope formed by profiles at stations 16, 24, and 28. Using the 
boundary layer at station 6 (x - 26.7 cm, 10.5 in.) as reference, it can be 
seen that as the flow moves downstream the inner part of the initial boundary 
layer profile is eroded by the viscous processes associated with separation 
and reattachment. A similar observation has recently been made in a widely 
different separated flow experiment. As reported by Paterson and Weingold 
(ref. 36), flow downstream of a simulated compressor airfoil blunt trailing 
edge displays velocity profile congruence over the axial extent of the trail- 
ing edge separation bubble. In that case profile similarity was directly 
evident in the measured data since streamline curvature and associated static 
pressure gradients were significantly less than that encountered in the 
present experiment. 

These results suggest that the outer flow can be considered basically an 
inviscid, rotational flow which convects around the viscously dominated 
recirculation zone. This is considered to be potentially important in the 
modeling of separated flow by viscous-inv isc id interaction procedures. 
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5.3.3 Boundary Layer Integral Properties 


5. 3. 3.1 Definitions of Integra^ ^roper_t_ie_ s 


In classical boundary layer theory in which the static pressure is 
assumed to be constant across the boundary layer, the displacement thickness, 
5*, and the momentum thickness, 0, are defined with the boundary layer edge 
velocity as the reference velocity 
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As long as the boundary layer thickness to the streamline radius of curvature 
is less than about 1 percent, the static pressure difference across the 
boundary layer will be less than 2 percent of the dynamic pressure and these 
definitions are adequate (ref, 38). However, in the present experiment, as 
noted above, streamline curvature within the test section was significant. 
Static pressure changes across the boundary layer exceeding 40 percent of the 
dynamic head occurred at several measurement stations in the vicinity of the 
separation bubble. 

Kiock (ref. 38) has considered the interpretation of boundary layer data 
when sizeable transverse pressure gradients exist across the boundary layer. 

He reviewed five definitions for both the displacement and momentum thickness 
and compared their calculated values in example flows having prescribed normal 
static pressure gradients. In the present experiment, boundary layer integral 
properties have been calculated from data using three of the five definitions 
reviewed by Kiock. 

The displacement thickness is defined as the thickness of the hypotheti- 
cal layer in which the fictitious inviscid flow is equal to the reduction of 
the real mass flow due to the boundary layer. Since the mass flux between the 
displacement surface and the edge of the boundary layer in the fictitious 
inviscid flow must equal the total mass flow in the real boundary layer, the 
displacement thickness in the current experiment can be defined by the 
following equation 
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(5-6) 
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Udy 


where is the local inviscid axial velocity defined above in Section 
5.3.2. 

Kooi (ref. 39) has shown that equation (5-6) can be solved to obtain the 
true definition of 6* in the presence of normal static pressure gradients as 
the sum of the displacement thickness evaluated from classical boundary layer 
theory, 6^ , and a correction term to account for the effects of normal static 
pressure gradients 
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static pressure and the inlet total pressure, equals U e when 3 p/9y = 0. 

Another definition of the displacement thickness, referred to by Kiock 
(ref. 38) as the "pressure-based" displacement thickness, 



(5-8) 


has been used by So and Mellor (ref. 40), Mayle, et al . (ref. 41), and others 
to account for normal static pressure gradients due to streamline curvature. 

Kooi has shown that the true definition of the momentum thickness can be 
expressed as the sum of the momentum thickness evaluated from classical 
boundary layer theory, 0 and a correction term to account for the effects 
of the normal static pressure gradients. 


6 6 



The "pressure-based" momentum thickness for incompressible flow has been 
defined (ref. 38) as 
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(5-10) 
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dy 


The liminHarv 1 avpr shanp factor. H. defined as the ratio between the disnlace- 

j j r - * * i 

merit thickness and the momentum thickness can be calculated for each of the 
three definitions of the integral properties defined above. 

Kiock (ref. 38) used examples to show that when U-: increases monotoni- 

"k "k ic 

cally across the boundary layer ^ 2 ^ ^ 3 ^ 0 p^and H 2 < Hj. 

When decreases monotonically across the boundary layer 6 2 < '^3 < '^l>^2 < ^ 

63 < 9 p and H 2 > H 3 > . The equations of classical theory hold when is 

constant across the boundary layer and 6 - S 2 "(S 3 , 0 ^ = 0 2 " ® 3 anc * ^ 1 = ^2 35 


5. 3. 3, 2 £ a l.c_ola_t_ed__V£lu.es_ ^f_In t egral^ ^ropj^rj^ie_s 

Each of the three values of the displacement thickness, momentum thick- 
ness, and the shape factor defined above were calculated at each LV traverse 
location within the test section. Plotted values of the displacement and 
momentum thicknesses on the wind tunnel centerline are shown in figure 5-6a. 
The displacement thickness increased rapidly downstream of the start of the 
intermittent separation zone (x - 40.6 cm, 16 in) until it peaked at x = 74.9 
cm (29.5 in.), the location of the maximum thickness of the separation bubble 
(see fig. 5-3). The shape of the displacement thickness curve is almost 
symmetric about the center of the bubble. The momentum thickness increased 
much more slowly as the boundary layer approached separation. Within the 
separation bubble the momentum thickness remained relatively constant before 
decreasing downstream of reattachment. 

Values of the shape factor are plotted in figure 5-6b. The true value of 

the shape factor, , increased from 1.38 at x = 26.7 cm (10.5 in.) to 2.09 at 

x = 40.6 cm (16.0 in.) near separation, then to a peak value of 16.1 at x = 
74.9 cm (29.5 in.) before declining to a value of 1.23 in. the accelerating 
flow in the outlet duct. Peak values of and also occurred at x = 74.9 
cm and equaled 11.8 and 7.5 respectively. 

The results follow the trends shown in Kiock 1 s example boundary layers. 
Upstream of separation (x < 48 cm, 19 in.) and downstream of reattachment 
(x > 113 cm, 44.5 in.), where 3 U^/ 3 y > 0 across the boundary layer, 2 < 6 3 < 
6 *, 0 2 < 0 3 < 0 L and H 2 > H 3 > Hj . Near the center of the separation bubble 
(48 cm < x < 90 cm), 3U^/9y < 0 across the boundary layer and 6 2 < 63 < <5^, 

0 2 < 0 3 < 9 1# and H 2 > H 3 > Hj. At x = 48 cm, 19 in., au^/ay “ 0 across the 

boundary layer and 6 ^ = 63 = 63 , 0 ^ = 03 = 03 , anc ^ ^1 = ^2 = ^ 3 * t ^ ie 
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remaining measurement stations the sign of SUf/Sy varied across the boundary 
layer. Although Kiock demonstrated the relationships between classical, true, 
and pressure-based integral parameters with simple mathematical models of the 
boundary layer, the agreement shown here demonstrates that these trends hold 
equally well in complex flowfields of interest. 

5.3.4 Skin Friction Coefficient on Test Surface 


Law-o f-the-wall velocity profiles were calculated from pitot data and LV 
velocity data in the attached flow upstream of separation and downstream of 
reattachment. The skin friction coefficient was determined from the friction 
velocity, U T 


C 


f 



(5-11) 


Within the separation bubble the skin friction coefficient was determined from 
backward facing pitot data. At all stations within the bubble Cf was 
evaluated from the normal velocity gradient at the wall estimated from pitot 
data. Where possible Cf was also calculated from friction velocities 
obtained from law-o f-the-wall plots of the backflow velocity profiles. 

Details of the calculations are given in Section 5.5.1. 


The skin friction coefficients in the test section are presented in terms 

of the upstream reference velocity (i.e., C^ = 2 x w /pU^ e £) anc * a ^- so terms of 

the local freestream velocity (i.e., = 2 t ^plJ^ ) . Values of the skin 

L 

friction coefficient based upon the upstream reference velocity, Cf, will be 
used to quantify wall shear stress variations in the test section. C^ , based 

on the local freestream velocity, will be used to compare to Swafford’s corre- 
lation of skin friction in attached and separated turbulent boundary layers 
(ref. 42). 


Values of Cf on the test plate centerline are plotted in figure 5-7a. 

At the reference station 24 cm upstream of the test section inlet Cf equaled 
0.00267, a typical value for a zero-pressure gradient turbulent boundary layer 
having Re 0 = 11,100. Cf decreased in the decelerating boundary layer 
upstream of separation until it reached a value of 0.0001 at x - 40.6 cm just 
upstream of detachment. The negative value of Cf calculated from the 
backflow velocity profiles within the separation bubble indicated that the 
wall shear stress within the separation bubble was approximately 1/10 as large 
as the wall shear in the attached flow wall upstream of separation. Down- 
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stream of reattachment Cf increased rapidly until a value of 0.00282 was 
reached at the last measurement station in the exit duct. 


Swafford (ref. 42) has proposed the following empirical correlation which 
uses measured values of shape factor, H, and momentum thickness Reynolds 
number, p.eg , to estimate in both attached and separated flows. 
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(5-12) 


The first term in Swafford’s correlation is the skin friction correlation 
given by White (ref. 43) for attached turbulent boundary layers. Swafford 
added the second term to correlate the negative skin friction coefficients in 
reversed flows. Experimentally determined values of Cf L on the test plate 
centerline have been plotted in figure 5-7b together with Swafford’s correla- 
tion. Experimentally determined values of 0 2 and ^ were used to evaluate 
eq. (5-12). The agreement between the correlation and the data is very good 
up to separation. The correlation correctly predicts the length of the 
separation bubble but significantly underpredicts the skin friction in the 
backflow boundary layer. The maximum negative value of skin friction coeffi- 
cient (i.e., Cf L = -0.00022) which can be obtained from eq, (5-12) limits 
its ability to correlate the wall shear in energetic backflow boundary layers. 
Downstream of reattachment the Swafford correlation underpredicts the rapid 
monotonic increase in skin friction by 20 to 40 percent. Swafford anticipated 
that his correlation would not be universally applicable to all flow situa- 
tions since it was developed largely from the open separation bubble data of 
Simpson et al. (ref. 1) and Alber et al. (ref. 44). The Swafford correlation 
has not been compared previously to separation bubble data through reattach- 
ment . 


5.3.5 Demonstration of Flow Two-Dimensionality 

The criterion to measure the spanwise uniformity of the flow at the 
initial measurement plane, test section inlet, just upstream of separation, 
and within the reversed flow region have been defined in Section 3.2.4 ’’Flow 
Two-Dimensionality and Steadiness Criteria”. The attainment of the spanwise 
uniformity criterion at the initial measurement plane, which was defined as 
deviation of less than one percent in the measured pitot pressure at tunnel 
midheight, has been described in Section 4.2 "Initial Plane Boundary Condi- 
tions". The attainment of a straight separation line which deviated from a 
line perpendicular to the tunnel centerline by ± 6 percent of the test section 
span excluding the sidewall regions has been described in Section 4.3 "Flow 
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Visualization". This ±6 percent deviation only slightly exceeded the ±5 
percent criterion. Satisfaction of the remaining criterion that the boundary 
layer momentum thickness at each spanwise measurement station (i.e., z = -15.2 
cm, 0, +15.2 cm) deviate from the spanwise average by less than 5 percent will 
be discussed in this section. 


First, however, the spanwise uniformity of the total pressure profiles 
within the test section should be noted. In figure 5-8 nondimensionalized 
total pressure data are plotted vs. transverse distance from the test plate at 
three spanwise locations for each of six axial stations: two stations 

upstream of separation, two stations within the separation bubble, two 
stations downstream of reattachment. The profiles at Stations 2, 11, 18, and 
24, encompassing the region from upstream of separation to the aft end of the 
bubble, show good spanwise uniformity with the maximum profile difference 
being on the order of .05 p /Q re f Downstream of reattachment at Station 31 
(fig. 5-8e) the total pressure profile on the centerline has a larger defect 
at 2.5 cm < y < 7.5 cm than the off-centerline profiles. This difference is 
almost eliminated at Station 34 (fig. 5-8f) in the exit duct. 
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ized by the boundary layer thickness. Minimal spanwise skewness is evident in 
the profiles at Stations 18 (fig. 5-9e) and 24 (fig. 5-9d) indicating that the 
skewness observed in figure 5-8c and 5-8d resulted from small spanwise varia- 
tions in the boundary layer thickness within the separation bubble. The vari- 
ation in the boundary layer thicknesses measured at the z = +15.2 cm and 
z = -15.2 cm locations at stations 18 and 24 were 11.6 percent and 7.0 percent 
of the boundary layer thicknesses on the centerline, respectively. 


The spanwise variations in the momentum thickness at the same six axial 
stations has been assessed. Values of 0 i, the momentum thickness derived 
from classical boundary layer theory, and 0 3 , the pressure-based momentum 
thickness (see Section 5. 3.3.1 definitions of Integral Properties") 
normalized by the spanwise average are plotted in figure 5-10. At stations 2 
and 11 upstream of separation and Stations 18 and 24 within the reversed flow 
region 0^ varied less than 5 percent across the span. Downstream of 
reattachment at stations 31 and 34 the momentum thickness peaked on the wind- 
tunnel centerline and varied from the spanwise average by 8 and 13 percent, 
respectively. The pressure-based momentum thickness, 0^ , showed similar 
trends to 0 1 except at station 18 where 0^ was 8 percent less than the 
spanwise average at z = 0. The slight three-dimensionality in the flowfield 
downstream of reattachraent , apparent from the spanwise variations in momentum 
thickness, is consistent with the inward skewness of the surface tufts shown 
above in figure 4-4. 
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5.4 Description of the Separation Process 
5.4.1 Nature of the Decelerating Boundary Layer 

The streamwise mean-velocity distribution in a turbulent boundary layer 
can be expressed as a linear combination of wall and wake functions as shown 
by Coles (ref. 45). 

u + = f(y + ) + — w (5-13) 

where k is a constant to be determined from experimental data and w(y/6) is a 
universal wake function. The near wall flow is dominated by viscosity whereas 
the wake flow represents the cumulative effect of the pressure gradient. When 
the inner layer (or wall) structure is in equilibrium with the outer layer (or 
wake) structure the turbulent boundary layer is said to be in equilibrium. 

Coles showed that the wake parameter H could be expressed as 


n 



(5-14) 


He further noted that when H is constant the turbulent boundary layer is in 
equilibrium and there is a certain balance between the constraints imposed by 
inertia (i.e. U e ) and viscosity (i.e. U T ). A balance also exists between 
the large-scale and small-scale mixing processes (ref. 45). This concept of 
equilibria is useful in zero, favorable, and slightly adverse pressure- 
gradient, turbulent boundary layer flows whose mean-velocity profile and 
turbulence structure are dominated by the effects of the wall (ref. 6). 

Another equilibrium parameter, which characterized the flow in the outer 
layer, was derived by Clauser (ref. 46). 


T w dx 


(5-15) 


Clauser showed that a turbulent boundary layer with variable pressure 
gradient, but g constant, is in equilibrium. 
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White (ref. 43) studied 13 near equilibrium flows presented at the 1968 
Stanford Conference (ref. 47) and proposed the following empirical relation- 
ship between II and 8 for equilibrium flows 


n - 0.8 (6 + o . 5)°- 73 


(5-16) 


To determine the nature of the boundary layer upstream of separation, 
velocity profiles were plotted in law-o f-the-wall coordinates (fig. 5-11). 

The velocity profiles at Stations 0-5, were determined from pitot traverses 
with the assumption that transverse static pressure gradients were negligible. 
At Station 6 (x - 26.25 cm) profiles were determined from both pitot and LV 
data. Only LV data is presented at Station 8 (x * 31.75 cm) because of 
significant transverse static pressure gradients at that station. 

Values of skin friction, Cf , and the wake parameter, H , were calculated 
at each station from the profiles. The skin friction decreased monotonicall y 
from a value of .00267 at Station 0 (x = -24.1 cm) to .00077 at Station 8. 

Thp wake parameter* IT* equaled 0.33 at Station Q where the boundary layer was 
nearly in equilibrium (see Section 4.2). (Coles suggested a value of 0.55 for 
a flat-plate zero-pressure-gradient boundary layer.) The wake parameter 
remained relatively constant through Station 3 (x = 13.34 cm) before 
increasing rapidly to a value of 3.4 at Station 8 (x = 31.75 cm) as the flow 
encountered an increasingly adverse pressure gradient as shown previously in 
figure 4-7. 

The Clauser equilibrium parameter, 8, was also evaluated at the same 
measurement stations after eq. (5-15) was transformed into the following 
form. 


6* dC 

e = — — (5-17) 

c f dx 


The calculated values of 8 and II are plotted in figure 5-12 together with the 
empirical curve, eq. (5-16), defining the relationship between IT and 8 for 
equilibrium flows. As shown in the figure, although the boundary layer was 
nearly in equilibrium at Stations 0, 1 and 2, neither H nor 8 remained 
constant between Sections 3 and 8 as the boundary layer lost its equilibrium 
structure when it encountered the strong adverse pressure gradient. 

Some insight into the reasons for this loss of equilibrium can be 
obtained by differentiating eq. (5-17) with respect to x 
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(5-18) 
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Thus, it should not be surprising that given the shape of the measured static 
pressure distribution upstream of separation (i.e. 3 2 Cp/3x 2 > 0) that the 
decelerating flow was not in equilibrium. Since 36*/3x > 0 and 9Cf/9x < 0 
upstream of separation, all the terms on the rhs of eq. (5-18) are positive, 
thereby precluding 3g/3x from equaling zero, the condition for equilibrium. 

In decelerating flows, equilibrium can only be achieved when3 2 Cp/3x 2 < 0. 

Cutler and Johnston (ref. 6) have noted that in strongly adverse pressure 
gradient and detaching flows the concept of equilibrium is not too useful. 

The boundary layer loses profile and structural self-similarity when the inner 
layer cannot adjust fast enough to changes imposed on the outer layer by the 
adverse pressure gradient. 

A measure of the strength of the gradient can be obtained from a ratio of 
the boundary layer thickness, 6, and a characteristic length scale of the 
adverse gradient, U e /(3U e /3x). 


6 

R = — — ^ (5-19) 

U 9x 


This non-dimensional gradient strength, R, usually can be calculated from 
experimental data when detailed profile information is not available to evalu- 
ate the wake parameter, II. In the following section values of R have been 
calculated for the experiments of Simpson et al. (ref. 2) and Sandborn and Liu 
(ref. 49) and are compared to the non-dimensional gradient strengths 
calculated from data measured in the current study. 

5.4.2 Analysis of Data in the Separation Region 

In the absence of a salient edge, turbulent boundary layer separation 
occurs over a region rather than at a single streamwise location as occurs in 
laminar boundary layer separation. This process of turbulent boundary layer 
separation on a faired surface has been described by Sandborn and Kline 
(ref. 23) as follows: 
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. . separations of a turbulent shear layer do not usually begin in a 
steady or two-dimensional fashion on a faired surface. Instead they more 
often commence with intermittent streaks of backflow extremely near the 
surface. At the present time the onset of these streaks of backflow is 
not well correlated; indeed, the onset apparently depends on the nature 
of the fluctuations in the free stream. However, this is probably not of 
crucial importance, since these streaks of backflow are now known to be 
only an evidence of the existence of islands of hesitation which form a 
portion of the inherent flow model of the wall layers of the turbulent 
shear layer described by Kline and Runstadler (ref. 48). In favorable or 
mild adverse pressure gradients these islands of hesitation merely move 
more slowly than neighboring fluid layers; backflows are not observed. 
However, as the adverse pressure gradient is increased, these layers of 
very small inertia begin, relatively soon, to show upstream motions and 
to create local intermittent flow separation with streaks of backflow 
near the surface of a very small length scale. As the adverse pressure 
gradient is increased still further, these local backflows begin to 
accumulate slow moving fluid near the wall faster than it can be removed 
by the through flow. This leads to large transient three-dimensional 
stall patterns wherein the slow moving fluid is removed by intermittent 
adjustments in the over-all flow pattern and pressure distribution. 
Finally, upon still further increase in adverse pressure gradient, this 
process leads to a complete breakdown of the boundary-layer potential 
flow pattern and a fully established stall." 

The following terminology has been defined (refs. 6 and 8) to quantify 
the process of turbulent boundary layer separation. 

Detachment (D) - the locus of points where the limiting streamline of the 
time-averaged flow leaves the wall; a location of zero mean skin fric- 
t ion. 


Incipient Detachment (ID) - the point upstream of detachment having back- 
flow at the wall 1 percent of the time (y p - 0.99). 


Intermittent Transitory Detachment (ITD) - the point upstream of detach- 
ment having 20 percent instantaneous backflow of the wall (y p = 0.80). 

*U 


Transitory Detachment (TD) - the point in 
50 percent instantaneous backflow at the 


the region of detachment having 
wall (y = 0.50). This point 


is usually coincident with or very close to the detachment point (D). 


The locations of incipient detachment, intermittent transitory detach- 
ment, and transitory detachment for the current study occurred at x = 32, 
35.5, and 42 cm, respectively, as shown in figure 5-13. Thus, the axial 
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extent of the separation region between incipient detachment and detachment 
was 10 cm which was less than the 11 cm thick boundary layer at Station 8 
(x = 31.75 cm) just upstream of incipient detachment. Sandborn and Liu 
(ref. 49) who studied a turbulent separation bubble on a curved surface also 
achieved separation over one boundary layer thickness whereas Simpson et al 
separated the flow more gradually over 4.4 boundary layer thicknesses. Cutler 
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Perry and Fairlie (ref. 5) did not measure the forward flow fraction. 


To compare the forward flow fraction data in the separation region with 
data obtained by other investigators, y p has been plotted vs. a shape factor, 


h = (H 2 -l)/H 2 , in figure 5-14. The shaded region represents the range of 
turbulent boundary layer separation and reattachment data surveyed by Kline et 
al (ref. 50). All the turbulent separation cases surveyed were near equili- 
brium. The data obtained upstream of detachment in the current study lies 
outside and to the left of the envelope of the equilibrium separation data 
surveyed. Transitory detachment was achieved in the current study at h = 

0.55, significantly less than the 0.69 < h < 0.77 range of the surveyed data. 
This result could have occurred because of the nonequil ibr ium nature of the 
decelerating boundary layer in the current study. The strong adverse pressure 
gradient caused detachment to occur before the profile could adjust to an 
equilibrium shape. 

Correlations have been developed in the H-A plane, where A * 6*/6, to 
demarcate the regions of unseparated flow, intermittent separation, and fully 
developed separation, Sandborn (ref. 51) defined the line of demarcation 
between unseparated flow and intermittent separation in a turbulent boundary 
layer by the following equation 


1 

H = 1 + (5-20) 

1-A 

Sandborn and Kline (ref. 23) demonstrated that eq. (5-20) correlated all known 
intermittent separation data available at that time. Sandborn (ref. 51) also 
defined a line in the H-A plane to denote the location of laminar boundary 
layer separation. Sandborn and Liu (ref. 49) later showed that the laminar 
separation line could be generalized as a zero wall shear separation line 
since it also demarcated the regions of intermittent and fully developed 
turbulent se par at ion. 

The correlations of Sandborn, the present experimental data, and the data 
of Simpson (ref. 1), Perry and Fairlie (ref. 5), Smith et al (ref. 25), and 
Sandborn and Liu (ref. 49), are shown in figure 5-15. Equations (5-7) and 
(5-9), v;hich include the effects of normal static pressure gradients within 
the boundary layer, were used to determine the boundary layer integral 
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properties from the current experimental data. The Sandborn curves correlate 
the experimental data of both Simpson and the current experiment very well. 

In the current experiment, intermittent separation (ITD) commenced just 
upstream of Station 11 and separation (i.e. detachment (D>) occurred down- 
stream of Station 13. The locus of data points from the present experiment 
also passes through the range of all experimental separation data surveyed by 
Sandborn. The location of intermittent and fully developed separation in 
Simpson* s experiment were also predicted by the Sandborn correlation. 

The data of Smith et al and Perry and Fairlie generally follow the trends 
of Simpson’s data which lie to the right of the data obtained in the current 
study. The data of Sandborn and Liu show intermittent and fully developed 
separations occurring near the boundaries predicted by Sandborn* s correla- 
tions. However, their data lie significantly to the left and outside the 
range of data surveyed by Sandborn and Kline. 


Values of R, the non-dimensional strength of the deceleration gradient, 
were calculated at the point of incipient detachment and transitory detachment 


•F /■A-*'* fr- V* t 


experiment and for the studies uf Simpson et al and Sandborn 


and Liu. At incipient detachment, R, was calculated using eq. (5-19) to be 
0.022, 0.052 and 0.061, for the data of Simpson, the current study, and 
Sandborn and Liu, respectively. At detachment, the same relationship existed. 
R equaled 0.026, 0.033, and 0.097 for the data of Simpson, the current study, 
and Sandborn and Liu, respectively. Thus the deceleration gradient imposed by 
Sandborn and Liu was larger than the gradient imposed in the current study, 
which in turn was larger than Simpson’s deceleration gradient. These data 
tend to indicate that increasing R causes the locus of data points in a given 
experiment to move up and to the left in figure 5-15. A consequence of this 
movement is that 6 /6 at separation varies inversely with R. 


Kline, Bardina, and Strawn (ref. 50) have developed an improved correla- 
tion for detachment for any boundary layer having a velocity profile which 
satisfies the modified Coles* wall-wake profile (ref. 30). 


U 


+ 





+ C 


(5-21) 


where Il(x) is the wake amplitude. Equation (5-21) has been extremely effec- 
tive in correlating the vast majority of turbulent boundary layers. The 
correlation becomes inadequate only for those boundary layers in which the 
inner and outer layers are far f rom e quilibr ium (ref. 50). Kline et al. cast 
eq. (5-21) into the following form 
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(5-22) 


U 

U„ 


1 + V T In (— ) - Vg cos 
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where Vg (the nond imensional wake amplitude) and V T (the dimensionless 
shear velocity) are defined as 


V B = 2U X n(x)/KU oo 


(5-23) 


V T = VkU. 


(5-24) 


respectively. The shear velocity, U T , is defined to be positive in both 
attached and separated flows. 

Using eqs. (5-4) and (5-5) to define the momentum and displacement thick- 
nesses, eq. (5-22) can be integrated directly to yield 


J - »T ♦ °' 5 V B - 2V T 2 


1.5895 V T V B - .375 V g 2 


(5-25) 


g V T + 0.5 V B 


(5-26) 


Solving eq. (5-26) for the wake strength, Vg , and substituting the result 
into eq. (5-25) yields the correlating equation 


h V T V T ^ # 

- = 1.5 + 0.179 - 1 + 0.321 (5-27) 

A A A 


where 


h = 




(5-28) 


A 



(5-29) 
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For equilibrium boundary layers the dimensionless shear velocity, Vp, is 
given implicitly by 


V 


T 


1 - 2A 



+ In (k Re^ *) + 0.05 


(5-30) 


The correlation of Kline et al. (eq. (5-27)) is plotted in figure 5-16 
along with the data of Cutler and Johnston, Sandborn and Liu, and the present 
experiment. The shaded region represents the range of near-equilibrium 
separation data surveyed by Kline et al. The data obtained in the current 
study upstream of separation ( i .e . , Stations 6, 8, and 11) and in the upstream 
half of the separation bubble (i.e., Stations 13-15) lie near the Cutler and 
Johnston data above the Kline correlation for equilibrium turbulent boundary 
layers and near the boundary of the shaded region. The data of Sandborn and 
Liu, which was obtained in a strong adverse pressure gradient on a convex 
surface, lie significantly above the correlation line and outside the shaded 
region. Data obtained in the downstream half of the separation bubble (i.e., 
Stations 18, 20, 22, 24 and 26) and downstream of reattachment (i.e.. Stations 
28 and 31) lie on or near the Kline correlation. 


Deviation of data obtained by Sandborn and Liu, Cutler and Johnston, and 
at Stations 6-15 in the current study from the correlation line is not 
surprising since the Kline correlation is restricted to the special case of 
equilibrium turbulent boundary layers. As discussed in Section 5.4.1 above, 
the decelerating boundary layer in the present experiment which was near 
equilibrium at the test section inlet lost its equilibrium structure when it 
encountered the strong adverse pressure gradient. To estimate the effect of 
boundary layer nonequilibrium on the Kline correlation, eq. (5-27) was 
rederived for a boundary layer having no wake component (i.e., II = V fi = 0). 

For a zero-wake boundary layer eq. (5-27) becomes 


h * 2A 


(5-31) 


As shown in figure 5-16 the data from Stations 6-15 generally lie between the 
equilibrium correlation (eq. (5-27)) and the zero-wake correlation (eq. (5- 
31)). As the flow moves downstream the data approaches the equilibrium 
correlation as the strength of the deceleration gradient diminishes. Ifote 
also in figure 5-16 that the Sandborn-Kline correlation not only predicts the 
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beginning of intermittent separation between Stations 11 and 13 but the 
reattachment location between Stations 26 and 28 also. 

In the fully developed separation region Simpson et al. (ref. 1) demon- 
strated that a similarity existed in the forward flow fraction profiles by 
normalizing and plotting 


% ■ Yp u • 

mm 


1 - Yi 


y 

VS . — 

M 


U 


mm 


(5-32) 


where ^P TT was the minimum value of y. in the profile near the wall and M 
U min P U 

was the distance from the wall to the peak in the axial turbulence profile. 

The data in the present experiment were normalized in the same manner and very 
good profile similarity was obtained. However, better similarity was obtained 
when the profiles were normalized with M chosen as the distance from the wall 
to the peak in the transverse turbulence profile. The profiles normalized by 
the distance from the wall to the peak in the transverse turbulence profile 
are plotted in figure 5-17. Data scatter is very small, being well within the 
bounds of Simpson's data which are based on M chosen as the distance from the 
wall to the peak in the axial turbulence profile. 


5.5 Description of the Separation Bubble 


Downstream of separation a large separation bubble, defined as Zone III 
in figure 5-1, having a maximum thickness, h, of 17 cm was formed between the 
mean dividing streamline and the upper wall test surface. The flow within the 
bubble was characterized by vigorous recirculation and high turbulence levels. 
As shown in figure 5-18, Zone III can be partitioned into three transverse 
subzones: Zone IIIA, the backflow boundary layer in which U < 0 and 3U/3y < 

0; Zone IIIB, the backflow shear layer in which U < 0 and 3U/3y > 0; and Zone 
IIIC, the forward flow shear layer, in which U > 0 and 3U/3y > 0. The 
backflow in Zones IIIA and IIIB was strong and steady. Mean backflow 
velocities exceeding 30 percent of the local freestream velocity have been 
measured within these subzones. Steady backflow (i.e., *Yp^ = was meas ured 


over a large portion of zones IIIA and IIIB. As indicated in figures 4-13, 
5-13, and Table II, the region throughout which the forward flow fraction was 
zero extended along the test surface from Station 18 (x = 67.3 cm) to Station 
24 (x = 90.2 cm). The maximum thickness of the steady backflow region was 
2.5 cm at Station 20 (x = 74.9 cm). Measurements of Yp = 0 in a separated 

boundary layer have not been reported previously (ref. 24). 
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As a result of the strong steady backflow within the separation bubble a 
well developed thin backflow boundary layer (Zone IIIA) grew along the test 
plate. The boundary layer, which had both laminar and turbulent characteris- 
tics, will be described in more detail in Section 5.5.1 below. 

Zone IIIB, the backflow shear layer, was approximately an order of magni- 
tude thicker than the backflow boundary layer. The backflow shear layer con- 
sisted of an inner layer having a relatively flat velocity profile (Zone 
IIIBl), which connected the viscous near wall region with the outer shear 
flow, and the outer backflow region (Zone IIIB2) which was part of the outer 
shear flow. Zones IIIA, IIIBl, and IIIB2 are components of the three-layer 
backflow model proposed by Simpson (ref. 1). Zone IIIB was characterized by a 
nearly constant turbulence level approximating 8 percent of the inlet free- 
stream velocity as shown in figure 4-18. 

Zone IIIC, the forward flow shear layer, defined as the region between 
the locus of U=0 points and the dividing streamline (^ = 0) contained higher 
normal velocity gradients, 3U/9y, and turbulence levels than Zone IIIB. In 
fact, the largest values of 9U/3y within the reversed flow region between 
Stations 13 and 26 occurred at the outer edge of Zone IIIC. As discussed 
below in Section 5.8, "Shear Stress Correlations" the largest values of the 
streamwise shear stress, aB , also occurred at the outer edge of Zone IIIC near 
the dividing streamline. As shown in figure 5-19, Zone IIIC was slightly 
thinner than the combined thicknesses of the backflow zones (Zone IIIA + Zone 
IIIB) at all stations. The thicknesses of the backflow zones and the forward 
flow shear layer (Zone IIIC) increased linearly over the upstream half of the 
separation bubble. The origin of the linear growth region was the point of 
transitory detachment (i.e.. Yd = 0.50 at the wall). 

5.5.1 Description of the Backflow Boundary Layer 

Although the existence of a steady backflow boundary layer was not 
anticipated at the outset of the test program, the measurement matrix was 
sufficiently dense to permit a limited analysis of the structure of the back- 
flow boundary layer. Backflow velocity profiles are plotted in figure 5-20. 
The profiles were calculated from backward- facing pitot data. The limited 
number of LV measurements made within the backflow boundary layer indicated 
that the backflow was steady (y D = 0) between Stations 23 and 18 (i.e. over 

the central 35 percent of the axial extent of the bubble) and only slightly 

intermittent (y^ < 0.12) at all other locations. LV velocity measurements 

U 

were in good agreement with the pitot data. 
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The boundary layer appeared to originate near the reattachment point 
(xr - 103 cm) just downstream of Station 27. The boundary layer grew until 
it reached a maximum thickness of approximately 0.76 cm at Station 15 near 
Zone II, the intermittent separation zone. The pitot data obtained at Station 
15 are not presented in figure 5-20 because values of as high as 25 per- 
cent measured near the wall at Station 15 could make the data unreliable. 


With the understanding that the accuracy would be less than desired 
because of the small number of data points obtained within the backflow 
boundary layer, boundary layer integral properties were calculated for each 
profile presented in figure 5-20. Calculated values of the momentum thickness 
Reynolds number, ReQ , the shape factor, H, and other backflow boundary layer 
parameters are listed in Table VII. Boundary layer parameters determined by 
Coles (ref. 30) for a zero pressure gradient turbulent boundary layer are also 
listed in Table VII for comparison. Although the maximum calculated value of 
ReQ , which was 184, was much lower than the value of 425 required for 
transition in low-f reestream-turbulence , zero-pressure-gradient boundary 
layers (Table 2 in Appendix A of ref. 30), the calculated value of H was 
transitional. The shape factor decreased from 1.98 at Station 25 to 1.48 at 
Station 16. The values are between the shape factor of 2.59 for a laminar 
boundary layer having a Blasius velocity profile (ref. 35) and the value of 
1.47 given by Coles (ref. 30) for a turbulent boundary layer for which Reg 
equals 855. 

To gain additional insight into the nature of the backflow boundary layer 

the local axial turbulence levels, Vv?7u, measured within Zone IIIA with the 
LV system were plotted in figure 5-21. The turbulence level was uniformly 
high across the boundary layer, varying from approximately 50 percent near the 
wall to 40 percent at the 11 freestream." The turbulence was dominated by low 
frequency large scale eddies (see Section 4.6). The turbulence profile 
measured by Klebanoff (ref. 31) for a flat plate turbulent boundary layer 
having a low level of freestream turbulence has been shown for comparison. 

The local turbulence levels in the Klebanoff profile are less than 10 percent 
over the outer 90 percent of the boundary layer. Because of the low free- 
stream turbulence level, the Klebanoff boundary layer is dominated by 
turbulence which is generated near the wall and diffused outward. If such 
turbulent transport exists in the backflow boundary layer it is masked by the 
large scale turbulent eddies. 

The turbulence levels within the backflow boundary layer exceed the 
turbulent boundary layer values measured by Klebanoff at all points in the 
profile. In this respect, the backflow boundary layer is clearly turbulent. 
However, it could be argued that the LV turbulence measurement resulted from 
averaging highly unsteady laminar flow at the wall which is driven by the 
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large scale turbulent eddies in the "freestream." Therefore, the backflow 
boundary layer can only be referred to as pseudo-turbulent until its exact 
nature can be determined from more detailed measurements. 

Simpson (ref. 26) has proposed that the backflow velocity profile in 
Zones IIIA and IIIB has a universal shape defined by the equations 



where N is the value of y at the transverse location of the maximum backflow 
velocity, U^. The present LV data are compared to the Simpson correlation 
in figure 5-22. Setting A = 0.3 as suggested by Simpson, eq. (5-33) 
correlated only the profile data at Station 15 just downstream of detachment. 
Setting A = 0.2 provided a better correlation for data in Zone III A (i.e., 
y/N < 1). Considering the scatter in the profiles in Zone IIIB when plotted 
in y/N vs. U/U N coordinates, eq. (5-33) could not correlate the Zone IIIB 
data for any value of A. Simpson also noted difficulties in using eq. (5—33) 
to correlate Zone IIIB profiles. 

The flow at Station 15 most closely approximated the flow conditions used 
in the test cases surveyed by Simpson in developing eq. (5-33). The velocity 

at Station 15 was intermittent having a minimum value of y p of 0.08. All 

U 

other profiles presented, except the Station 16 data, had steady backflow 
(i.e., Yp^= 0) throughout the profile. None of the profiles surveyed by 

Simpson had steady backflow at any point in the profile. 

Because of the highly turbulent nature of the backflow boundary layer, an 
attempt was made to plot the backflow data in U + vs. y + law-of-the-wall 
coordinates. Simpson had suggested (ref. 1) that such a scaling was not 
possible in the backflow boundary layer because the law-of-the-wall length 
scale v/U T varies inversely with its velocity scale U T which differs from 
the behavior of N and which both increased in the downstream direction. 

In the present data, the behavior noted by Simpson occurred over a short 

region just downstream of detachment (Stations 13-15) where > 0 as shown 

U 

m figure 5-14. In the steady backflow region (Stations 18-24), however, N 
increased in the direction of the backflow while remained relatively 
constant (see Table VII) which is exactly the behavior attributed to boundary 
layers which can be modeled in law-of-the-wall coordinates. 
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Backflow boundary layer data plotted in law-o f-the-wall coordinates are 
presented in figure 5-23. The data are characterized by a short log-linear 
region (11 < y + < 60) and a negative wake. Data scatter is minimal. The 
absence of a wake component in the backflow boundary layer is not surprising. 
Coles has plotted the strength of the wake component in equilibrium turbulent 
flow vs. Reg in figure 12 of reference 30. Coles's plot shows that the 
magnitude of the wake component is zero at Reg - 425 and increases to its 
equilibrium value at Reg - 5000. The maximum value of Re 0 in the backflow 
boundary layer was 184. Coles also noted in figure 15 of reference 30 that 
increasing freestream turbulence decreased the strength of the wake component. 
Kline et al . (ref. 52) showed that the wake component in an equilibrium 
turbulent boundary layer, which was normal at a freestream turbulence level of 
1.8 percent, decreased to zero when the freestream turbulence level increased 
to 4.5 percent. As shown above in figure 5-21 the backflow boundary layer had 
a "freestream" turbulence level of 40 percent. 

The friction velocity, U T , and the skin friction coefficient, C f = -2 

(U /U r) 2 , were determined from the 1 aw^o f-the-wall plots. The skin friction 
Tret # , 

values, shown in figure 5-23, were approximately an order of magnitude lower 

than the values measured in the attached boundary layer upstream of separa- 
tion. However, when localized by the local maximum backflow velocity, Uy, 
rather than the inlet reference velocity, the skin friction coefficient, C f ' = 
2 (U t /U n ) 2 , was quite large. C f ' has been plotted vs. Re x i, in figure 5-24 
where Re x i = U R x'/v and x' = x R -x where x R , the axial location of reattach- 
ment, equaled 103 cm. 

Also plotted in figure 5-24 are values of Cf' estimated from the normal 
velocity gradient at the wall at each station within Zone III using the 
following equation 


C 


f 


V If ly=0 Re 9 ^ 9U N 

p U n 2 /2 2 9 9x 


(5-34) 


The second term on the rhs of eq. (5-34) is a correction term to be used in 
flows having axial pressure gradients (ref. 53). The correction term did not 
exceed 12.5 percent at any station within Zone III. 

Cf' calculated from eq. (5-34) decreased from 0.022 at Re x ? = 2. 1 x ICf* 
near reattachment to 0.0074 at Re x i = 1.1 x l(P near detachment. Very good 
agreement was obtained between values of ' calculated using eq. (5-34) and 
the values determined from the law of the wall plots. Values of Cf* 
correlated reasonably well with the equation 
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C' f = 3.29 Re x , 


- 1/2 


(5-35) 


. - 1/2 

which has the same Re x » ' dependence but which is 400 percent higher than 
zero pressure gradient laminar skin friction coefficients at the same Reynolds 
number . 

Two possible explanations for the high skin friction coefficients 
measured in the backflow boundary layer could be the favorable pressure 
gradient and the high turbulence levels impressed upon the backflow boundary 
layer by the external flow field. The effect of the favorable pressure 
gradient on the shear at the wall can be estimated using the Pohlhausen 
parameter, A, (ref. 35, p. 245) 


where 


w 


= y 


3U 

3y 



A 


dp 

dx 1 yU N 


(5-36 ) 


The effect of the pressure gradient on the skin friction can be estimated from 
eq. (5-36) as 


C 


i 



C' (A=0) 
t A 



(5-37) 


The maximum value of A within the backflow boundary layer occurred at Station 
23 (x = 86.4 cm) where A = 13.9. Thus, the effect of the favorable pressure 
gradient could account for as much as 100 percent of the 400 percent increase 
in skin friction at some points in the backflow boundary layer. 

The effect of freestream turbulence in reducing the transition Reynolds 
number and increasing the skin friction coefficient is well known but not well 
quantified. Dyban et. al. (ref. 54) have measured a 56 percent increase in 
skin friction in a laminar boundary layer having a 25 percent freestream 
turbulence. The increase in the backflow boundary layer skin friction 
coefficients due to the 40 percent 11 freestream" turbulence within the 
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separation bubble would be expected to be large but cannot be estimated at 
this time. 

5.6 Description of the Reattachment Process 

Re attachment of the separated shear layer is an inherently unsteady 
process which occurs within Zone IV, the Reattachment Zone. Smoke flow 
visualization movies showed that the Reattachment Zone was characterized by 
random streamline flapping (shown schematically in figure 4-6) with large 
scale eddies being convected alternately upstream and downstream from the 
impingement point on the test surface. The reattaching flow, impinging 
steeply onto the test surface, caused a static pressure peak which drove the 
strong backflow within the separation bubble described above. 


The flow within Zone IV was characterized by frequent flow reversals and 
high turbulence levels. The flow reversals have been quantified by the 
forward flow fraction, which is the fraction of the time the flow is moving in 
the downstream direction. Values of the forward flow fraction, y p , near the 


test surface (y = 0.10 cm) at the four LV measurement stations within the 
reattachment zone have been plotted in figure 5-25 versus x*, a non- 
dimensional axial coordinate 


* 


x 


x " X R 
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x 6* 


max 


(5-38) 


where xr is the reattachment location defined as the point on the test 

surface where Cr = 0 and Xr* is the axial coordinate at the location of the 

° max 

maximum separation bubble thickness. The definition of the non-dimensional 
coordinate, x*, has been chosen to permit direct comparison of data in the 
reattachment zone to data obtained in backward facing step flows where the x- 
coordinate has its origin at the step and x* = (x-xr)/xr (ref. 17). 

Note in figure 5-25 that y p rises from 0.02 at x = -0.44 to 1.0 at x = 

U 

+ 0.40. If the reattachment zone is defined as the region between 0.01 < y p < 

“ F U “ 

0.99 and x - Xr* is approximately 1/2 the bubble length, the reattachment 
0 max 

zone can be calculated to be 40 percent as long as the separation bubble. A 
correlation developed by Westphal et al. (ref. 17) for reattaching flow down- 
stream of a backward facing step shows excellent agreement with the forward- 
flow fraction data obtained in the present study. 
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The skin friction coefficient measured within the reattachment zone has 
been plotted vs. the forward flow fraction in figure 5-26. The skin friction 
increased monotonically with forward flow fraction throughout the reattachment 
zone. Also plotted in figure 5-26 are data from two cases investigated by 
Westphal et al. (ref. 17) for reattaching flow behind a backward facing step. 
The backward facing step data showed similar trends to the data in the present 
study. In particular, for both cases studied by Westphal et al . and for the 
present data the location of Cr = 0 occurred near the location at which y p = 

0.50. U 

Downstream of the reattachment location the reattached flow recovers the 
structure of a flat-plate turbulent boundary layer. This process occurs in 
two stages. In Zone V, a new boundary layer develops along the test surface 
and quickly achieves equilibrium conditions. The relaxation of the outer 
shear layer in Zone VI occurs much more slowly, however. 

Streamwise velocity and turbulence profiles downstream of reattachment 
are plotted in figure 5-27. The two-layer structure of the reattached bound- 
ary layer can be seen quite clearly in the figure as the new boundary layer 
grows rapidly toward equilibrium under the relaxing outer shear layer. The 
effective "f reestream” turbulence level impressed upon the developing boundary 
layer was quite high. The turbulence level was approximately 10-14 percent of 
U re f anc * on the order of 25 percent of the streamwise velocity at the edge of 
the new boundary layer at all measurement stations downstream of reattachment. 
Recovery of the two-layer profile to a flat-plate turbulent boundary layer 
profile was not complete at x = 148 cm, the final measurement station within 
the test section. 

To study the structure of the recovering flat-plate turbulent boundary 
layer in more detail, axial velocity profiles at these measurement stations 
downstream of reattachment have been plotted in law^of-the wall coordinates in 
figure 5-28. At Station 28 (x* = 0.09) at reattachment the boundary layer 
structure had no log-linear region. The profile was dominated by the wake- 
like structure of the reattaching outer shear layer. At Station 30 (x* - 
0.36) a short log-linear region had become established at y + < 100 but the 
profile dipped below the standard log-linear line before reverting to a 
strong, albeit somewhat diminished from Station 28, wake component at the edge 
of the boundary layer. At Station 34 (x* = 1.36) the dip below the log- 
linear line persisted and the wake component had diminished to the extent that 
the profile had a negative wake component. 

The dashed lines in figure 5-28 represent data obtained by Chandrsuda and 
Bradshaw (ref. 18) in a reattached fully developed turbulent mixing layer 
downstream of a backward facing step. The profiles obtained by Chandrsuda and 
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Bradshaw show the same characteristics as the profiles obtained in the present 
study. Chandrsuda and Bradshaw explained that the dip below the log-linear 
line occurred because the standard logarithmic law is based upon the assump- 
tion that the length scale of the turbulence is proportional to y, while at, 
and just downstream of, reattachment the length scale will be nearly constant, 
as it was in the separated shear layer, except very near the surface. 

Bradshaw and Wong (ref. 11) have shown that the dip persists for 50 step 
heights (x* = 7.5) downstream before equilibrium is reached in the 
turbulence energy exchange between the small length scale turbulence being 
generated near the wall and the large length scale turbulence being fed into 
the developing boundary layer from the reattached mixing layer. In the 
present study x* = 7.5 would correspond to an axial location approximately 
160 cm (i.e., 3.75 bubble lengths) downstream of reattachment. 

The displacement thickness, 6*, the momentum thickness, 0, and the 
shape factor, H, have been calculated throughout the reattachment region and 
are plotted in figure 5-29. The open circles represent data calculated using 
eqs. (5-4) and (5-5) based on classical boundary layer theory and the solid 
circles represent pressure-based integral properties calculated using eqs. 
(5-8) and (5-10) (see Section 5.3.3 "Boundary Layer Integral Properties"). 

Also plotted in figure 5-29 are integral property data obtained by Driver and 
Seegmiller (ref. 14) in flow reattaching behind a backward facing step within 
an adverse pressure gradient. The displacement and momentum thicknesses have 
been non-dimensionalized by h which equaled the maximum separation bubble 
thickness of 17 cm in the present study and the step height in the Driver and 
Seegmiller study. 

The data obtained by Driver and Seegmiller for their reference case in a 
nondiverging tunnel having zero pressure gradient in the reattachment zone are 
plotted as squares in figure 5-29. The trends shown in their data, namely the 
raonotonic decrease in displacement and shape factor, and the generally 
increasing value of momentum thickness with axial distance are similar to the 
trends observed by Kim et al. (ref. 13) and Pronchick and Kline (ref. 16) for 
zero pressure gradient backward facing step flows. For flows reattaching in 
an adverse pressure gradient (a = 6 deg) Driver and Seegmiller have shown that 
the displacement thickness decreases less rapidly and the momentum thickness 
increases more rapidly with axial distance than the zero pressure gradient 
case while the behavior of the shape factor is relatively unaffected by the 
pressure gradient. Conversely, in the present study the separated shear layer 
reattached within a converging channel (a = -38 deg) having a very favorable 
pressure gradient downstream of reattachment. Thus, the rapidly decreasing 
displacement and momentum thicknesses downstream of reattachment in the 
present study follow the trends identified by the Driver and Seegmiller study 
of backward facing step flows. The behavior of the shape factor calculated 
from classical boundary layer theory was uninfluenced by the pressure 
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gradient (ref. 14 data), whereas the more accurate pressure-based shape factor 
calculated in the present study displayed a more rapid decrease with axial 
distance within the favorable pressure gradient. 


5.7 Streamwise Turbulent Shear Stress Measurements 

As noted by Cutler and Johnston (ref. 6), when the mean streamlines are 
inclined at a significant angle to the test surface, turbulent shear stress 
data are more meaningful when presented in an orthogonal coordinate system 
comprised of the mean streamline direction and the normal to the streamline 
than in a coordinate system orthogonal to the test surface. This is because 
flow history effects are important in boundary-layer flows (except near the 
wall) and the streamline orthogonal coordinate system can track the flow 
development along a given streamline. Cutler and Johnston concluded that the 
behavior of a turbulence quantity is more likely to be modeled in "universal" 
form if streamline orthogonal coordinates are used. In the current study, 
flow within the boundary layer was inclined more than 30 degrees relative to 
the test surface near separation and reattachment, significantly more than the 
10 degree inclinations reported by CJutLer and Johnston. 

The Reynolds stress, uv, data in the wall orthogonal (x,y) coordinate 
system can be rotated into the streamwise orthogonal (s,n) coordinate system 
using the following relation (ref. 6) 


ag = uv (cos 2 <(> - sin 2 <j>) + (v^ - u^) sin$ cos<J> 


(5-39) 


where ag is the streamwise turbulent shear stress and <J> = tan ^ (V/U), the 
angle of the local mean-flow streamline to the wall. The wall orthogonal 
shear stress data presented in figure 4-21 have been rotated into streamwise 
orthogonal coordinates. The resulting values of otg are tabulated in Table III 
and are plotted in figure 5-30. 

Trends are much more evident in the streamwise shear stress profiles than 
in the uv profiles plotted in figure 4-21. At Station 6 upstream of incipient 
detachment the shear stress peaks near the wall. At Station 11 as the bound- 
ary layer approaches detachment the peak increases in amplitude and moves away 
from the wall. Downstream of detachment, between Stations 13 and 18, the peak 
grows much more rapidly and is located near the edge of the backflow region 
(i.e., locus of U = 0 points) indicating strong turbulence generation in this 
region. Between Stations 22 and 28 aft of the maximum bubble thickness the 
peaks broaden and decay and move farther away from the U = 0 line into the 


119 


wake of the bubble. This decaying trend is reversed as the shear stress 
begins to grow again as the flow is accelerated into the exit duct at Stations 
30 and 31. No data has been presented near reattacheraent where highly 
unsteady flow and steep flow angles make the rotated data of questionable 
value . 
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Figure 5-1 Separation Bubble Flow Zones 




Mean Velocity Measurements (z = 0) 



AXIAL LOCATION, x (cm) 





Figure 5-5 Normalized Axial Velocity Profiles Overlayed to Demonstrate 
Invariant Outer Boundary Layer Profile 
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Figure 5-7a Skin Friction Coefficient on Test Surface (z=0) 


























Figure 5-9 Normalized Total Pressure Profiles Showing Spanwise Uniformity of Flow 
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Figure 5-10 Spanwise Variation in Momentum Thickness 
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Figure 5-15 Comparison of Mean Velocity Profile Parameters at Separation with 
Sandborn Correlation 
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Figure 5-16. Comparison of Experimental Data to Correlation of Kline et al (ref. 50) for 
Attached and Separated Flows 
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Figure 5-17 Normalized Forward Flow Fraction Profiles in the Separation Region 
(Note Log-Linear Abscissa) 
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ZONE III 


Figure 5-18 Definition of Zone III (Separation Bubble) 
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Figure 5-19 Growth of Forward-Flow and Backflow Regions within Separation Bubble 


85 - 2 - 23-3 


143 






I 


STATION 

NO. 

X 

(cm) 

N 

(cm) 

o 

16 

59.7 

0.76 

o 

18 

67.3 

0.76 

L 

20 

74.9 

0.50 

A 

22 

826 

0.50 

D 

24 

90. 2 

0.15 



Figure 5-21 Local Axial Turbulence Levels in Backflow Boundary Layer 
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Figure 5-23 Backflow Velocity Prefiles Plotted In Law-of-the-Wall Coordinates 
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Figure 5-24 Skin Friction Coefficient in Backflow Boundary Layer 


85 - 2 - 23-4 







Figure 5-25 Forward Flow Fraction at y = 0.10 cm in Reattachment Region 
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Figure 5-26 Skin Fri 
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Figure 5-30 Nondimensional Streamwise Shear Stress Profiles Determined from Laser 
Velocimeter Data (z = 0) 




CHAPTER 6 


CONCLUSIONS 


Overall conclusions and observations regarding the nature of the 
separated flowfield investigated in the current study are summarized here . 
Following these general remarks are more specific comments regarding flow in 
the various zones which comprised the overall flowfield. 

The wind tunnel test configuration employed here produced a large scale 
(55 cm. axial extent, 17 cm. maximum thickness) closed separation bubble on a 
flat plate test surface with an incident near-equilibrium turbulent boundary 
layer. Based on flow visualization and spanwise flowfield measurements, the 
separated flowfield is concluded to be a good approximation to an idealized 
two-dimensional flow with separation and reattachment and hence a suitable 
reference case for assessing computational procedures used to predict two- 
dimensional separated flow phenomena. With regard to flowfield steadiness, an 
outer inviscid zone was defined in which velocity fluctuation levels were low 
(approximately two percent of the reference inlet velocity). From this low 
level it was concluded that the overall wind tunnel test section flow was 
steady although high levels of unsteadiness were obtained within certain 
viscously dominated flow zones. Lack of periodicity in the unsteady velocity 
signal suggests that the intermittent separation and unsteady reattachment 
process were not linked by a discrete frequency feedback loop mechanism. 

As discussed in more detail below, the separated flow region at the 
forward end of the bubble was found to be generally similar to that obtained 
in previous open separation studies (turbulent separation without reattach- 
ment). While this could be expected, an unanticipated result was that the 
reattaching flow at the aft end of the bubble bore many similarities to 
reattaching flows downstream of backward facing steps. In addition to these 
results regarding the forward and aft end of the separation bubble, important 
features were identified in the near wall and outer flow regions. A reverse 
flowing boundary layer was identified at the base of the bubble and the flow 
around the bubble was found to be basically inviscid and rotational. Con- 
gruence of velocity profiles was obtained in this outer flow region supporting 
a similar observation in a previous separated flow study with a different test 
geometry. 

For the purpose of providing a more detailed analysis of the various 
flowfield phenomena encountered here it was found that the overall flow could 
be conveniently divided into eight flow zones. Observations and conclusions 
for these various zones are given below: 


155 



Zone I - Attached Flow Region - Flow in this incident boundary layer 
region consisted of a near equilibrium turbulent boundary layer which lost its 
equilibrium structure when it encountered a strong adverse pressure gradient. 

Zone II - Intermittent Separation Region - Flow in this region exper- 
ienced an increasingly adverse pressure gradient causing detachment to occur 
within a distance equal to one incident boundary layer thickness downstream 
from incipient detachment. Detachment occurred at a shape factor, H, of 2.2 
in line with established separation criterion, although the boundary layer did 
not track the (H-l)/H vs. 6*/<$ (i.e. displacement thickness/boundary layer 
thickness) curve established by Kline et. al. (ref. 50) for equilibrium 
separation. 

Zone III - Separation Bubble - This zone, defined as the region between 
the mean dividing streamline and the test surface, consisted of three sub- 
regions: a backflow boundary layer (Zone IIIA) at the base of the bubble, a 

backflow shear layer (Zone IIIB), and a forward flow shear layer (Zone IIIC) . 
The backflow shear layer which is approximately an order of magnitude thicker 
than Zone IIIA, consists of an inner layer (Zone IIIBl), having a relatively 
flat velocity profile and a constant turbulence level approximating 8 percent 
of the inlet freestream velocity, which connected to viscous near wall region 
to the outer backflow region (Zone IIIB2) which was part of the outer shear 
flow. The forward flow shear layer. Zone IIIC, defined as the locus of U = 0 
points and the dividing streamline, contained higher velocity gradients and 
turbulence levels than Zone IIIB. Zone IIIC was slightly thinner than the 
combined thicknesses of the backflow zones (Zones IIIA and Zone IIIB). The 
thicknesses of the backflow zones and the forward flow shear layer (Zone IIIC) 
increased linearly over the upstream half of the bubble from the point of 
transitory detachment. 

Although the three-layer backflow structure defined by Zones IIIA, IIIBl, 
and IIIB2 agrees with the model proposed by Simpson (ref. 1), the backflow 
velocity profiles only matched his "universal 11 backflow profile (ref. 26) in 
the intermittent flow region just downstream of detachment. This lack of 
correlation can be attributed to the vigorous backflow within the closed 
separation bubble driven by the static pressure peak at reattachraent . A large 
region of steady backflow, wherein the LV data indicated flow in the upstream 
direction 100 percent of the time, existed near the test surface over the 
central 35 percent of the axial extent of the bubble. 

Significant normal static pressure gradients existed within the boundary 
layer in the vicinity of the separation bubble due to streamline curvature 
induced by bubble blockage. Changes in the blockage, measured by the 
displacement thickness, which exceeded 20 cm at the maximum bubble thickness, 
caused convex streamline curvature at separation and at reattachment and 
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concave curvature over the central part of the bubble. The resulting static 
pressure field was calculated from measured total pressures and LV velocity 
data. 


The inclusion of the normal static pressure changes, which were as large 
as 50 percent of the dynamic head between the test surface and the lower wall 
near the location of the maximum bubble thickness, was found to be very 
important. Integral properties were calculated three ways: (1) the classical 

definition based on the boundary layer edge velocity, (2) a pressure-based 
definition which calculated the retardation of the boundary layer relative to 
an equivalent inviscid field determined from the measured upstream total 
pressure and the local derived (from measurements) static pressure, and (3) 
the true definition (refs. 38 and 39) which was equivalent to the classical 
definition with the addition of correction terms to account for the normal 
static pressure gradient effects. These effects were especially evident at 
the location of the maximum bubble thickness where the shape factor, H, was 
calculated to be 16.1, 11.8, and 7.5 based upon the true, pressure-based, and 
classical definitions, respectively. 

Because of the large streamline inclination angles trends in the turbu- 
lent shear stress measurements were not evident until the measurements 
obtained in a coordinate system orthogonal to the test surface were rotated 
into a coordinate system orthogonal to the local streamline direction. High 
values of the resulting streainwise shear stress, indicative of strong turbu- 
lence generation, were measured near the edge of the backflow region (i.e., 
near the locus of U = 0 points) over the forward half of the bubble. Down- 
stream of the maximum bubble thickness, the peaks broadened, decayed, and 
moved away from the U = 0 line into the wake of the bubble. 

Zone IIIA - Backflow Boundary Layer - The region at the base of the 
separation bubble is a highly unsteady laminar boundary layer having pseudo- 
turbulent characteristics. The backflow velocity profiles had a limited log- 
linear region with no wake component when plotted in law-of-the-wall 
coordinates. Skin friction coefficients derived from the law-of-the-wall 
plots (and corroborated by shear stress estimated from the normal velocity 
gradient at the wall) were approximately 100 percent higher than skin friction 
coefficients previously reported for intermittent backflow under open separa- 
tion bubbles. The measured skin friction coefficient showed the Re x 7 1/2 
dependence (where x 1 is distance upstream from reattachment) typical of 
laminar boundary layers but was 400 percent higher than values anticipated for 
flat-plate zero-pressure-gradient boundary layers. Some of the increased 
shear stress is due to the favorable pressure gradient experienced by the 
backflow boundary layer as it moves upstream from the static pressure peak at 
reattachment. It is conjectured, however, that the majority of the increase 
can be attributed to the effect of low frequency (i.e., length scales on the 
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order of the bubble size) turbulence in the bubble which is typically 40 
percent of the maximum backflow velocity. 

Zone IV - Reattachment Zone - In this zone the reattaching flow impinges 
steeply onto the test surface causing a static pressure peak which drives the 
strong backflow within the separation bubble. The highly unsteady flow in 
this zone is characterized by random low frequency streamline flapping with 
large scale eddies being convected alternately upstream and downstream from 
the impingement point on the test surface. The forward flow fraction measured 
near the test surface, when plotted vs. non-dimensional distance from 
reattachment, matches the correlation of Westphal et al. (ref. 17) for 
reattaching backward-f ac ing step flows. The appropriate length scale for 
closed separation bubble data is the distance from the maximum bubble 

thickness, x** , to reattachment. In addition, skin friction correlated with 

°max 

forward flow fraction in Zone IV in the present study in the same manner as 
the backward- fac ing step data measured by Westphal et. al. (ref. 17). 

Zone V - Redeveloping Near Wall Flow - This zone encompasses the near 
wall region downstream of reattachment where a new boundary layer begins to 
develop along the test surface. The boundary layer, which has high turbulence 
levels impressed upon it from the bubble wake, quickly develops a log-linear 
velocity profile near the wall and begins to merge with Zone VI, the relaxing 
outer shear layer, to obtain the wake-like structure of an equilibrium tur- 
bulent boundary layer. 

Zone VI - Relaxing Outer Shear Layer - The free shear layer in this zone 

originated near the mean dividing streamline surrounding the bubble and was 

shed into the bubble wake. It contains relatively large length-scale, high 

intensity turbulence. This relaxing outer shear layer merges with the 

redeveloping near wall flow (Zone V) and begins an exchange of energy which 

eventually would lead to the development of an equilibrium turbulent profile 

given a sufficiently long test section. In the current study, law^o f-the-wall 

plots of the merged boundary layer (Zones V and VI) contained negative wake 

components characteristics of reattaching flows. The data correlated well 

with the profiles of Chandrsuda and Bradshaw (ref. 18) when x* * - x R was 

w max 

used to non-dimensionalize the axial distance downstream from reattachment 
where x^ was the axial location at reattachment. It is estimated that 3.75 
bubble lengths downstream of reattachment, which exceeded the test section 
length, would have been required to eliminate the negative wake. 

Zone VII - Invariant Outer Flowfield - The flow in this zone is inviscid 
but rotational. Velocity profiles in this region are congruent over the axial 
extent of the bubble when static pressure variations associated with stream- 
line curvature are accounted for. In effect, the outer 80 percent of the 
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incident boundary layer is convected unaltered over the blockage associated 
with the inner viscously dominated bubble recirculation region. The identifi- 
cation of this feature of the velocity profile development appears to have 
important implications for the numerical modeling of the flow. 
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0 0118 
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0 187 
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0 37 4 
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0 74 8 
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0.218 

0 02 6 

0.241 

0 .014 
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0 014 

0 .364 

0.014 

0.409 
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0 

59 3 

0 
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0 
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0 
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1 

000 
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1 
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1 
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.000 

1 
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l 

.000 
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0.0160 
-0 . 01 10 
0 04 9 0 
0.0 300 


0 0140 
0 0 500 
-00120 
-00110 
-0 0020 


0 027 1 
0 0399 
0 0 4 43 
0.02 95 
0.0 25 1 
0 0157 
0 0113 
0.00 63 
0 0041 
0 0022 
0 00 1 4 
0.0013 
0.0014 


-0 1364 
-0 1585 
-0 .2 0 0 3 
-0 1599 
-0 1274 
-0 1295 
- 0.1056 
- 0.0225 
-0 0444 
-00149 
-0 .0 0 2 5 
-0 0061 
0 .0007 
-0 .0 0 30 



TABLE III REYNOLDS STRESS DATA MEASURED WITH LASER VELOCIMETER 








































































































OrciGSiwiL PAGE io 

OF POOR QUALITY 


STATION NO. Ci 3 


Cl 3 


«= 4 B 3 
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1 r.m 
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\ 

j=, 

V v Z /U, e( 
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0 5 113 




3 810 
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10.160 
1 2 70 0 
15.240 
1 7 7 8 0 
20.320 
22 160 


27 94 0 
30.480 


0 0 2 0 
0.041 
0 0 6 1 
0 0 8 2 
G 1 02 
0 122 
0.163 
0 2 0 4 
0 2 45 
0 2 85 
0.326 
0 4 C 6 
0 4 8 9 
0 5 7 1 
0.653 
0 7 3 4 
0 8 16 
0 8 9 7 
0 9 79 


0 .068 
0 07 ? 
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0 112 
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0 113 
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-0.0 82 0 
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0 . 2280 

0.023 
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0 .017 
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0.0069 
0.007 1 
0 0 0 8 7 
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0 0 5 6 0 
0.0 2 6 8 


0 . 0023 
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0 0 5 4 1 
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0.0618 
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0 . 007 6 
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0.0189 
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0 0 3 3 4 


0 0106 
0 0 6 5 6 
0.0311 
0 0 07 ? 
0.003 6 
0.0 02 4 


TABLE III REYNOLDS STRESS DATA MEASURED WITH LASER VELOCIMETER (CONT.) 









































































































































































STATION NO 


( = jjl 2 
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700 

15 

.240 

17 
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20 
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-0 06 7 
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0.085 
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0 0 8 4 
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-0 068 
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0 0 3 4 
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-0.081 

0 3 03 

0.075 
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0 .347 

0 126 
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0 0 8 8 

-0.142 

0 433 

0 217 
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0.520 
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0 52 3 
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0.0456 
0.0 38 8 
0 0172 
0 0 2 6 1 
0.0187 
0 007 ? 
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0.0017 


a0/U fc fe j 


-00749 

-01184 

-01781 

-01163 
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0.1731 
0 07 1 5 
0.1583 
0.2 7 9 0 
0.2 4 3 0 
0 0 0 3 1 
0.1366 
0.0 9 5 2 


FFF1 


-0.0177 

- 0.0010 


0 014? 
0.6174 
0 02 1 0 


0.0 3 8 8 
0 03 0 4 
0.0 2 3 6 
0 0102 
0.007 7 


STATION NO 


10.160 
12.700 
1 5 2 40 
17.710 


0034 -0078 
0 067 -0 032 
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0134 0022 


22 

160 

25 
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0.1 6 1 0.0 7 7 -0 *030 

0375 0.087 -0048 

0.459 0.089 -0082 

0 5 8 3 0 09 0 - 0 . 078 


-0.163 
- 0.202 
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-0 248 
-0 29 2 
-0.316 
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-0.346 
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-0.371 


0 

063 

0 
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0 
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0 
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0 
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0 
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0 
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0 7 6 40 
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0.043 
0.030-0 .0370 


0 017 1 
0 0 2 4 6 
C 0 2 7 5 
0.0 3 3 8 
0 0 4 2 8 
0.0 4 7 4 
0 .052? 
0.05 4 1 
0.0 4 7 2 
0 . 04 0 0 


0.0192 
0 02 1 7 
0 .0 20 5 
0.00 9 6 
0 00 6 3 
0 0 02 4 


-01652 
0 0 4 6 1 
0 3 3 1 7 
0 3 6 0 2 
0.5 3 5? 
0 6 3 1 1 
0.4 7 0 0 
0.0416 
-0.0410 
-0 2575 
-0.264c 
-0.2481 
-0.126? 
-0.0467 
- 0.0068 
-0.0024 
-0 . 02 30 


TABLE III REYNOLDS STRESS DATA MEASURED WITH LASER VELOCIMETER (CONT.) 



















































































OF POOR QUALITY 



y.cm 
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U'U, el 

0 4 35 

0 0 3 2 

0 111 

1 270 

0 0 6 4 

0.150 

1 705 

0 0 75 

0 160 

2 5 4 0 

0.127 

0 178 

3 175 

0.157 

0 2 08 

3.610 

0 17 1 

0 2 3 2 





STATION NO 
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0 

635 

0 

.039 

0 
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0 

077 

0 

367 
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705 
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.116 
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0 
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0 
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0 
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0 

4 1 6 
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0B0 

0 
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0 
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0 
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0 
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0 
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0 
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0 

6 18 

0 
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12 
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0 

773 

0 
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.240 
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727 

0 

7 17 

10 0 0 
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6 3 5 0 

7 620 

8 8 7 0 
10 14 0 
1 2 7 0 0 
13.770 


0.087 -0.5200 

0076 -0 5030 

0.054 -0.3300 

0.042 -0.1430 

-0.0320 
-00110 


0.0 3 2 8 
0 03 7 1 
0 . 0307 
0.026 1 
0.0166 
0.0125 
0 0035 
0.002 1 
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-0 4 6 55 
-0 64 1 1 
-04517 
-0.4101 
-0 . 26 7 2 
-0.1732 
-0 0577 
0 . 0066 
0 .00 00 ; 


0.0210 
0 02 4 4 
0 0 27 3 
0 .0 3 57 
0.0 35 7 
0.0 3 6 7 
0 026 2 
0 0 32 1 
0.0211 
0 0114 
0.004 6 
0 002 2 
0.0000 


TABLE III REYNOLDS STRESS DATA MEASURED WITH LASER VELOCIMETER (CONCLUDED) 
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TABLE IV PRESSURE TRAVERSE DATA (MEASURED TOTAL PRESSURES AND INFERRED STATIC PRESSURES) 
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TABLE IV PRESSURE TRAVERSE DATA (MEASURED TOTAL PRESSURES AND INFERRED STATIC PRESSURES) (CONT.) 
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TABLE IV PRESSURE TRAVERSE DATA (MEASURED TOTAL PRESSURES AND INFERRED STATIC PRESSURES) (CONT. 






Wr.1* ‘ U 



TABLE TV PRESSURE TRAVERSE DATA (MEASURED TOTAL PRESSURES AND INFERRED STATIC PRESSURES) (CONT. 
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TABLE IV PRESSURE TRAVERSE DATA (MEASURED TOTAL PRESSURES AND INFERRED STATIC PRESSURES) (CONT. 
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TABLE IV PRESSURE TRAVERSE DATA (MEASURED TOTAL PRESSURES AND INFERRED STATIC PRESSURES) (CONT. 
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TABLE IV PRESSURE TRAVERSE DATA (MEASURED TOTAL PRESSURES AND INFERRED STATIC PRESSURES) (CONT.) 
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TABLE IV PRESSURE TRAVERSE DATA (MEASURED TOTAL PRESSURES AND INFERRED STATIC PRESSURES) (CONCLUDED) 
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TABLE V HOT-FILM TPAVERSE DATA (Cont'd) 
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TABLE VI 

STREAMLINE COORDINATES DETERMINED FROM LASER VELOCIMETER MEAN 

VELOCITY DATA 



X LOCATION 


26.7 

31.8 

40.6 

48.3 

52.1 

55.9 

59.7 

67.3 

74.9 

82.6 

90.2 

97.8 

105.4 

113.0 

119.4 

140.3 

(cm) 

1 

10.5 

12.5 

16.0 

19.0 

20.5 

22.0 

23.5 

26.5 

29.5 

32.5 

35.5 

38.5 

41.5 

44.5 

47.0 

55.25 

(in.) 

1 

0.00 

0.00 

0.00 

3.10 

4.62 

6.99 

8.94 

14.35 

16.99 

15.95 

11.76 

4.37 

0.00 

0.00 

0.00 

0.00 


2 

1.27 

1.55 

2.77 

5.56 

7.14 

9.37 

11.20 

16.54 

19.08 

18.16 

14.45 

8.66 

3.84 

1.85 

1.32 

1.07 


3 

2.54 

2.95 

4.45 

7.49 

9.19 

11.43 

13.54 

18.82 

21.26 

20.45 

16.89 

11.51 

6.65 

3.86 

2.90 

2.26 

< 

4 

3.81 

4.32 

6.05 

9.25 

11.05 

13.46 

15.77 

21.03 

23.39 

22.61 

19.08 

14.07 

9.17 

5.89 

4.55 

3.61 

n c 

5 

5.08 

5.69 

7.65 

10.97 

12.95 

15.44 

17.93 

23.19 

25.53 

24.77 

21.18 

16.48 

11.51 

7.92 

6.22 

5.13 


6 

6.35 

7.06 

9.17 

12.70 

14.73 

17.35 

19.96 

25.30 

27.64 

26.92 

23.24 

18.75 

13.72 

9.83 

7.87 

6.68 

8 

7 

7.62 

8.41 

10.67 

14.35 

16.48 

19.23 

22.02 

27.46 

29.82 

29.08 

25.32 

20.90 

15.80 

11.58 

9.37 

8.23 

V 

8 

8.89 

9.73 

12.14 

15.98 

18.21 

21.13 

24.16 

29.72 

32.11 

31.32 

27.53 

23.04 

17.83 

13.31 

10.80 

9.73 


9 

10.16 

11.02 

13.59 

17.60 

19.94 



32.16 



29.87 

25.20 

19.84 

14.99 

12.17 

11.18 
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TABLE VII 

BACKFLOW BOUNDARY LAYER PARAMETERS 


STATION 

NO. 

x, cm 

N, cm 

U N /U ref 

NU 7 /p 

<y 

Be d 

C, Re e 

H 

25 

94.0 

0.097 

0.156 



0.0221 

38 

0.84 

1.98 

24 

90.2 

0.15 

0.191 

— 

0.0163 

57 

0.93 

1.90 

23 

86.4 

0.25 

0.210 


0.0149 

62 

0.92 

1.87 

22 

82.6 

0.50 

0.202 

— 

0.0090 

125 

1.13 

1.70 

21 

78.7 

0.63 

0.194 

60 

0.0124 

111 

1.38 

1.59 

20 

74.9 

0.50 

0.192 

125 

0.0105 

124 

1.30 

1.61 

19 

71.1 

0.51 

0.178 

110 

0.0081 

184 

1.49 

1.60 

18 

67.3 

0.76 

0.154 

155 

0.0076 

174 

1.32 

1.54 

17 

63.5 

0.76 

0.139 

147 

0.0074 

158 

1.17 

1.55 

16 

59.7 

0.76 

0.126 

132 

0.0059 

174 

1.03 

1.48 


PARAMETERS FOR 
TURBULENT BOUNDARY LAYER 
UNDER CONSTANT PRESSURE 
COLES (REF. 30) 

240 

0.00590 

425 

2.51 

1.535 

500 

0.00426 

1150 

4.88 

1.445 


mmsm 

2650 

8.94 

1.390 

2000 


5650 

16.36 
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APPENDIX A 

LASER VELOCIMETRY SYSTEM 


A.l System Description 

The LV system (fig. 3-8) consisted of a 2W argon-ion laser, a backscatter 
optical system employing a Bragg cell and a 3.73X beam expander, a TS1 1990B 
counter type signal processor, a PDP 11/10 computer for on-line data reduc- 
tion, and a 7.8 MB hard disk with integral 8 in floppy disk backup for data 
storage (and subsequent off-line data reduction). Traversing capability was 
achieved by mounting the laser and optical system on a lathe bed. All system 
components are commercially available. 

The LV was operated in a dual beam or "fringe" mode in which light from 
the intersection of two incident beams is heterodyned to detect the Doppler 
shift from an injected seed particle moving, at the local, instantaneous fluid 
velocity. In this mode, the LV measures the velocity component in the plane 
of the incident beams that is perpendicular to the bisector of the beams. The 
effective shape of the resultant measurement volume is an ellipsoid with major 
axis in the direction of the bisector of the beams. Sketch A (p. 21) shows 
these features and the theoretical measurement volume dimensions assuming that 
the ellipsoidal surface is defined by the locus of points where Doppler signal 
amplitude is 1/e 2 of its centerline value. No direct measurement of the 
effective measurement volume size was made; it can be effected by signal 
amplitude and the signal processor threshold level setting. To obtain maximum 
resolution, the major axis was aligned in the spanwise direction. Resolution 
of the LV system was quite high when compared to the relevant dimensions of 
the experiment. The minor axis of the measuring ellipsoid, d m , was 1/850 of 
the inlet boundary layer thickness and the major axis was 1/560 of the test 
section width . 


A. 2 Optical System 

The optical system, used a 3. 7 5X beam expander, a relatively wide angle 
lens (k = 4.9 deg), and 152 mm dia optics to produce a sufficiently high 
signal/noise ratio (SNR) to permit measurements to be made in the backscatter 
mode. The SNR was enhanced by a factor of 45 compared to an equivalent 
optical package using standard 50 mm dia optics without a beam expander. 

Focal length (762 mm) of the transmitting and collecting lens was set by the 
transverse dimensions of the test section. A 10 MHz Bragg shift was used to 
eliminate directional ambiguity and to provide 360 deg acceptance angle for 
all particle velocities. The 10 MHz Bragg shift was accomplished by 
upshifting the frequency of one incident beam by 40 MHz and downshifting the 
collected Doppler signal by 30 MHz. 
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A. 3 Signal Processor 


The TS1 1990B counter-type signal processor featured user selectable 
measurement modes, fringe count (N f = 2 N , N = 1 to 7) , validation accuracy 
(1% to 20% in 9 increments), input gain, and large amplitude signal rejection 
level (i.e. large particle discriminator). Four measurement modes were avail- 
able: (1) continuous mode (CONT) in which a measurement is made each time the 

minimum number of cycles is satisfied; (2) single measurement per burst mode 
(SM/B); (3) total burst count mode (TBC) in which the number of cycles in the 

burst are the Doppler frequency of the first Nf cycles are measured; (4) 
total burst mode (TBM) in which the duration of the burst and number of cycles 
in the burst are measured. The CONT mode was used for all measurements 
reported herein and the SM/B mode was used at selected points for redundancy. 
The CONT mode was chosen to minimize individual realization bias errors (see 
discussion below). Note also the CONT mode data is a superset of SM/B data 
because CONT data can be processed using the measured time between data points 
as a discriminator to eliminate multiply counted samples per burst to yield 
SM/B data. 


A. 4 Data Processing and Storage 

The time period and the time between samples for each CONT sample 
obtained in this test program has been stored and retained on floppy disks to 
permit future additional processing. 

After tr ial-and-error testing the 32 fringe count setting was found to 
yield an optimal combination of high data rate and "clean" histograms having 
minimal spurious data points in the tails. In this setting data validation is 
performed by making a 5/8 comparison between the time period for 20 cycles and 
the time period for 32 cycles. Validation accuracy was chosen to be 2% which 
permitted signals from seed particles following turbulence fluctuations up to 
approximately 25 kHz to be accepted. The signal gain and large amplitude 
signal rejection level were adjusted to reject 10 to 25% of the Doppler bursts 
to eliminate particle by biasing. 

Based upon the uncertainty analysis in Appendix C and the large 
turbulence levels generated by the separation bubble, 2000 samples were 
acquired per data point to generate velocity histograms for mean velocity and 
turbulence intensity measurements. For Reynolds stress measurements the 
previously measured local turbulence intensity was used to determine whether 
2000, 4000, or 8000 samples per data point were required to obtain the desired 
accuracy. Additional details concerning LV data processing are given in 
Appendix C. 
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A. 5 Seeding System 


Titanium dioxide powder having a nominal particle size of 1.0 pm was 
dispersed into the tunnel airstream to provide highly reflective light scat- 
tering centers for LV measurements using a fluidized bed seeding system driven 
by a low pressure dry air supply. A mixture of titanium dioxide powder and 
Cab-o-Sil, a silicon deaglomerating agent, were fluidized in the seeder to 
produce a densely seeded airstream which was directed to the tunnel. The seed 
was injected into the tunnel just downstream of the aftermost turbulence 
suppression screen using a cylindrical seeder probe, shown in figure A-l , 
which was designed to minimize the disturbance to the tunnel airstream. A 
similarly configured probe, designated the "zero-wake seeder" has been 
developed independently by Simpson (ref. 55). In principle, the flowrate 
through the seeder is adjusted until the momentum of the seeded air injected 
into the base region of the cylinder equals the cylinder drag and eliminates 
the wake deficit. Optimal seeder probe operation was achieved using a dry air 
supply pressure of 7 psi. Residual turbulence in the seeder wake was reduced 
by the 4 to 1 tunnel contraction. Total pressure probing revealed no 
discernible wake deficit at the test section inlet during seeder operation. 
Similarly, hot-wire measurement showed no difference in freestream turbulence 
level at the test section inlet during seeder operation compared to the clean 
tunnel operation with the seeder probe out of the tunnel. As expected, the 
turbulence level was increased when the probe was in the flow but not 
operating. 

The seeder probe produced a seed cloud having an approximately circular 
cross-section with a 12.7 cm dia at the test section inlet. During LV 
traversing upstream of separation the seeder probe was frequently repositioned 
to center the measuring volume within the seed cloud. Traverses downstream of 
separation required infrequent probe repositioning because of the rapid diffu- 
sion of the seed cloud within the separated flow. For off-centerline 
traverses the seeder probe was moved to an off-centerline location. 



Figure A-1 Zero-Wake Seeding Probe 



APPENDIX B 


LASER VELOCIMETER DATA REDUCTION EQUATIONS 

Figure B-l below shows the four laser beam orientations used to extract 
axial and transverse velocity data with a vertical traverse of the LV 
measurement volume. , and represent the magnitudes of the 

instantaneous velocity in the directions indicated by the sketch during a 
validated particle count. Veloc ities and U 2 were used to measure the axial 
and transverse velocity fields respectively. Ve loc i t ies and were used 
primarily to determine the Reynolds stress (uv) and secondarily to provide 
redundant axial and transverse mean velocity measurements. The angle £ was a 
small 5 deg laser beam inclination angle required to obtain transverse 
velocity and Reynolds stress measurements near the test surface. The angles 
cig and were set at 45 deg for Reynolds stress measurements. 

The desired axial and transverse velocities are related to the measured 
velocities by the following equations: 

/V/ ~ 


u x = u 

(Bl) 

/v /V 

U 2 “ Vcos£ + Wsin£ 

(B2) 

U 3 = Ucosa^ - Vsina^cos? - Wsing 

(B3) 

^ 'N* rsj rs* 

U 4 = Ucosa^ - Vsina^cos? + Wsin£ 

(B4) 

Taking each velocity, U, as the sum of a 
u(t), yields the mean velocity eqs . (B5) 

mean, U, and a fluctuating part, 
through (B 8 ) 

u = u x 

(B5) 

u 9 

V = * W tan 5 

cosg 

(B 6 ) 
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u = 


l^sirn^ 


- U/.sina': 


+ W 


sing (sinot 3 + sinct^) 


(B7) 


cosa^sina^ - sina^ cosa^ 


cosa oSina/, - sinctocosa 


V = 


U4COSCX3 - l^cosa^ 


tangfcosao + coscca) 

-W — 3 ^ (B8) 


cosgf cosa^sinctA - sina^cosa^) cosa^sina^ - siix^cosa/. 


The measurement of U and V are redundant because eqs. (B5) and (B6) are 
independent from eqs. (B7 ) and (B8). 

Evaluating eqs. (B7) and (B8) for a 3 = +45 deg and a 4 = -45 deg yields 
considerably simplified eqs. (B9) and (BIO). 

u = ( U 3 + U 4^ (B9) 

V = — — fUi - U,J - /T Wtan£ (BIO) 

2cos? J 4 


time averaging, and subtracting the mean flow equations yields 
for the axial and transverse turbulence intensities. 


(Bll) 


— 2 vw tan g - w 2 tan 2 g (B12) 

cos z g 


and the Reynolds stress 


uv = ( u , 2 - u^ 2 ) - / 2 uw tang (B13) 

2cosg 4 J 


after setting a 3 « +45 deg and = -45 deg. 


In evaluating equations (B6) and (BIO-13) only the first terms on the 
right hand side were used. The additional correction terms, which contained 


Squaring, 

equations 
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unmeasured spanwise velocity components, were very small because of the two- 
dimensionality of the flow (W - 0) and the small inclination angle (5 = 5 
deg). Away from the immediate vicinity of the wall (5 = 0) the correction 
terras disappeared. 
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Figure B-1 LV Measurement Geometry 


APPENDIX C 


LV ERROR ANALYSIS 


A procedure for quantifying LV measurement errors has been developed 
(ref. 37) in conjunction with the present study as part of UTRC’s independent 
research activities in gas dynamics. This uncertainty analysis includes the 
consideration of fixed (bias) and precision (random) errors and the methods 
for calculating the propagation of measurement errors through the system. 
Errors in the LV measurement system have been categorized as: (1) data 

processing errors, (2) laser beam geometrical errors, (3) processor errors, 
and (4) errors associated with seeding. Data processing errors arise from 
averaging a finite number of data samples per data point. Processor errors 
are the clock synchronization error, the quantizing error, the threshold limit 
error, the pedestal removal filter error, and the electronic noise induced 
errors. Laser beam geometrical errors include positioning uncertainty of the 
probe volume, angular sensitivity of the probe volume, fringe spacing uncer- 
tainty, and beam orientation errors, as well as, limitations imposed by a 
finite-sized probe volume. Seeding errors include flow distortion caused by 
seed injection, errors associated with the arrival rate of seed passing 
through the probe volume (individual realization bias), and particle lag 
errors in accelerating (or decelerating) flowfields. 

In this Appendix the error analysis presented in reference 37 will be 
summarized briefly. The error analysis methodology including the definition 
of terms, the determination of measurement uncertainty, and the propagation of 
measurement errors will be explained. Then, each error source will be 
described briefly. Finally, values of the precision and bias errors for the 
current experiment will be calculated and combined to determine estimates of 
the total uncertainty of LV measurements throughout the flowfield. 


C.l Error Analysis Methodology 


C.1.1 Definition of Terms 


An error is the difference between the measurement and the true value of 
the parameter being measured. Uncertainty is the maximum error which reason- 
ably might be expected in a parameter which is measured or which is computed 
from measured data. Measurement errors are of two types, random (or preci- 
sion) errors and fixed (or bias) errors. 

Random errors occur because of variations in repeated measurements of the 
same parameter. The precision index, S, is an estimate of the random (or 
precision) error. For N measurements (X-j^ , . .., X^ , ... X^) of the 

parameter X the precision index is defined as 
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s = 


N 

I (x r x) : 

i=L 

N - 1 


(C-l) 


where X, the mean of the measured values , is defined as 


X 


1 


N 


N Jj. 


(C-2) 


Bias is the fixed or systematic error. In repeated measurements, each 
measurement has the same bias, as illustrated in figure C-l. Examples of bias 
errors encountered in LV measurements are laser beam geometrical errors, 
process errors, and errors associated with seeding. 

C.1.2 Determination of Measurement Uncertainty 

The uncertainty is determined by combining the precision and bias errors 
in the following manner: 


U = +_ (B + t 95 S) 


(C-3) 


where B is the bias limit and tg^ is the 95th percentile point for the two- 
tailed Student's "t" distribution (ref. 56). When the bias limit is non- 
symmetric, upper (U + - B + + tg^ S) and lower (U = B - tg^ S) uncertainty 
limits are calculated. The tg^ value is a function of the number of degrees 
of freedom used in calculating S as shown in Table C-l. In a sample, the 
number of degrees of freedom is the size of the sample. For a statistic 
calculated from the sample the number of degrees of freedom is reduced by one 
for every estimated parameter used in calculating the statistic. For example, 
the mean, X, calculated from the sample (eq. (C-2)) has N degrees of freedom 
and the precision index, S, calculated from the mean (eq. .(C-l)) has N-l 
degrees of freedom. 

C.1.3 Propagation of Measurement Error 

Although only the uncertainty parameters, U, need be presented with the 
data, the bias and precision values are required to determine the propagation 
of the measurement error. The proper method for combining elemental measure- 
ment uncertainty values is to determine the root-sum-square values of the 
elemental bias limits and the elemental precision indices separately using a 
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“f 
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FREEDOM 

FREEDOM 

1 

12.706 

17 

2.110 

2 
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18 

2.101 

3 

3.182 

19 

2.093 

4 

2.776 

20 

2.086 

5 

2.572 

21 

2.080 

6 

2.447 

22 

2.074 

7 

2.365 

23 

2.069 

8 

9 

10 
11 
12 

13 

14 

15 

16 

2.306 

2.262 

2.228 

2.201 

2.179 

2.160 

2.145 

2.131 

2.120 

24 

25 

26 

27 

28 

29 

30 

2.064 

2.060 

2.056 

2.052 

2.048 

2.045 

2.042 
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Table C-1 Two-Tailed Student’s “t” Table 
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Taylor series expansion of the function being calculated from the elemental 
measurements. For example, the streamwise Reynolds stress, aB , is determined 
from measured values of the Reynolds stress, axial and transverse turbulence 
components, and the flow angle 


ctg = uv (cos <|> - 


o 

sin <(> ) 


(V 2 - 


u 2 ) 


sin<J> cos<j> 


(C-4) 


The precision index for the streamwise Reynolds stress, , can be determined 
as 


s °i6 - 


nUV 

2 

3 («B ) 2 

r* V 

2 

3 (aB ) 2 

oU 

2 

3 (aB ) ' 

- - pf 1 

2 ) 

b 

9 (uv) 

+ 

b 

_3 (v 2 ) 

T 

b 

_3 (u 2 ) 


b T 

3<}> 

1 


(C-5) 


Performing the indicated differentiation yields 



r 

2 

/ — \2 / — \2l 

s aB = 

(c.os 2 ()i - sin 2 <|>) S uv 

. 9 2 

+ sin <|> cos <J> 

[\s u2 j + \s v 7 J 


+ 





- 4uv sin<(> cos<J> 



(C-6) 


Similarly, the bias limit of the calculated streamwise Reynolds stress, 
B a ^ , is propagated from the bias limits of the measured variables . 
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3 (aB ) 


3 u 


3(ctB ) 


L 9<j> 


\1 / 2 


(C- 7) 


To propagate non symmetrical bias limits, the bias limit portion of the 
analysis must be completed for both the upper and the lower limits. The 
uncertainty in the calculated value of U is then determined by combining the 
bias and precision indices as shown in eq. (C-3) . 
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C.2 Data Processing (Precision) Errors 


Data processing errors are precision errors which arise from averaging a 
finite number of data samples per data point. In general, the precision error 
for any finite number of measurements of a fixed quantity can be calculated 
using eqs. (C-l) and (C-2). However, in LV velocity measurements, the 
velocity being measured does not remain constant during the sampling period 
but fluctuates due to turbulence. Thus, for LV measurements, equations (C-l ) 
and (C-2) can be written in the form 


S 



(u r u ) 2 


N-l 


(C-8) 


U 


- I 

N i=l 


U i 


(C-9 ) 


where is the velocity of the ith sample and U is the sample mean velocity 
which is an unbiased estimate of Up, the population mean velocity. The 
precision (or random) error calculated by eq. (08) is an estimate of the rms 
turbulence level of the flow being sampled. In the limit of an infinite 
number of samples, S becomes an exact measurement of the turbulence level. 

For a finite number of samples both the indicated rms turbulence level 
calculated in eq. (C-8) and the mean sample velocity calculated in eq. (C-9) 
will deviate from the true turbulence level and mean velocity of the flowfield 
by precision errors and respectively. 

C.2.1 Precision Errors in Turbulence Measurements 


Since the mean square turbulence velocity has a Chi-squared distribution 
(ref. 57) S ? can be obtained from tables by entering the number of degrees 
of freedom, N-l. For a large sample size, however, the inconvenient tables 
are unnecessary since 


s£ _ 1 

£ /2N 


N > 50 


(C-10) 


Values of the precision error in turbulence measurement are calculated for 
various sample sizes typical of LV measurements in Table C-2. In the present 
study 2000 samples were taken to acquire each turbulence measurement. Thus, 
the precision error in the measurement of the turbulence component is ± 1.58 

percent of the turbulence magnitude. 
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TABLE C-2. PRECISION ERROR IN TURBULENCE 
MEASUREMENT, S ? /e 


N 


500 

.0316 

1000 

.0224 

2000 

.0158 

4000 

.0112 

8000 

.0079 


C.2.2 Precision Errors in Mean Velocity Measurements 

Relative to mean velocity measurements, the sampling distribution of U is 
normal about U p as a mean with a standard deviation u p / /N. Thus, the preci- 
sion error in the mean velocity measurement can be estimated as 



(C-ll) 


which is a function of the turbulence level. Up and £ are unknown but can 

be approximated by the measured quantities U andV^?~/U. Values of the preci- 
sion error in mean velocity measurement for turbulence levels and sampling 
sizes typical of LV measurements are given in Table C-3. Note that both the 


TABLE C-3. PRECISION ERRORS IN MEAN VELOCITY 

o U 9 

MEASUREMENT, -jf — x 10 z 
U ref 



^ — x 100 
U ref 

N 

1% 

5% 

15% 

25% 

500 

.045 

.224 

.671 

1. 118 

1000 

.032 

.158 

.474 

.791 

2000 

.022 

.112 

.335 

.559 

4000 

.016 

.079 

.237 

.395 

8000 

.011 

.056 

.168 

.280 


218 



precision error and the turbulence level have been normalized by U re f to avoid 
extreme local turbulence levels which occur in regions of the flow where 
U « 0, Since 2000 samples were acquired to obtain each mean velocity 
measurement, the precision error in mean velocity varies from ± 0.02% in 

regions of low turbulence intensity ( V^"/ U re f = 0.01) to ± 0.56% in regions 
of high turbulence ( Vu^/U re f = 0.25) . 

The precision error in Reynolds stress measurements obtained with a one- 
component LV system can be shown to be (ref. 37) 

S? (C-12) 

The precision error in the streamwise shear stress component, a£ , can then be 
calculated using eq. (C-6). Calculated values of the precision error in the 
Reynolds stress and the streamwise shear stress measurements for each data 
point are presented in Table III. 


C .3 Bias Errors 


C.3.1 Laser Beam Geometrical Errors 


C . 3 . 1 . 1 ^ini^e^ro^e^ol\me_ B^a_s_(P\TB) 

Finite probe volume bias occurs because LV measurements represent a 
spatial average of the velocity throughout the probe volume rather than at a 
point in space. This bias error tends to increase the turbulence measurement 
and may increase or decrease the mean velocity measurement depending on the 
velocity gradients within the probe volume. 

C . 3 . 1 . 2 Be^am^Lo^c^tJ-Oji Bijas__ (BLB) 

Beam location bias is caused by the uncertainty in positioning the probe 
volume within the test section. It can be calculated as the product of the 
estimated velocity or turbulence gradients and the positional uncertainty in 
the direction of the gradient. 

C.3.1.3 Beam Orientation Bias (BOB) 


Two types of bias errors arise from beam orientation errors. The first 
type, B £ , which effects all measurements, is caused by misalignment of the 
laser beams relative to the rig axis. The second type, B^ 5 , is caused by the 
uncertainty in setting ± 45 degree beam orientation during shear stress 
measurements . 
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C . 3 . 1 . 4 Fr c_i_ng_Un c_er_t aj_n_t y_ (FSjB ) 

This bias is due to the uncertainty in the fringe spacing caused by the 
error, Ak, in setting the angle, k, between the intersecting laser beams. 

C . 3 . 1 . 5 N e gat^i v e_V ocj^t y_ Be am_( NVB ) 

The laser velocimeter measures the relative speed between the seed 
particles in the flow and the fringes in the measurement volume. If the 
fringe pattern is stationary the LV system cannot discriminate between seed 
particles traveling in the positive or negative directions. The LV processor 
thus rectifies the velocity signal in the same manner as a hot-wire. Conse- 
quently, in flow fields containing regions or periods of flow reversal the 
ensemble-averaged velocity will be biased high and the turbulence level will 
be biased low. When a Bragg shift is used to move the fringes in the negative 
direction at a velocity greater than the maximum negative velocity, NVB is 
eliminated. 


C .3 . 1 . 6 ^n^ompJ L e_te_Sj 1 griaJ 1 j3i_as_C[SJB ) 

Incomplete signal bias is caused by the angular sensitivity of the 
measurement volume. The probability that a seed particle passing through the 
measurement volume will cross the minimum number of fringes required to regis- 
ter a valid velocity sample decreases as the angle between the seed particle 
velocity vector and the normal to the fringe plane increases. Bragg shifting 
can be used, however, to achieve near isotropic response to eliminate ISB. 

C . 3 . 1 . 7 Fre_qu_en.cy_ Broa^enAn-g^B^as^ _(_FBB)_ 

Frequency broadening bias occurs when the laser beams do not intersect at 
the beam waists, where the wavefronts are plane, but elsewhere where the wave- 
fronts have a finite radius of curvature. The resulting measurement volume 
has spatial variations in the fringe spacing. In a properly designed and 
aligned LV system the beams will intersect at the beam waists and frequency 
broadening will be eliminated. 

C.3.2 Processor Bias Errors 


C . 3 . 2 . 1 BLa s €^ s _( C_A_B ) 

Comparison accuracy biases occur in flows in which individual particle 
velocity changes within the probe volume normal to the fringe pattern are 
significant. When changes in velocity occur within the probe volume, the 
comparator on the LV processor may reject valid signals causing a bias. Two 
types of comparison accuracy biases are particle acceleration bias (PAB) which 
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can occur in strong acceleration or deceleration fields and comparator toler- 
ance bias (CTB) which can occur in highly turbulent flows due to changes in 
flow direction within the probe volume and/or particle velocity changes 
induced by swirling turbulent eddies. 

C . 3 . 2 . 2 Cl_ock _Sync_hronJ.z_atj.on _Erro_r (CSB)_ 

The clock synchronization error is the bias caused by the mismatch 
between the randomly occurring Doppler burst and the start of the clock cycle 
in the processor. The uncertainty ranges from 0 to -1 clock pulse with all 
values in between being equally probable. 

C . 3 . 2 . 3 ^uant_i^ing_E£r o^r_( QBj^ 

The quantizing error is the bias arising from the determination of the 
frequency of the analog Doppler signal using the digital reference clock in 
the processor. The quantized time interval determined by the clock is 
accurate to within ±1 clock pulse with all values in between being equally 
probable . 

C . 3 . 2 . 4 Thr_es_hold_ L im i_t E rr or _(TLB_J 

The threshold limit error is the bias caused by the use of a non-zero 
volt Schmitt trigger in the processor to digitize the Doppler burst. The 
processor used in the current study used a true-zero-crossing detector which 
eliminated this bias. 

C . 3 . 2 . 5 EN B) 

The electronic noise induced error is the bias which occurs at low 
signal-to-noise ratios (SNR) when electronic noise mixed with the Doppler 
signal causes incorrect LV measurements. This bias can be significantly 
reduced by data validation circuitry within the processor and by operating at 
SNR > 3. 


C . 3 . 2 . 6 Pejd e_s t_a l__F_i l_t er Remov_a_l Erjror _(PFB_)_ 

This error is caused by the improper setting of the cutoff frequency of 
the high-pass pedestal removal filter. If the cutoff frequency is set too 
high, the Doppler burst caused by a low velocity particle will be distorted by 
the filter or even dropped out causing the velocity estimate from the LV 
processor to become biased high. If the cutoff frequency is at too low, some 
high velocity bursts will pass through the filter with residual pedestals 
causing a biased velocity estimate of the distorted frequency packet. 
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C.3.3 Seeding Bias Errors 


C . 3 . 3 . 1 F^low^D^s^or^tJ^ori B^ia_s_(FDB^) 

Flow distortion bias occurs when the mean flowfield is distorted by the 
seed injection process, Small seed particles, which have been sized to follow 
the flow, cannot be dispersed across the tunnel cross section except by 
natural diffusion. They cannot be projected across flow streamlines but 
instead must be convected with the seed injection jet. 

C . 3 . 3 . 2 P. a £_t_icJ L e_Lj_g_Bj 1 as^ ^PLB)_ 

The accuracy of LV flow velocity measurements is limited by the ability 
of the seed particles, which reflect the laser light, to move at the same 
velocity as the surrounding fluid. When the seed particles are too large to 
accurately follow rapid oscillations in the flowfield, particle lag bias 
occurs . 


C . 3 . 3 . 3 i^vi^ual_Re al^i_zat^i^n_B^as 

If the seed particles are distributed uniformly in space the number of 
particles passing through a unit volume in the flowfield varies directly with 
the particle velocity. High velocity regions of the flow contribute more LV 
samples per unit time than low velocity regions. In a shear layer or in 
turbulent flow the LV probe volune will span a region of the flow having both 
high and low flow velocities. In such cases the ensemble average velocity 
will be larger than the time-mean velocity of the flow because of the 
increased number of high velocity samples in the ensemble. A detailed 
discussion of this individual realization bias and correction schemes to 
reduce it are presented in reference 37. 

C . 3 . 3 . 4 Br ag£ Bia_s_(BB_)_ 

When the LV processor is operated in a continuous sampling mode while a 
Bragg shift is used to move the fringe pattern, slower particles will be 
multiply sampled more often than faster particles because they spend more time 
in the probe volume. This effect has been termed Bragg bias and it partially 
offsets individual realization bias. 


C.4 Determination of Total Uncertainty 

Before the uncertainty of LV measurements could be evaluated using the 
equations and graphs given in reference 37, the parameter values listed in 
Table C-4 were estimated. These parameter values included probe volume dimen- 
sions, processor settings, seed particle characteristics, estimates of the 
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TABLE C-4 

PARAMETER VALUES USED IN LV ERROR ANALYSIS 


FLOW PARAMETERS: 


Umox = 


AU 


Ax 


max 


27.4 m/sec , v max 
= lOO/sec, 


= 12.2 m/sec, w=0, <£ =30 deg 
, AU 

= 1920/sec , — = 0 

max 


av 

Ax 


max 


, AV 
= 200/sec, -gj 

max 


AV 

= 100/sec, ^ 


= 0 


AW _ AW _ AW 
Ax Ay AZ 


1 d -±A 

C ax 

i ac 

— x 

C ay 



55/m 

l 

2b 

= 25/m 
max 

1 


= 45/m 
max 

max 

I 

dx 

7~ 

d, 

max 

160/m 

1 

7 

ay 

= 32/m 
max 

i 

7" 

ay 

= l/m 
max 


i ?is. - 1 li i - L th 

i dz ~ C <*z " C dz 


BEAM GEOMETRY: M = 32, Nf = 44 , f B = 10. 1 6 MHz, SNR = 40 

d m = 133 /j.m, i m = i.6mm, d f = 3^m 

Ax = ± 2 mm , Ay = ± 025mm, Az = ±imm, a<£ = ± .25 deg 
fc = 4 .9 deg , Ak = ± .01 deg 


PROCESSOR PARAMETERS f c = 500 MHZ , COMP = 0.02, N = 2000 


SEED PARAMETERS: SEED MATER I AL , T i 0 2 : d p = l^m 
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maximum mean flow and turbulence gradients in each of the three coordinate 
directions, and uncertainties in laser beam positioning, orientation, and 
intersection angle. Each of the precision and bias errors discussed in the 
previous sections has been evaluated to determine estimates for the total 


uncertainty of LV measurements of U, 




The uncertainties for U 


and measurements have been estimated for regions of the flowfield having 

moderate ( V^/U ref - 0.05) and high ( V^/u ref = 0.20) axial turbulence 

levels. The uncertainties for V and VI? measurements have been estimated for 
regions of the flowfield having an intermediate transverse turbulence level 

= O.IO). Values for each error source and the total uncertainties 


( V? 


^ U ref 

are listed in Table C-5. 


When evaluating the uncertainties several biases were assumed to be 
negligible. They are the pedestal removal filter error (PFB), the threshold 
limit error (TLB), the finite probe volume bias (PVB), particle acceleration 
bias (PAB), negative velocity bias (NVB) , frequency broadening bias (FBB) , and 
flow distortion bias (FDB) . The negligible values for these biases follow 
respectively from the following assumptions: the pedestal removal filter has 

been set correctly, the processor has a true-z ero-c rossing detector, the probe 
volume is small compared to any spatial velocity gradients in the flowfield 
(PVB and PAB), Bragg shifting was used in regions of reversed flow, the laser 
beams intersected at the beam waists, and the zero-wake seeder probe caused 
negligible flowfield distortion at the measurement location. 


As shown in Table C-5, none of the processor bias errors exceeds 
percent. For mean velocity measurements beam geometry biases are also very 
low, the highest being the ± 0.7% bias due to uncertainty in probe volume 
location (BLB). The combined effect of seed bias errors on the mean velocity 
uncertainty is also very small. Particle lag bias, estimated from figure 8 in 
reference 37, is approximately 0.1 percent in both U and V components. The 
individual realization bias in the U component measurement would have been 

large in regions of high turbulence level ( V^7^ref = 0-20) but it 

was offset by Bragg bias because the processor was operated in the continuous 
sampling mode. Thus, the total bias from all sources is ± 0.0080 and ^0*0087 
for U and V measurements, respectively. When combined with the precision 
error estimated for 2000 sample measurements, the uncertainty in LV measure- 
ments of U is ± 0,010 and ± 0.017 in regions of moderate ( V?7u ref = 0.05) 
and high ( \/l?/U ref = 0.20) turbulence, respectively. The estimated 


uncertainty in the transverse velocity is -q re §i° ns intermediate 

transverse turbulence (V?/ U ref = °- 10 ) . This uncertainty value means that 
19 out of every 20 mean velocity measurements will be within these quoted 
accuracies as noted above in Section C.1.2. 
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TABLE C-5 

ESTIMATED BIAS ERRORS, PRECISION ERRORS, AND UNCERTAINTIES OF LV MEASUREMENTS 
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±0.0169 ±0.0102 + 0.0162 +0.087 +0.087 +0.068 

TOTAL UNCERTAINTY (U) - 0.0119 -0.106 -0.087 -0.059 


For turbulence measurements, the oias errors are expressed as a fraction 
of the maximum turbulence intensity (^ max = 0.25 U ref ) . Beam geometry biases 
are small except for the beam location bias which is estimated to be ± .055 
and ± .025 for measurements of axial and transverse turbulence components, 
respectively. This is a worst-case estimate, assuming the maximum positioning 
uncertainty occurs at the location of the maximum turbulence gradient. The 
effect of particle lag on turbulence measurements is negligible since the 1.0 
ym dia. TiC> 2 seed particles can track at 2500 Hz sinusoidal oscillation to an 
accuracy better than 1% and most of the energy in the turbulence spectrum was 
in the 0-100 Hz band (fi^. 4-15). The individual realization bias is esti- 
mated to be as high as nsn f° r axial turbulence component measurements in 

the high turbulence ( Vu^*/U re f = 0.20) environment but throughout most of the 
flowfield IRB is less than 1%. It is assumed that Bragg bias did not effect 
the turbulence measurements. 

When combined with the precision error (± 0.0158) estimated for a 2000 
sample measurement, the uncertainty in the measurement of the axial turbulence 
component is ± 0.087 and ^q"i 06 re gi° ns of moderate and high turbulence 
levels, respectively. The estimated uncertainty in the transverse turbulence 
component measurement is ^(j’o59 *' n re g*- ons of intermediate transverse 
turbulence level. As noted above, these uncertainty estimates are expressed 
as fractions of the maximum turbulence intensity (r = 0.25 U raf ) and are 
worst-case estimates. 
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MEASUREMENT 



Fig. C-1 Measurement Uncertainty, Nonsymmetrical Bias 
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APPENDIX D 


ERROR ESTIMATES FOR TOTAL PRESSURE AND HOT-FILM MEASUREMENTS 


D.l Accuracy of Total Pressure Measurements 

The accuracy of the total pressure measurements is dependent on several 
factors including the accuracy of the pressure transducers, the angular 
sensitivity of the pitot probes, and the effect of turbulence on the pitot 
reading* Each of these factors will be discussed below. 

D.1.1 Transducer Calibrations 


Pressure transducers were calibrated to an accuracy oft 0.0005 Q re f with 
a micromanometer (fig. 3-9) at the beginning of each test day. Transducer 
output voltages were zeroed prior to each traverse to compensate for thermal 
drift. Zeroes were also checked after each traverse to ensure that drift 
during the traverse did not exceed ± 0.002 Q re f otherwise the data was 
retaken. 

D.1.2 Angular Sensitivity of Probes 

The total pressure probes were calibrated for angular sensitivity in 
pitch and yaw over a ± 30 deg. range at three dynamic pressures: Q re f> 0.45 
Q re f, and 0.09 Q re f* The boundary layer probes were accurate to within 1% of 
the dynamic head pressure over ± 7 deg. relative to the normal to the probe 
shaft. The aspirated kielhead probes were accurate to within 0.005 Q re f over 
± 26 deg. relative to the normal to the face of the probe. The angular 
sensitivity of all probes used in the current study were essentially 
uneffected by flow speed over the range tested. 

Since the boundary layer probes were only used to obtain data near the 
test surface (y < 2 cm) where the mean flow inclination angle was within the 
± 7 deg. acceptance angle range of the probe, the angular sensitivity bias 
error in boundary layer pitot measurements is less than to.01 Qref exce P t near 
separation and reattachment. In those regions, the mean flow inclination 
angles exceeded the probe acceptance angle and frequent flow reversals added 
to measurement uncertainty. The angular sensitivity bias error in kielhead 
measurements is better than *q qq^ Q re j throughout the flowfield because of 
the large acceptance angle of the aspirated kielhead probe and because the 
probe was bent when necessary prior to traversing to optimize its alignment 
with the mean flow direction known beforehand from smoke flow visualization 
(see figure 4-6). 


receding page blank not filmed 
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D.1.3 Effect of Turbulence on Pitot Readings 


The instantaneous total pressure at a point in the flowfield can be 
expressed as 



(U 2 - 


W 2 ) 


(D-l) 


where l), V, and are the magnitudes of the instantaneous velocities in the 
axial, transverse, and spanwise directions, respectively. These instantaneous 
velocities can be assumed to be comprised of a time-mean component plus a 
sinusoidally oscillating unsteady component 


U = U + 2u sin oat (D-2a) 
V = V + 2v sin nit (D-2b) 
W = W + 2w sin oat (D-2c) 


Substituting eqs . (D-2) into eq. (D-l) and time averaging yields 


P = p + - p (U 2 + V 2 + W 2 + u 2 + v 2 + w 2 ) 
2 


(D-3) 


Thus, when a pitot tube is used in a turbulent stream its reading can 
exceed the total pressure corresponding to the mean flow by an amount propor- 
tional to the mean kinetic pressure of the turbulent velocity fluctuations, 
1/2 p (u z + + w^). Goldstein (ref. 58) suggested that the pitot reading 

should be corrected by the entire turbulent kinetic pressure. However, the 
magnitude of any correction should depend on the scale as well as the inten- 
sity of the turbulence fluctuations (ref. 59). 

In the current study, the turbulence scale is large compared to the 
diameter of the pitot probes employed. The flow approaches the probe at 
randomly varying inclination angles due to turbulence fluctuations. Such 
inclination angles can be quite large for turbulence intensities at which the 
mean kinetic pressure, a P t , is significant as shown in Table D-l. For 
example, at a 7.5 percent turbulence level having a mean kinetic pressure of 
0.0169 Q re f , the flow angle occasionally exceeds the ±30 deg. acceptance of 
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TABLE D-l. MAXIMUM FLOW ANGLES INDUCED BY TURBULENT FLUCTUATIONS 


(u = .5 u f , v = w = o, Vu 2 ” = vv 3 " = w") 



AP t 

max 

(deg) 

U ref 

Q ref 

0 

0 

0 

.025 

.0019 

±8.9 

.050 

.0075 

±19.5 

.075 

.0169 

±31.2 

.100 

.030 

±43.3 

.125 

.047 

±54.7 


the aspirated kielhead probes. At such large inclination angles the pitot 
probe would read low causing the time- aver aged pitot reading to exceed the 
total pressure corresponding to the static pressure and the mean flow by an 
amount which is less than the mean kinetic pressure of the turbulent velocity 
fluctuations . 

Since a rigorous procedure is not presently available, the total pressure 
data tabulated in Table IV have not been corrected for the effect of the 
turbulent kinetic pressure. However, it is the author's opinion that such a 
correction would be much less than the magnitude of the mean kinetic 
pressure . 

D.1.4 Estimate of Total Pressure Accuracy 

The accuracy of the total pressure measurement varied throughout the 
flowfield. In the freestream and near the test surface (except near separa- 
tion and reattachment) overall accuracy was better than ± 0.01 Q re f In other 
parts of the flowfield measurements were less accurate but high turbulence 
levels and flow reversals preclude a quantitative estimate of the error. 


D.2 Accuracy of Hot-Film Measurements 

In many flow experiments the accuracy of hot-film measurements can be 
easily quantified and is limited only by the accuracy of the calibration, the 
linearizer settings, and the sensitivity of the sensor to moderate flow 
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angularity. In the current study flow reversals, high turbulence levels, and 
extreme flow angularity combine to limit the accuracy of the hot-film measure- 
ments. Each of these factors will be discussed in this section. 

D.2.1 Probe Calibration and Linearization 


Hot-film probes were calibrated in a 3.8 cm dia low turbulence calibra- 
tion jet. The probes were aligned such that the sensor was perpendicular to 
the jet axis and the tips of the needle supports were parallel to the jet 
axis. Mean velocity and bridge output voltage were recorded for 20 jet speeds 
ranging from 0.01 Q re f to 1.20 Q re f The mean response equation for each 
sensor was assumed to conform to a modified King's law in the following form 


E 2 = A + B U?g£ 5 (D-4) 

where E is the bridge output voltage. 

The calibration data was then used to calculate coefficients for a 
fourth-order polynomial linearizer. Maximum deviation of the calibration data 
from the linearized fit was 2.5%. To account for hot-film aging and day-to- 
day tunnel temperature changes, a single point in-place calibration was con- 
ducted at the start of each test shift to ensure that a 10 volt linearizer 
output corresponded to 100 ft/sec (30.48 m/sec). This was accomplished by 
positioning the probe within the test section at a specific location (x = 7.6 
cm, y = 10.2 cm, z = 0) where the streamwise velocity was known to equal 
0.9925 U re £ and adjusting the linearizer span accordingly. 

D.2.2 Effect of Turbulence and Flow Angle Variations on Hot-Film 
Measurements 


Jorgensen (ref. 60) has shown that the effective cooling velocity of a 
hot-film/wire sensor can be expressed as 


U eff = U N 2 + h2 U B 2 + fc2 U T 2 (D-5) 

where TJg, and are the normal, binormal, and tangential components of 
the instantaneous velocity vector as shown in the sketch. These velocity 
components correspond to the instantaneous axial, T) , transverse, 'V , and span- 
wise, \5, velocity components in the test section. Thus, the effective cooling 
velocity can be written as 
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(D— 6 ) 


~ 9 ~ 9 ~9 2~? 

U eff = u 2 + h v + k w 



The coefficients h and k are functions of the probe geometry and should be 
determined from calibration. Typical values for h and k are 1.05 and 0.18, 
respectively (ref. 60). Replacing the instantaneous velocity components in 
eq. (D-6) with the time-mean and fluctuating velocity components (eqs. D-2) 
and time averaging yields the following expression for the effective cooling 
velocity 


U 


eff 


-if 


(U + u) 2 + h 2 (V + v) 2 ] 


1/2 


(D-7) 


assuming k2 is small and W ^ 0. After expanding the rhs of eqs. (D-7) with a 
binomial expansion and performing the indicated time average, the bias in 
single-sensor hot-film measurements of axial velocity due to transverse 
velocity and axial and transverse turbulence can be estimated. 


j U eff . , , h 2 (V 2 +v 2 > 
U 2 U 2 


(D-8) 


A comparison of hot-film and LV data acquired at Station 18 (x = 67.3 cm) 
is shown in figure D-l to indicate the magnitude of this bias error in the 
present study. Data at Station 18 was chosen for illustration because it is 
the location of the maximum bias error. Except near the wall, the difference 
between the hot-film data and U/U re f measured with the LV is approximately 0.1 
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U f. The hot-film data also measures approximately 0.05 U re £ higher than the 
streamwise velocity component , + V^/U^f. However, the profile 

calculated from LV data using eq. (D-7) approximates the hot-film data, thus 
verifying the magnitude of the bias error given in eq. (D-8). 

D.2.3 Effect of F low Reversals on Hot-Film Output 

Since a single-component hot-film probe cannot distinguish the direction 
of the effective cooling velocity, the bridge output voltage is always posi- 
tive, even in regions or periods of reversed flow. This rectification of the 
velocity signal can cause the measured time- averaged effective cooling 
velocity to be significantly larger than the streamwise mean velocity and the 
indicated turbulence level to be lower than actual when flow reversals occur 
during the averaging period (see Section D.3 below.) 


D.3 Comparison of Hot-Film, LV, and Pitot Data Near Separation 

LV and hot- film measurements of the mean streamwise velocity and turbu- 
lence profiles at Station 11 (x = 40.6 cm) near separation are plotted in 
figure D-2 . In addition, the mean streamwise velocity profile calculated from 
total pressure measurements assuming constant static pressures across the 
traverse is plotted for comparison. Several aspects of these profiles are 
interesting. 

Near the wall, where the forward flow fraction, Yp^ , is approximately 

0.5, the mean velocity measured with the LV approaches zero as y approaches 
zero. However, the pitot and hot-film data approach a constant positive 
value. This positive value occurs because the hot-film probe and, to a lesser 
extent, the pitot probe are directionally insensitive and tend to rectify the 
velocity fluctuations causing the t ime- averaged velocity to have a positive 
bias. Signal rectification also effects the turbulence measurements made with 

the hot-film in regions of the flow where Yr> < 1.0. Near the wall, the 

F U 

turbulence measured with the hot-film is half the LV measurement because the 
peak-to-peak oscillation sensed by the hot-film are half the actual value 
because of signal rectification. Near the freestream, where the turbulence 
levels are low, LV and hot-film mean velocity and turbulence measurements are 
nearly identical. The mean velocity inferred from pitot readings in signifi- 
cantly lower than actual, however, because of the negative transverse 
pressure gradient near separation. 
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DATA OBTAINED AT STATION 1 8, x = 67.3 cm, z = 0 
□ HOT-FILM DATA 



- 0.2 0 0.2 0.4 0.6 0 . 


^eff/Uref 

Figure D-1 Comparison of Hot-Film and LV Data 


85 — 6 — 36 — 10 





3. Recipient’s Catalog No. 

5. Report Date 

March 1987 

6. Performing Organization Code 


8. Performing Organization Report No. 

None (E-3254) 

10. Work Unit No. 

535-05-01 

11. Contract or Grant No. 

NAS3-22770 

13. Type of Report and Period Covered 

Contractor Report 
Final 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Project Manager, Thomas F. Gelder, Propulsion Systems Division. 


16. Abstract 

The separation and reattachment of a large-scale, two-dimensional turbulent 
boundary layer at low subsonic speed on a flat plate has been studied experimen- 
tally. The separation bubble was 55 cm long and had a maximum bubble thickness, 
measured to the height of the mean dividing streamline, of 17 cm, which was twice 
the thickness of the inlet boundary layer. A combination of laser velocimetry, 
hot-wire anemometry, pneumatic probing techniques, and flow visualization were 
used as diagnostics. Principal findings were that an outer inviscid rotational 
flow was defined which essentially convected over the blockage associated with 
the Inner, viscously dominated bubble recirculation region. A strong backflow 
region in which the flow moved upstream 100 percent of the time was measured near 
the test surface over the central 35 percent of the bubble. A laminar backflow 
boundary layer having pseudo-turbulent characteristics including a log-linear 
velocity profile was generated under the highly turbulent backflow. Velocity 
profile shapes in the reversed flow region matched a previously developed univer- 
sal backflow profile at the upstream edge of the separation region but not in the 
steady backflow region downstream. A smoke flow visualization movie and hot-film 
measurements revealed low frequency nonperiodic flapping at reattachment. How- 
ever, forward flow fraction data at reattachment and mean velocity profiles in 
the redeveloping boundary layer downstream of reattachment correlated with 
backward-facing step data when the axial dimension was scaled by the distance 
from the maximum bubble thickness to reattachment. 

17. Key Words (Suggested by Author(s)) 18. Distribution Statement " 

Separation bubble; Turbulent separation. Unclassified - unlimited 
Reattachment; LV error analysis STAR Category 34 


19. Security Classif. (of this report) 20. Security Classif. (of this page) 1 21. No. of pages 1 22. Price - 

Unclassified Unclassified 245 All 


1. Report No. 2. Government Accession No. 

NASA CR-4052 

4. Title and Subtitle 

Flowfield Measurements in a Separated and Reattached 
Flat Plate Turbulent Boundary Layer 

7. Author(s) 

William P. Patrick 


9. Performing Organization Name and Address 

United Technologies Research Center 
Experimental Gas Dynamics Group 
East Hartford, Connecticut 06108 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


i 


For sale by the National Technical information Service, Springfield, Virginia 22161 



